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ABSTRACT 

A  water  flow  visualization  facility  was  constructed  in  order  to  investigate  the  miving 
processes  inside  of  four-inlet  ducted  rocket  engines  (DRE)  with  various  flow  rates  and 
geometries.  The  observ'ed  flow  fields  were  documented  with  video  and  conventional 

I 

photography  so  that  relationships  between  various  geometric/flow  parameters  and  the  flow  field 
behavior  could  be  obtained.  Flow  visualization  was  achieved  by  seeding  the  flow  field  with 
polystyrene  spheres  and  illuminating  them  with  laser  sheet  lighting.  Images  revealed  complex 
mixing  patterns  within  the  dome  head  region  of  the  combustors  which  exhibited  a  strong 
dependence  on  dome  height,  momentum  ratios,  and  inlet  flow  angles.  For  large  inlet  flow  angles 
on  the  same-station  combustors,  it  was  shown  that  the  geometries  intrinsically  direct  a  significant 
portion  of  the  inlet  flow  towards  the  dome  head  region  without  the  need  for  diverter  devices. 
The  resulting  dome-region  flow  field  contains  two  distinct  recirculation  patterns  shown  to  be 
dependent  on  the  momentum  ratio  and  inlet  flow  inlet  angle.  These  parameters  also  strongly 
affect  the  existence  and  strength  of  transport  mechanisms  between  the  fore  and  aft  combustor 
segments.  By  designing  a  combustor  accordingly,  it  may  be  possible  to  take  advantage  of  one  or 
both  of  the  operating  conditions  resulting  in  a  combustor  with  improved  combustion  properties. 
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INTRODUCTION 

Ducted  rocket  engine  systems  have  recently  received  a  revived  interest  due  to  their  high 
specific  impulse,  throttlability.  and  suitability  for  long  range  applications.  They  obtain  the  high 
specific  impulse  by  ducting  the  inlet  flow  in  which  they  fly  through  inlet  ducts  and  using  it  as  the 
oxidizer  which  is  mixed  with  a  fuel  rich  gas  Irom  a  gas  generator.  The  mixture  combusts  and  is 
then  accelerated  through  a  converging/diverging  nozzle  to  generate  thrust.  A  diagram  showing 
the  general  layout  of  a  DRE  is  shown  in  Figure  1.  Much  of  the  flow  visualization  woiic 
performed  in  the  past  has  been  with  one  and  two-inlet  geometries.  These  geometries  offer  an 
aerodynamic  advantage  since  the  inlet  ducts  may  be  tailored  to  generate  lift  for  the  vehicle  but 
often  require  some  type  of  inlet  flow  deflecting  device  to  be  installed  in  the  ducts  so  that 
sufficient  inlet  flow  enters  the  dome  region  of  the  combustor  where  the  flameholding  process 
typically  occurs.  These  combustors  have  often  been  investigated  using  flow  visualization  with 
water  for  qualitative  information  and  direct-connect  hot-fire  testing  [Schadow  and  Chieze,  1981] 
for  quantitative  results. 

The  work  described  in  this  article  has  been  performed  on  a  facility  constructed  at  The 
University  of  Alabama  in  Huntsville.  This  facility  was  used  to  study  different  mixing  scenarios 
for  four-inlet  ducted  rocket  engine  geometries  tlirough  flow  visualization  with  water  flows 
simulating  both  the  “‘fuel  ’  flow  and  the  “air  inlet  flow’  flow.  Four-inlet  geometries  ofiBsr 
advantages  over  the  one  and  two-inlet  geometries  such  as  steer-to-tum  capability  and  the 
tendency  to  passively  divert  necessary  inlet  flow  into  the  dome  region  of  the  combustors  for 


nameholding  requirements  without  the  use  of  diverter  devices  commonly  needed  on  other 
combustors.  The  effects  of  inlet  flow  angle,  dump  station,  air/fuel  mass  ratios,  air/fuel 
momentum  ratios,  dome  height,  and  injector  performance  were  investigated  over  a  wide 
operational  range.  The  wide  range  allowed  the  investigation  of  flame  holding  and  performance 
potential  of  the  various  geometries  over  a  large  number  of  operating  conditions. 

FACILITY  DESCRIPTION 

The  flow  visualization  facility  (Figure  2)  has  overall  dimensions  of  6. 1  meters  by  4.6  meters. 
Not  shown  in  the  figure  are  the  3800  liter  reservoir  tank  and  the  centrifugal  pump  located  below 
the  deck.  The  power  plant  for  the  facility  is  a  22.4  kW  variable  RPM  centrifugal  pump  capable 
of  delivering  3000  liters  per  minute  at  207  kPa  to  the  combustor  model.  This  pump  was 
dedicated  to  delivering  high  volume  at  low  pressure  and  was  used  to  supply  water  as  the  "air" 
source  for  the  combustor.  A  supplemental  517  kPa  water  supply  arm  was  used  to  deliver  the 
water  for  the  "fuel"  source  at  up  to  132  liters  per  minute.  The  amount  of  water  delivered  to  the 
model  was  controlled  both  by  the  pump  RPM  and  by  throttling  the  supply  arms  with  globe 
valves.  A  PC  monitored  Blancett  401  flowmeters  used  on  each  arm,  two  pressure  transducers, 
and  provided  remote  pump  control  if  desired.  The  combustor  models  were  placed  inside  a  0.92m 
X  1 .00m  X  0.60m  viewing  tank  which  contained  water  to  minimize  refraction  effects. 

The  dome  head  (Figure  3)  is  a  concave  piston  head  with  0-rings  that  contains  the  fuel  injector 
and  which  allows  positioning  at  any  axial  location  in  the  front  end  of  the  combustors.  The  dome 
height  is  defined  as  the  axial  distance  between  the  beginning  of  the  inlet  flow  dump  station  and 
the  perimeter  of  the  dome  head.  A  five  orifice  injector  was  used,  one  axial  orifice  and  four  off- 
axis  orifices  angle  at  45  degrees  to  the  combustor  axis.  See  figure  3  for  injector  dimensions. 


Visualization  of  the  flow  field  was  achieved  by  seeding  the  reservoir  tank  with  0.0254  cm 
(0.010")  polystyrene  spheres.  The  fuel  jet  was  left  unseeded.  -‘Outgassing"  of  air  in  the  “fuel” 
supply  after  experiencing  the  pressure  drop  across  the  injector  resulted  in  a  mist  of  tiny  inlet  flow 
bubbles  which  provided  excellent  visualization  ot  the  fuel  flow.  .A  5  watt  Argon  Ion  laser  was 
used  to  illuminate  the  flow  field.  The  laser  output  was  first  focused  and  then  expanded  through  a 
cylindrical  lens  to  produce  a  light  sheet  approximately  1.0  mm  thick.  The  sheet  was  capable  of 
being  oriented  either  vertically  or  horizontally.  The  images  observed  were  then  captured  on  Hi8 
video  and  by  conventional  photography.  Figure  4  is  a  sample  photograph  of  a  45  degree  same- 
station  combustor  test  and  is  the  flow  field  depicted  in  Figure  12. 

RESULTS  AND  DISCUSSION 

Six  combustors  have  been  evaluated  (Figure  5).  Three  of  the  combustors  had  the  inlet  flow 
inlets  at  the  sarne  axial  station.  The  remaining  three  had  one  pinlet  flow  of  inlet  flow  inlets 
staggered  1.5  inlet  diameters  downstream  of  the  upstream  pinlet  flow.  The  three  six-inch  I.D. 
combustors  of  each  group  consisted  of  45.  60,  and  90  degree  inlet  flow  inlet  arm  dump  angles 
each  with  an  internal  diameter  of  3.81  cm.  All  were  investigated  up  to  a  Reynolds  number  of  5.5 
X  10  ,  but  no  Reynolds  number  effects  were  observed  above  Reynolds  numbers  of  1.5  x  10  .  The 
mass  flow  ratio  (rh,  /rhf)  and  momentum  ratio  (rii^v,  /ih^Vf )  were  varied  between  10  to  60 
and  0.43  to  2.58  respectively. 

Dome  Height 

Dome  height  was  the  first  parameter  investigated.  The  first  combustor  tested  was  the  90 
degree  same-station  combustor.  A  dome  height  ot  0.25D  (3.81  cm)  resulted  in  symmetric,  stable 


recirculation  regions  tor  the  90°  same-station  combustor  over  a  wide  range  of  momentum  ratios 
and  appeared  to  be  the  preferred  height  for  that  combustor  and  the  fuel  injector  used  based  on  the 
size  and  strength  of  the  recirculation  zones.  The  preferred  height  was  determined  after  testing 
each  dome  height  setting  over  the  range  of  momentum  ratios  specified  earlier.  Figure  6  shows 
the  operating  zones  for  momentum  ratio  versus  dome  height.  The  preferred  dome  height  was 
then  chosen  as  the  height  which  minimized  the  fluctuating  behavior  between  the  impinging  and 
passive  recirculation  regimes  (to  be  discussed  next).  Dome  heights  greater  than  0.27D  resulted 
in  the  fuel  jets  impinged  on  the  combustor  wall  and  not  the  incoming  inlet  flow  columns.  The 
preterred  dome  height  tor  the  60°  same-station  combustor  was  determined  to  be  «0.083D  (1.27 
cm)  by  the  same  manner.  An  acceptable  dome  height  for  the  45°  combustor  was  not  found,  so  a 
dome  height  of  0  cm  was  used  in  order  to  provide  the  strongest  recirculation  possible.  The 
alternating-station  combustors  did  have  problems  with  two  of  the  fuel  jets  impinging  on  the  wall 
because  the  same  fuel  injector  was  used  on  all  models.  This  situation  is  discussed  further  in  the 
alternating-combustor  section. 

Same-Station  Combustors 

The  90  degree  same-station  combustor  was  investigated  first  with  a  low  “fuel”  flow  rate.  The 
flow  field  contained  symmetric  counter-rotating  zones  on  both  sides  of  the  inlet  flow  inlet  station 
for  plane  A-A.  The  recirculation  zones  increased  in  strength  as  the  “inlet  flow”  flow  rate  increase, 
but  the  downstream  zone  appeared  to  fluctuate  in  size  and  shape.  Figure  7  shows  the  orientation  of 
the  sheet  lighting  relative  to  the  combustor.  Figure  8  shows  the  pairing  of  two  counter-rotating 
regions  found  between  each  of  the  inlet  arms  for  plane  B-B.  These  regions  appeared  to  be  tube¬ 
like  and  were  seen  to  contain  substantial  amounts  of  “fuel”.  The  endpoints  of  the  tube  could  not  be 


isolated,  but  appeared  to  begin  in  the  dome  head  region  and  terminate  approximately  1.0  inlet 
diameter  downstream  from  the  inlet  flow  inlet  station.  1  he  structures  seen  in  the  flow  field  did  not 
appear  to  develop  further  for  Reynolds  numbers  above  5.5  x  10^,  but  fluctuations  in  size  and  sh^ie 
were  observed. 

The  low  “fuel”  flow  rate  resulted  in  a  momentum  ratio  of  2.39.  For  the  high  momentum  ratios, 
the  fuel  jet  was  seen  to  “bend  back”  and  follow  the  governing  “inlet  flow  flow  (Figure  9).  This 
was  termed  a  “passive”  flameholding  scenario.  There  was  no  significant  change  in  this  flow  field 
until  a  momentum  ratio  between  1 .90  and  1 .6 1  was  reached.  In  this  range,  the  fuel  jet(s)  fluctuated 
between  the  “passive”  scenario  and  a  straight  fuel  jet,  termed  “impinging  (Figure  10),  which 
impinged  on  the  incoming  “inlet  flow”  column.  Below  a  momentum  ratio  of  1.61,  the  fuel  jet 
would  partially  penetrate  the  inlet  flow  column  and  then  bifurcate  into  two  branches.  One  would 
be  deflected  towards  the  wall  and  the  other  towards  the  stagnation  point  of  the  four  jets.  The  result 
was  two  recirculation  zones  split  and  driven  by  the  fuel  jet. 

The  60  degree  combustor  was  investigated  at  a  dome  height  setting  of  0.083D  and  over  the  same 
momentum  ratio  range.  The  recirculation  regimes  observed  near  the  dome  head  on  the  90  degree 
combustor  also  occurred  in  the  60  degree  combustor  but  were  seen  to  transition  at  a  higher 
momenmm  ratio.  The  “passive”  regime  (Figure  1 1)  existed  down  to  a  momentum  ratio  of  2.05, 
fluctuated  in  the  range  between  2.05  and  1.72,  and  became  “impinging”  below  1.72  (Figure  12). 
The  critical  momentum  ratio  values  are  defined  as  the  value  at  which  the  fuel  jet  transitions  fitnn 
passive  to  impinging.  This  wtis  taken  as  the  center  value  in  the  fluctuating  region.  The  critical 
momentum  ratios  for  the  60  and  90  degree  combustors  are  shown  to  be  similar  when  the  only 
lateral  momentum  component  of  the  “inlet  flows”  were  compared. 


The  downstream  recirculation  regions  and  reattachment  lengths  remained  relatively  constant  over 
the  momentum  ratio  range  investigated,  but  did  reduce  in  size  once  the  center  fuel  jet  penetrated 
the  stagnation  point  of  the  four  inlet  flow  columns.  The  transport  "tubes"  existing  in  between  the 
"inlet  flow"  inlet  arms  appeared  to  be  weaker  and  shorter  in  length  than  observed  in  the  90  degree 
combustor.  No  fuel  accumulation  in  the  dome  head  region  was  observed  for  this  combustor, 
indicating  that  the  transport  mechanisms  were  indeed  functioning. 

The  45  degree  same-station  combustor  exhibited  the  weakest  recirculation  zones.  This  was 
likely  due  to  the  reduced  amount  of  inlet  flow  momentum  normal  to  the  combustor  axis  which  was 
seen  to  produce  much  of  the  mixing  and  transpon  that  occurs  between  the  "fuel"  and  "air  inlet 
flow"  jets.  The  center  fuel  jet  was  seen  to  penetrate  the  "inlet  flow"  jets  easily  while  the  fuel  from 
the  four  off-axis  jets  seemed  to  accumulate  in  the  dome  head  region.  The  flow  field  did  not  changff 
significantly  over  the  entire  momentum  ratio  range  (Figure  13).  The  transport  “tube”  mechanisms 
existing  in  between  the  inlet  arms  for  the  60  and  90  degree  combustors  could  not  be  located  for  the 
45  degree  combustor. 

Alternating-Station  Combustors 

The  stagger-station  combustors  generated  a  very  different  flow  field  compared  to  the  same- 
station  models.  The  investigation  began  with  the  90  degree  combustor.  The  “inlet  flow”  and 
“fuel”  flow  rates  were  varied  as  in  the  same-station  combustor  runs.  The  dome  heights  were 
fixed  at  0.2D  so  that  two  of  the  fuel  jets  would  impinge  on  two  of  the  “inlet  flow”  columns 
effectively.  The  result  was  that  the  other  two  fuel  jets  would  impinge  directly  onto  the  wall  of 
the  combustor.  Since  this  behavior  would  increase  wall  erosion  greatly,  it  was  seen  that  an 
alternate  injector  would  have  to  be  fabricated  for  further  investigation  of  the  alternating-station 


combustors.  The  impingement  point  ot  the  two  upstream  inlet  tlow  columns  resulted  in  a 
■fanning''  effect.  .A.  portion  of  this  fanning  of  the  upstream  jets  would  then  divert  the 
downstream  jets  further  downstream.  The  fanning  effect  witnessed  on  all  alternating-station 
combustors.  Since  two  of  the  "fuel"  jets  impinged  on  the  model  wall  for  these  configurations 
and  the  level  of  mixing  appeared  low.  further  testing  was  not  performed. 

The  low  level  of  mixing  and  “fuel"  transport  with  these  configurations  was  believed  to  be  due 
to  the  absence  of  the  transport  tube  mechanisms  seen  in  the  same-station  combustors. 

CONCLUSIONS 

The  tlow  fields  observed  in  the  four-inlet  geometries  have  the  intrinsic  characteristic  of  diverting  a 
substantial  amount  of  the  inlet  flow  into  the  dome  region  of  the  combustors  resulting  in  two 
recirculation  regions.  Characteristics  of  the  recirculation  patterns  and  geometrical  limitations 
resulted  in  preferred  dome  heights  of  0.25D,  0.083D,  and  0  for  the  90,  60  ,  and  45  degree  same- 
station  combustors  respectively. 

The  90°  same-station  combustor  contained  strong  and  stable  recirculation  zones  both  upstream 
and  downstream  of  the  inlet  flow  inlet  locations.  The  upstream  recirculation  zone  interacted  with 
the  off-axis  fuel  jets  in  two  distinct  patterns.  It  was  seen  that  above  a  critical  momentum  ratio,  the 
fuel  jet  “bent  backwards”  and  followed  a  pattern  dictated  by  the  “inlet  flow”  flow  resulting  in  one 
primary  recirculation  zone  per  jet.  Once  the  momentum  ratio  was  lowered  below  a  critical 
momentum  ratio,  the  same  fuel  jet  impinged  and  penetrated  the  incoming  inlet  flow  column 
resulting  in  two  counter-rotating  recirculation  zones  located  on  either  side  of  the  fuel  jets.  In 


between  each  of  the  “inlet  flow”  inlet  columns  existed  two  counter-rotating  tubes  aligned  with  the 


combustor  axis.  These  tubes  were  believed  to  be  the  primary’  transpon  mechanisms  responsible  for 
transferring  the  dome  region  products  to  the  aft  combustor  chamber. 

The  60  degree  same-station  combustor  contained  similar  flow  structures  as  seen  in  the  90  degree 
same-station  combustor,  but  they  were  observed  to  be  less  organized  and  weaker  in  strength.  The 
dome  height  for  the  60  degree  combustor  was  reduced  to  0.083 D  in  order  to  maintain  strong 
recirculation  zones  in  the  dome  head  region.  The  axially  aligned  transfer  tubes  were  not  easily 
located,  but  effects  of  the  transfer  mechanisms  were  seen. 

None  of  the  45  degree  same-station  combustors  generated  a  favorable  flow  field.  The  shallow 
inlet  flow-inlet  angle  resulted  in  directing  most  of  the  "inlet  flow"  momentum  downstream, 
thereby  greatly  reducing  the  amount  of  mixing  which  occurred.  This  produced  a  fuel-rich  dome 
head  region.  After  viewing  the  three  same-station  combustors,  it  was  concluded  that  the  strength 
and  organization  of  the  recirculation  regions  and  tubes  were  highly  dependent  upon  the  amount  of 
momentum  directed  normal  to  the  combustor  axis.  Therefore  it  is  recommended  that  only  inlet 
flow  inlet  angles  of  60  degrees  or  greater  be  given  further  consideration. 

Although  actual  combustor  performance  cannot  be  determined  from  flow  visualization 
studies,  it  is  likely  that  the  conclusions  reached  from  the  observed  flow  fields  will  be  applicable 
in  designing  a  combustor  with  favorable  combustion  characteristics. 
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NOMENCLATURE 


D  ^  combustor  model  diameter 

U  =  mean  combustor  flow  velocity 

h  =  dome  height 

X  =  axial  location  along  combustor 

tha  =  mass  flow  of  air 

rhf  =  mass  flow  of  fuel 

rh^Va  =  momentum  of  air 

rhfVf  =  momentum  of  fuel 

Rep  =  Reynolds  number,  UD/u 

V  =  Kinematic  viscosity 
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Figure  i;  A  ducted  rocket  engine  schematic 


Figure  2:  Flow  Visualization  Facility 
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SIDE  VIEW 


Figure  3:  “Fuel”  injector  used  in  study 
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Figure  4:  Model  Combustors  Used 


Figure  5:  45  Degree  Same-Station  combustor  test 


Figure  6:  Recirculation  regimes  for  90  degree  Same-Station  combustor 


Figure  7:  Laser  Sheet  Orientation 


Figure  8:  View  B-B  of  90°  same-station  combustor 


Figure  10:  View  A-A  of  90°  same-station  combustor 


Figure  12:  View  A-A  of  60°  same-station  combustor 


re  13;  View  A-A  of  45°  same-station  combustor 
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Abstract 

Fuel  concentrations  existing  in  the  flameholding  region  of  a  four-inlet  ducted  rocket  engine 
with  a  five-orifice  showerhead  fuel  injector  are  presented  for  a  range  of  flow  rates  and  fiiel  injector 
orientations.  Methane  was  used  as  the  fuel  simulant  and  vibrational-Raman  spectroscopy  was  used 
to  image  the  methane  distribution  in  the  dome  head  region  of  an  optically  accessible  four-inlet 
combustor.  A  fuel  injector  orientation  which  produced  a  fuel  distribution  and  internal  flowfield 
desensitized  to  the  air/fiiel  momentum  ratio  was  found. 

Observations  made  during  the  gas-gas  mixing  study  agreed  well  with  conclusions  from 
previous  water  flow  visualization  studies  which  investigated  the  internal  flow  fields  of  similar  fuel 
injector  orientations  and  combustor  geometries. 


Nomenclature 


D 

Internal  Diameter  of  Combustor 

L 

Length  of  Main  Combustor 

/ 

Fuel  Mass  Fraction 

Wa 

Air  Mass  Flow  Rate 

OTf 

Air  Mass  Flow  Rate 

Va 

Air  Injection  Velocity 

Vf 

Fuel  Injection  Velocity 

0 

Air  Inlet  Arm  Angle 
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Introduction 


Ducted  rocket  engine  (DRE)  systems  have  recently  received  a  revived  interest  due  to  their 
high  specific  impulse,  throttlability,  and  suitability  for  long  range  applications.  They  obtain  high 
specific  impulse  by  inducting  the  air  through  which  they  fly  through  inlet  ducts 
and  using  it  as  the  oxidizer  which  is  mixed  with  liquid  fuel  or  fuel  rich  gas  from  a  gas  generator. 
The  mixture  combusts  and  is  then  accelerated  through  a  converging/diverging  nozzle  to  generate 
thrust.  A  diagram  showing  the  general  layout  of  a  DRE  is  shown  in  Figure  1 .  Much  of  the  work 
performed  in  the  past  has  been  concentrated  on  one-  and  two-inlet  geometries.  These  geometries 
offer  an  aerodynamic  advantage  since  the  inlet  ducts  may  be  tailored  to  generate  lift  for  the  vehicle 
but  often  require  some  type  of  inlet  flow  deflecting  device  to  be  installed  in  the  ducts  so  that 
sufficient  inlet  flow  enters  the  dome  region  of  the  combustor  where  the  flameholding  process 
typically  occurs.  The  complex  internal  flow  fields  of  these  combustors  have  often  been 
investigated  using  both  flow  visualization  with  water  for  qualitative  information  [Kennedy,  1974] 
and  direct-connect  hot-fire  testing  for  quantitative  results  [Schadow,  et  al.  1981]. 

Four-inlet  geometries  offer  advantages  over  the  one  and  two-inlet  geometries  such  as  steer-to- 
tum  capability  where  the  vehicle  may  turn  in  any  direction  while  keeping  at  least  two  air  inlet  arms 
exposed  to  the  windward  side  in  order  to  maintain  combustor  pressure.  One-  and  two-inlet 
geometries  must  bank-to-tum  in  order  to  accomplish  the  same  task.  The  four-inlet  geometry  also 
presents  the  potential  to  passively  divert  sufficient  inlet  flow  into  the  dome  region  of  the 
combustors  to  satisfy  flameholding  requirements  without  the  use  of  active  diverter  devices 
commonly  needed  on  other  combustors,  thereby  potentially  increasing  the  total  pressure  recovery 
and  improving  performance. 
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Figure  1:  A  Solid-Fuel  Ducted  Rocket  Engine  Schematic 


Figure  2:  Ducted  Rocket  Engine  Geometry  Definitions 


Ducted  rocket  geometries  have  been  successfully  developed  and  flown  demonstrating  their 
excellent  long-range  intercept  performance  [Meinkohn  and  Bagmen,  1982].  Although  some  data 
exist  for  a  selected  four-inlet  geometry  [Zetterstrom  and  Sjoblom,  1985],  additional  information  is 
needed  for  systems  with  conventional,  showerhead-type,  fuel  injectors.  Information  such  as 
flameholder  fuel  distribution,  fuel  residence  time,  and  flow  pattern  characteristics  for  this  fuel 
injection  scenario  are  needed  to  facilitate  further  development  of  these  systems.  The  work 
performed  in  this  study  increases  the  level  of  understanding  of  the  dome  head  flowfields  and  fuel 
distribution  existing  in  four-inlet  DRE’s  with  showherhead  fuel  injectors.The  geometrical 
definitions  used  in  this  paper  for  the  four-inlet  DRE  system  are  defined  in  Figure  2. 

Experimental  Setup 

The  cold  gas  and  reacting  flow  experiments  were  performed  on  a  direct-connect  blow  down 
facility  located  at  the  Propulsion  Directorate  of  the  U.S.  Army  Missile  Command.  The 
experimental  layout  is  shown  in  Figure  3.  The  facility  utilizes  a  500  cubic  feet  air  storage  system 
at  a  maximum  pressure  of  2500  psi  for  supply  air  which  can  be  delivered  at  flow  rates  up  to  10 
Ibm/sec.  The  large  air  capacity  affords  long  run  times  exceeding  6  minutes  in  duration,  even  at  the 
largest  flow  rate. 

The  vitiating  heater  uses  gaseous  hydrogen  as  a  fuel  and  can  provide  inlet  temperatures  up 
to  1500  °R  at  the  maximum  air  flow  rate.  The  replenishment  oxygen  is  added  upstream  of  the 
vitiating  heater  in  order  to  promote  more  complete  combustion  in  the  vitiator. 

A  high  pressure  air  ejector  is  used  to  provide  subatmospheric  exhaust  pressure  conditions. 
The  high  pressure  air  ejector  allows  nozzle  base  pressures  as  low  as  4  psia.  The  resulting  flight 
simulation  envelope  includes  Mach  numbers  from  1.9  to  3.1  at  sea  level  increasing  to  a  range  of 
2.55  to  3.7  at  a  maximum  simulated  altitude  of  65,000  feet  [Blevins,  et  al.  1995]. 
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Figure  4:  Optical  Combustor  Schematic 
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The  facility  was  supported  and  monitored  by  a  Keithley  12-bit  A/D  data  acquisition  system  which 
monitored  and  array  of  Validyne  and  Omega  pressure  transducers  and  6  oven-referenced  Omega 
thermocouples. 

The  methane  delivery  rates  were  obtained  by  monitoring  the  static  pressure  differential 
existing  across  mass  metering  venturies  fabricated  and  calibrated  by  Fox  Metering  Systems. 
Methane  static  pressure  and  temperature  measurements  were  made  across  the  venturi  bank  and 
immediately  prior  to  methane  injection  to  allow  precise  mass  flow  rate  measurements. 

Optically  Accessible  Combustor 

The  heart  of  the  optical  diagnostics  hardware  was  the  optically  accessible  combustor.  A 
schematic  of  the  combustor  is  shown  in  Figure  4.  The  combustor  was  composed  of  three  segments; 
the  dome  head,  the  main  body,  and  the  aft  section.  All  the  combustor  segments  were  made  of  317L 
stainless  steel  except  the  inlet  arms,  which  were  made  of  304L  stainless  steel.  The  dome  head  had 
a  center  hole  cut  in  it  with  an  accompanying  #26  0-ring  groove  for  the  fuel  injector.  A  MS33649 
standard  hydrogen  torch  port  was  cut  adjacent  to  the  hole  and  was  used  for  the  ignition  torch.  The 
dome  head  also  contained  a  1.5mm  by  30mm  slit  through  which  the  laser  sheet  was  directed.  A 
fused  silica  window  was  placed  into  the  cutout  recession  which  had  grooves  cut  into  the  base 
through  which  nitrogen  was  purged.  The  nitrogen  purge  was  used  to  keep  the  window  clean  and 
was  fed  through  a  pressure  sealing  face  which  bolted  to  the  dome  head  and  holds  the  fuel  injector 
and  window  in  place. 

The  main  body  had  an  internal  diameter  of  101.6  mm  and  a  wall  thickness  of  25.4  mm.  It 
contained  four  air-inlet  arms  with  an  internal  diameter  of  25.4  mm  and  a  wall  thickness  of  3.175 
mm  that  had  been  welded  in  place.  It  bolted  to  the  dome  head  by  eight  Allen  bolts  and  formed  a 
pressure  seal  through  a  combination  o-ring/ceramic  gasket  seal.  The  ceramic  gasket  had  an  I.D.  of 
101 .6  mm  and  an  O.D.  of  121.3  mm  and  protected  an  O-ring  seal  with  a  3.1  mm  thickness  and  an 
inner  diameter  of  121.6  mm.  The  main  body  contained  two  pinhole  camera  cavities,  of  which  one 
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housed  the  primary  collection  optics.  The  primary  collection  lens  was  a  biconvex  fused  silica  lens 
with  a  focal  length  of  17.0  mm  at  1^=400  nm  and  had  a  anti-reflective  coating  optimized  for  a 
wavelength  of  400  nm. 

The  cavity  was  capped  with  a  25.4  cm  viewport  window  which  had  a  thickness  of  10  mm. 
|N|i^i*ogen  was  purged  directly  below  the  window  and  was  allowed  to  flow  around  the  biconvex  lens 
in  order  to  keep  it  cool  and  clean.  Care  was  taken  to  generate  a  uniform  Nitrogen  flowfield  around 
the  lens  in  order  to  inhibit  any  form  of  circulation  to  occur. 

The  aft  section  of  the  combustor  had  a  wall  thickness  of  6.3  mm  and  was  bolted  to  the  main 
body  in  a  similar  fashion  as  the  dome  head.  It  contained  two  static  pressure  transducers  located 
within  50  mm  of  either  end  of  the  section.  Thermocouples  were  also  placed  on  this  section  to 
monitor  the  combustor  surface  temperature.  Flow  exited  the  aft  section  through  a  graphite  sonic 
nozzle. 

The  fuel  injectors  used  with  the  combustor  were  geometrically  similar  to  the  one  used  for 
the  water  flow  visualization  study  [Brophy  and  Hawk,  1996].  They  were  both  five-orifice  injectors 
made  of  304  stainless  steel  with  one  axial  orifice  and  four  off-axis  at  45  degrees  to  the  combustor 
axis  and  equally  spaced  at  90  degrees.  The  orifices  of  the  injectors  were  scaled  so  that  the 
momentum  ratio  range  investigated  in  the  water  flow  visualization  study  could  be  properly 
matched  in  the  cold  gas  and  reacting  flow  studies. 

Optical  Illumination  and  Collection  Optics 

The  laser  source  used  in  the  tests  was  a  Quantel  Nd;YAG  system.  Using  a  pair  of  INRAD 
mixing  ciystals,  230  mJ  at  the  third  harmonic  (355nm),  was  deliverable  to  the  combustor  at  a 
repitition  rate  of  5  Hz  and  a  pulse  width  of  approximately  10  ns.  All  of  the  mixing  crystals, 
harmonic  separators,  and  sheet  lighting  lenses  were  coated  with  anti-reflective  coatings  in  order  to 
provide  greater  than  99%  transmission  at  the  optimized  wavelength. 


S 


A  number  of  check-out  tests  were  performed  with  the  combustor  located  in  a  laboratory  on  the 
UAH  campus.  The  tests  were  performed  to  investigate  the  Raman  imaging  quality.  The  laser 
output  was  formed  into  a  sheet  with  a  height  of  about  25  mm  and  an  approximate  thickness  of  0.5 
mm  at  the  imaged  location  through  the  use  of  one  plano-convex  spherical  lens  and  two  cylindrical 
lenses.  The  sheet  was  directed  into  the  combustor  through  the  dome  head  and  illuminated  the 
center  plane  of  one  of  the  off-axis  fuel  jets. 

A  Princeton  Instruments  ICCD  576s  Intensified  CCD  (ICCD)  camera  was  used  to  collect  the 
CH)  Vibrational  Raman  image  signal  at  a  wavelength  of  396  nm.  The  camera  was  controlled  by  a 
Princeton  Instruments  ST  130  controller  and  a  Princeton  Instruments  PG  200  programmable  pulser. 
The  image  from  the  primary  collection  lens  was  turned  90°  after  passing  through  the  window  by  a 
fused  silica  prism.  The  image  was  then  filtered  so  only  the  selected  wavelengths  were  transmitted 
to  the  camera.  An  ORIEL  filter  centered  at  396nm  was  used  for  the  study  and  was  50.8  mm  in 
diameter.  The  filter  was  a  3  cavity  design  and  had  a  FWHM  of  10  nm  which  was  required  to 
minimize  the  attenuation  of  off-axis  rays  collected  in  the  image.  The  filter  was  attached  to  the  end 
of  a  105mm  f  2.8  NIKON  focusing  lens  attached  to  the  ICCD.  The  imaged  area  was  approximately 
30  mm  in  diameter. 


Results 

The  knowledge  gained  from  the  water  flow  visualization  work  was  used  to  direct  the  cold  flow 
tests  performed.  The  reader  is  directed  to  Brophy  and  Hawk  [1995,1996]  for  additional 
information  on  the  water-flow  visualization  results. 

The  images  presented  in  this  work  are  for  a  60  degree  four-inlet  same-station  combustor  with 
a  dome  height  of  0.22D.  The  “same-station”  description  implies  that  all  of  the  air  inlet  arms  are 
located  at  the  same  axial  location.  Figure  5  shows  a  90  degree  same-station  combustor  and  depicts 
the  location  of  the  imaged  area.  The  imaged  location  is  the  same  for  the  60  degree  combustor. 
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Figure  5:  Imaged  Area  Location  and  Sheet  Lighting  Orientation 


Cold  Flow  Experiments 

The  cold  flow  experiments  investigated  the  mixing  of  vitiated  air  and  gaseous  methane  (CH4) 
over  similar  momentum  ratios  and  fuel  injector  orientations  as  in  the  water  visualization  study.  A 
total  of  fourteen  tests  were  run  with  the  optically  accessible  combustor  which  had  60°  air  inlet 
arms,  a  dome  height  of  0.22D,  and  a  five  orifice  fuel  injector  modeled  after  the  water  flow 
visualization  injector.  Two  injectors  with  different  orifice  sizes  were  fabricated  for  the  combustor 
in  order  to  maintain  subsonic  injection  and  still  cover  the  momentum  range  investigated  in  the 
water  flow  visualization  study.  The  conditions  simulated  and  flow  rate  settings  for  each  run  are 
listed  in  Table  1 .  The  term  <j)  describes  the  overall  equivalence  ratio  of  the  combustor  and  is 
defined  as: 


Ft  A 

(r'A^ 


(1) 


All  combustor  configurations  had  sonic  graphite  nozzles  installed  with  various  area  ratios  to 
obtain  the  targeted  chamber  pressures.  Thrust  and  chamber  pressure  measurements  were  made  for 
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each  configuration  as  well  as  inlet  total  pressure,  inlet  static  pressure,  and  inlet  static  temperature 
for  the  vitiated  air.  Fuel  flow  was  controlled  by  a  dome  loaded  regulator  which  supplied  methane 
service  pressures  from  350  to  1000  psi  to  a  bank  of  three  mass  metering  venturies.  Any 
combination  of  venturies  could  be  opened  or  closed  which  provided  a  wide  range  of  fuel  delivery 
rates.  Two  fuel  injectors  were  fabricated  in  order  to  properly  cover  the  desired  mass  flow  and 
momentum  ratio  ranges.  Both  injectors  were  modeled  after  the  water  flow  visualization  fuel 
injector  and  had  one  center  fuel  jet  and  four  equally  spaced  off-axis  fuel  jets  angled  at  45°.  The 
“high  velocity”  fuel  injector  had  a  center  fuel  orifice  diameter  of  5.32  mm  and  four  off-axis  holes 
with  diameters  of  3.1  mm.  The  “low  velocity”  fuel  injector  had  a  center  fuel  orifice  diameter  of 
10.65  mm  and  four  off-axis  holes  with  diameters  of  6.53  mm. 

The  calibration  of  the  Raman  signal  was  performed  by  first  sealing  and  evacuating  the 
combustor  to  IkPa  and  then  pressurizing  it  with  methane  to  a  pressure  of  200.0kPa.  Laser  shots 
were  directed  through  the  quartz  window  in  the  dome  head  and  into  the  combustor  as  thin  laser 
sheets.  Image  intensities  were  obtained  as  a  function  of  pressure  thereby  obtaining  a  molar 
concentration  versus  image  intensity  calibration.  Since  the  methane  distribution  could  be  assumed 
to  be  uniform,  the  images  were  spatially  corrected  for  the  sheet  light  intensity  distribution, 
pincushioning  effects,  and  attenuations  associated  with  filtering  off-axis  rays. 

Most  of  the  images  obtained  from  the  cold  gas  study  are  an  accumulation  of  20  to  30  sequential 
laser  shots  which  were  acquired  over  a  5  to  6  second  time  period.  This  was  needed  in  order  to  obtain  signal 
to  noise  ratios  of  acceptable  levels  (>10).  The  ICCD  camera  was  synchronized  to  the  laser  Q-switch  output 
and  was  gated  for  a  shutter  speed  of  40ns.  The  camera  was  also  set  to  collect  a  specified  pixel  definition 
which  only  included  the  imaged  area.  This  reduced  the  scanning  time  and  allowed  the  camera  to  collect  an 
image  for  every  laser  pulse  when  the  laser  was  set  to  a  repetition  rate  of  5  Hz.  The  pixel  resolution  at  the 
center  of  the  images  was  0.17mm  x  0.17mm  x  0.5mm  and  decreased  at  the  perimeter  to  0.47mm  x  0.47mm  x 
0.5mm  due  to  the  pincushioning  effect.  The  presence  of  concentric  circles  near  the  perimeter  of  the  cold 
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Table  1:  Cold  Gas/Reacting  Flow  Test  Conditions  and  Settings 


*With  combustion 

flow  images  is  a  result  of  the  post-processing  of  the  images  to  correct  for  the  pincushioning  effect. 

Of  the  12  successful  tests  listed  in  Table  1,  four  produced  images  which  clearly  depict 
flowfield  structure  and  fuel  concentration  levels  for  the  cold  flow  tests.  A  “ghost”  image  of  the 
combustor  geometry  is  shown  with  the  images  presented  to  provide  a  point  of  reference.  The  fuel 
concentration  levels  are  represented  by  the  grey-scale  legend  to  the  right  of  the  images  which 
asscociates  an  f-parameter  value  to  each  portion  of  the  imaged  area.  The  f-parameter  is  also  known 
as  the  stoichiometric  mixture  fraction  parameter  and  is  defined  as 

A  (2) 

F  F 
- 1 - 

A  A  ST 
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t?  T7 

where  —  is  the  local  fuel  to  air  ratio  and  —  is  the  stoichiometric  fuel  to  air  ratio.  The  term  f  is 
A  A  ST 

convenient  since  it  can  only  vary  from  0  to  1.  When  f  is  between  0  and  0.5,  fuel  lean  conditions 
exist  and  a  value  between  0.5  and  1  therefore  corresponds  to  fuel  rich  conditions.  A  value  of  0.5 
corresponds  to  a  local  equivalence  ratio  of  1 .  The  flammability  limits  for  methane  and  air  at  STP 
conditions  are  between  equivalence  ratios  of  0.53  and  1.6.  These  values  correspond  to  f  parameter 
values  of  0.36  to  0.615  respectively  and  are  labeled  on  the  figures  for  convenience. 

The  following  sections  will  refer  to  two  different  recirculation  regimes  observed  for  the  in-line 
fuel  injector  orientation.  Both  of  the  recirculation  regimes  were  observed  and  labeled  during  the 
previous  water  flow  visualization  portion  of  this  work.  They  were  termed  passive  and  impinging. 
The  passive  regime  was  found  to  exist  at  high  air/fuel  momentum  ratios  when  the  off-axis  fuel  jets 
were  observed  to  “bend  back”  and  passively  follow  the  incoming  air  stream.  The  impinging 
regime  was  seen  to  occur  at  low  air/fuel  momentum  ratios  where  the  off-axis  fuel  jets  dominated 
the  flowfield  and  impinged  on  the  incoming  air  stream.  The  passive  and  impinging  regimes  are 
depicted  in  Figures  6  and  7  respectively. 


High  Air/Fuel  Momentum  Runs 

The  best  image  of  a  high  air/fiiel  momentum  ratio  condition  is  from  Run  3  and  is  shown  in  Figure 
8.  This  run  had  an  air/fuel  momentum  ratio  of  29.98  and  Figure  8  reveals  evidence  of  the  passive 


recirculation  condition.  Evidence  of  the  incoming  air  jet  can  be  seen  in  the  lower  right  comer  of 


Figure  8.  The  strong  fuel-deflecting  power  of  the  high  momentum  air  is  seen  to  force  the  fuel  to 
follow  the  large  recirculation  zone  in  the  dome  head.  The  disappearance  of  fuel  from  the  right  side 
of  the  imaged  area  is  believed  to  be  due  to  deflection  of  the  fuel  jet  out  of  the  illumination  plane  by 


the  incoming  air  column. 
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Figure  7:  Impinging  Recirculation  Regime 


14 


The  higher  air  flow  rate  runs,  2  and  3,  resulted  in  higher  overall  combustor  pressures  and 
therefore  yielded  higher  S/N  ratios.  The  lower  combustor  pressure  runs  (*174  kPa)  appeared  to 
reveal  similar  flow  field  features  and  methane  concentration  levels,  but  at  lower  S/N  ratios.  Signal 
to  Noise  levels  for  the  low  equivalence  ratio  runs  ((|)<0.4)  bordered  on  acceptable  at  times, 
obviously  due  to  the  overall  reduced  fuel  levels  present. 

Low  Air/Fuel  Momentum  Runs 

Run  2  provided  a  good  image  of  the  impinging  fuel  jet  regime  shown  in  Figure  9.  The  image 
clearly  shows  the  fuel  jet  reaching  the  air  inlet  jet.  Evidence  of  the  center  fuel  jet  can  be  seen  at  the 
top  of  the  image.  This  run  had  an  air/fuel  momentum  ratio  of  1 1 .3 1  and  generated  very  fuel-rich 
conditions  in  the  dome  head  region. 

Similar  pressure  and  overall  equivalence  ratio  effects  on  signal  quality  were  seen  for  the  low 
air/fiiel  momentum  ratio  runs  as  in  the  high  momentum  runs. 

Fuel  Injector  Orientation  Effects 

Figures  10  and  11  display  images  of  the  same  plane  as  the  previous  images,  but  for  these 
images  the  fuel  injector  was  rotated  45°  and  placed  in  the  interdigitated  orientation.  It  was  seen 
that  much  less  fuel  accumulated  on  plane  A-A  compared  to  the  in-line  fuel  injector  orientation. 
Figure  10  shows  the  fuel  concentration  levels  for  a  low  air/fuel  momentum  ratio  condition  in  Run 
8.  Figure  1 1  shows  a  similar  flowfield  for  a  higher  momentum  ratio  condition  in  Run  7.  Only  a 
slight  increase  in  fuel  concentration  for  the  fuel  rich  runs  was  seen  over  the  fuel  lean  condtions. 
This  reveals  that  the  interdigitated  fuel  injector  orientation  desensitizes  the  dome  head  flow  field  to 
the  overall  fuel  flow  rate  and  momentum. 
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Figure  8:  High  Air/Fuel  Momentum  Ratio  Image  of  Run  3.  In-Line  Fuel  Injector. 


Figure  9:  Low  Air/Fuel  Momentum  Ratio  Image  of  Run  2.  In-Line  Fuel  Injector. 
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Figure  10:  Low  Air/Fuel  Momentum  Ratio  of  Run  8.  Interdigitated  Fuel  Injector. 


Figure  11:  High  Air/Fuel  Momentum  Ratio  of  Run7.  Interdigitated  Fuel  Injector. 
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Discussion 


The  cold  flow  studies  involved  the  imaging  of  CEi  using  Raman  spectroscopy.  Since  the 
imaged  signal  strength  was  a  function  of  molar  concentration  (pressure),  the  higher  pressure 
conditions  yielded  higher  signal  to  noise  ratios  than  the  lower  pressure  runs.  This  was  seen  to  be 
very  important  for  the  lower  equivalence  ratios  investigated  since  lower  overall  methane 
concentrations  existed.  Figure  8  shows  results  from  Run  3  with  a  high  air/fuel  momentum  ratio  of 
29.98.  For  this  condition,  the  air  inlet  jet  dominantly  controlled  the  off-axis  methane  jet.  Notice 
that  the  fuel  jet  is  forced  to  follow  a  pattern  suggested  by  Figure  6  of  a  similar  momentum  ratio 
condition  simulated  by  the  water  flow  visualization  study.  The  large  well-defined  fuel 
concentration  gradient  seen  in  Figure  8  implies  that  the  flow  field  was  reasonably  stable  and  steady 
over  the  20  image  collection  period.  Evidence  of  the  air-inlet  jet  can  be  seen  on  the  right-hand 
side  of  the  Figure  8.  The  disappearance  of  CH4  from  the  right  side  of  the  image  was  believed  to  be 
due  to  some  degree  of  flow  directed  out  of  the  imaged  plane.  This  test  condition  generated  a  fuel 
rich  yet  combustible  mixture  in  the  dome  head  recirculation  region. 

Figure  9  shows  an  image  from  Run  2  which  had  a  much  lower  air/fuel  momentum  ratio  of 
11.13.  The  region  had  a  very  large  amount  of  methane  present  and  was  predominantly  fuel  rich, 
and  revealed  f-parameter  values  well  above  the  flammability  limits.  The  image  of  the  fuel  jet 
appears  to  indicate  an  impinging  type  of  recirculating  regime,  depicted  in  Figure  7.  Conclusions 
from  previous  water  flow  visualization  experiments[Brophy,Hawk  1997]  indicated  that  a 
momentum  ratio  of  11.13  should  result  in  the  recirculating  regime  fluctuating  between  impinging 
and  passive,  but  the  image  seems  to  depict  a  purely  impinging  regime.  This  is  believed  to  be  a 
result  of  the  compressible  gas  flow  versus  incompressible  water  flow  visualization. 

The  effects  of  compressibility  in  the  cold  flow  tests  were  believed  to  be  the  reason  for  the 
definitive  impinging  flow  regime  occurring  at  such  a  high  momentum  ratio  of  11.13.  Further 
testing  needs  to  by  performed  in  order  to  isolate  this  behavior  and  draw  firmer  conclusions. 
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Interdigitated  Fuel  Injector 


Cold  flow  images  for  the  inter-digitated  fuel  injector  are  shown  in  Figures  10  and  11.  Figure 
1 1  shows  the  results  from  run  7  which  was  for  an  air/fuel  momentum  ratio  of  28.0  and  Figure  10 
shows  the  results  from  run  8  which  had  a  momentum  ratio  of  12.64.  Although  Figure  10  represents 
a  much  higher  fuel  flow  rate,  there  is  not  much  difference  in  fuel  concentration  levels  on  plane  A-A 
between  the  two  images.  This  observation  implies  that  the  recirculating  regime  on  plane  A-A  is 
apparently  insensitive  to  the  overall  air/fuel  momentum  ratio  of  the  combustor.  This  is  believed  to 
be  due  to  the  majority  of  the  fuel  contained  in  the  off-axis  jets  being  injected  around  or  “in- 
between”  the  four  air-inlet  jets  and  into  the  aft  end  of  the  combustor.  The  slight  increase  in  fuel 
concentration  levels  was  expected  for  the  lower  momentum  ratio  conditions  and  was  believed  to  be 
due  to  the  portion  of  the  center  fuel  jet  which  was  redirected  back  into  the  dome  head  region. 

The  fuel  concentration  levels  seen  for  the  interdigitated  configuration  fell  within  the 
flammability  limits  and  were  seen  to  be  higher  at  the  extreme  front  end  of  the  combustor.  A 
similar  out-of-plane  motion,  as  seen  in  earlier  figures,  was  believe  to  be  responsible  for  the 
diminishing  levels  of  methane  on  the  right  side  of  the  images.  This  fuel  injection  scheme  appeared 
to  be  very  favorable  from  a  flameholding  standpoint  due  to  the  insensitivity  of  the  dome  head  fuel 
concentration  levels  to  the  momentum  ratios  evaluated  and  overall  combustor  equivalence  ratios. 

Figure  12  summarizes  observations  of  the  previous  water  flow  visualization  tests  for  three 
combustor  geometries  with  the  current  cold  flow  study  of  the  60-degree  same-station  combustor. 
All  information  presented  in  Figure  12  is  for  the  in-line  fuel  injector  orientation,  depicting  the  two 
recirculation  regimes  present. .  Select  cold  gas  test  conditions  are  shown  as  “hollow”  circles  on  the 
frgure.  Since  the  interdigitated  configuration  did  not  appear  to  have  any  distinct  recirculation 
regimes,  a  figure  for  that  configuration  is  not  appropriate. 
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Figure  12:  Recirculation  Regimes  for  Combustors  Evaluated  with  In-Line  Fuel  Injector 


Summary  and  Conclusions 

The  work  presented  in  this  paper  further  characterizes  the  internal  flow  fields  of  four-inlet 
ducted  rocket  engine  geometries.  Vibrational  Raman  spectroscopy  of  gas-gas  mixing  was  used  to 
provide  quantitative  information  on  the  fuel  distribution  for  various  configurations.  This  work 
compliments  previous  water  flow  visualization  experiments  to  provide  additional  internal  flow 
field  characteristics  of  four-inlet  geometries  investigated.  While  the  previous  water  flow 
visualization  tests  provided  valuable  qualitative  characteristics  on  such  flowfields,  the  cold  gas 
flow  tests  provide  quantitative  information  which  supports  prior  qualitative  conclusions  and  reveals 
local  fuel  concentrations  to  be  expected  on  such  systems  . 
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The  important  conclusions  reached  from  this  work  are: 


1)  The  use  of  Raman  spectroscopy  to  image  the  CH4  distribution  was  successful  and  yielded  the 
spatial  distribution  of  fuel  and  information  on  the  dome  head  flowfield  structure.  The 
technique  was  successfully  applied  to  cold  gas  mixing  flow  fields  of  the  optical  combustor. 
The  successful  application  of  the  technique  to  a  reacting  flowfield  is  possible  and  should  be 
considered. 

2)  Observations  of  flowfield  structures  for  the  cold  flow  gas-gas  mixing  study  agreed  well  with 
the  observations  made  from  a  previous  water  flow  visualization  study  and  indicated  that  a 
substantial  amount  of  inlet  air  is  passively  diverted  into  the  dome  head  for  four-inlet  DRE 
combustors  when  compared  to  one  and  two  inlet  geometries. 

3)  The  60-degree  same-station  combustor  was  observed  to  have  two  distinct  dome  head 
recirculation  regimes  with  the  in-line  fuel  injector  configuration  as  seen  in  previous  water  flow 
visualization  tests.  These  regions  were  termed  passive  and  impinging  and  were  seen  to  be  a 
strong  function  of  the  air/fiiel  momentum  ratio.  Fuel  rich  conditions  typically  existed  on  the 
imaged  plane  for  this  injector  orientation.  The  combination  of  flow  rate  dependent 
recirculation  zones  and  generally  fuel  rich  conditions  existing  over  a  majority  of  the 
investigated  range,  resulted  in  the  in-inline  fuel  injector  orientation  to  be  categorized  as 
undesirable. 

4)  The  interdigitated  fuel  injector  produced  combustible  fuel-air  mixture  concentrations  on  the 
imaged  plane.  The  fuel  concentrations  observed  in  this  configuration,  appeared  to  be 
desensitized  to  the  air/fuel  momentum  ratios  for  fuel  lean  and  fuel  rich  conditions.  The 
interdigitated  fuel  injector  was  determined  to  be  the  most  favorable  injector  orientation  when 
operating  the  combustor  over  a  wide  range  of  air/fuel  momentum  ratios. 
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Abstract 

Water  flow  visualization  and  preliminary 
combustion  studies  have  been  performed  in  order  to 
investigate  the  mixing  and  combustion  processes  inside 
of  four-inlet  side  dump  combustors  also  known  as 
ducted  rocket  engines  (DRE).  The  flow  fields  of  four- 
inlet  side  dump  combustor  geometries  were  visualized 
using  water  flow  visualization  and  documented  with 
video  and  conventional  pbotogrt^hy  so  that 
relationships  between  various  geomecric/flow 
parameters  and  the  flow  field  behavior  could  be 
obtained.  It  was  shown  that  for  inlet  angles  greater 
than  60  degrees,  the  geometries  intrinsically  direct  a 
significant  portion  of  the  inlet  flow  towards  the  dome 
head  region  without  the  need  for  diverter  devices.  The 
resulting  dome-region  flow  field  contains  two  distinct 
recirculation  patterns  shown  to  be  dependent  on  the 
momentum  ratio  and  inlet  flow  inlet  angle.  These 
parameters  are  believed  to  strongly  affect  the  existence 
and  strength  of  transport  mechanisms  between  the  fore 
and  aft  combustor  segments.  The  flow  visualization 
woik  concluded  that  the  60  degree  same-station 
combustor  produced  a  favorable  flow  field  over  the 
range  investigated  and  became  the  chosen  geomeuy  for 
the  reacting-flow  study.  The  geometry  was  used  to 
construct  an  optically  accessable  stainless  steel 
combustor  for  a  reacting  flow  investigation.  The 
reacting  flow  woik  utilizes  laser  induced  flutvescence 
to  view  the  flame  front  (CH)  distribution  and  will 
potentially  utilize  stimulated  Raman  spectroscopy  to 
view  the  cold/hot  flow  CH4  behavior.  Initial 
calibration  work  has  been  performed  on  campus. 
Direct-connect  testing  will  take  place  at  the  U.S.  Anny 
Missile  Command  this  summer. 
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Nomenclature 


D 

combustor  model  diameter 

u  = 

mean  combustor  flow  velocity 

h 

dome  height 

X  = 

axial  location  along  combustor 

liia  = 

mass  flow  of  air 

ihf  = 

mass  flow  of  fuel 

maVa  = 

momentum  of  air 

riifVf  = 

momentum  of  fuel 

Rbd  = 

Reynolds  number,  UD/\) 

V  = 

Kinematic  viscosity 

Introduction 

Ducted  rocket  engine  systems  have  recently 
received  a  revived  interest  due  to  their  high  specific 
impulse,  throttlability,  and  suitability  for  long  range 
applications.  They  obtain  the  high  specific  impulse  by 
ducting  the  inlet  flow  in  which  they  fly  through  inlet 
ducts  and  using  it  as  the  oxidizer  which  is  mixed  with 
a  fuel  rich  gas  from  a  gas  generator.  The  mixture 
combusts  and  is  then  accelerated  through  a 

converging/diverging  nozzle  to  generate  thrust.  A 
diagram  showing  the  general  layout  of  a  DRE  is  shown 
in  Bgure  1.  Much  of  the  work  performed  in  the  past 
has  been  concentrated  on  one-  and  two-inlet 
geometries.  These  geometries  offer  an  aerodynamic 
advantage  since  the  inlet  ducts  may  be  tailored  to 
generate  lift  for  the  vehicle  but  often  require  some  type 
of  inlet  flow  deflecting  device  to  be  instaUed  in  the 
ducts  so  that  sufficient  inlet  flow  enters  the  dome 
region  of  the  combustor  where  the  flameholding 

process  typically  occurs.  These  combustors  have  often 
been  investigated  using  flow  visualization  with  water 
for  qualitative  information’*^  and  direct-connect  hot- 
fire  testing  for  quantitative  results^*’. 

Four-inlet  geometries  offer  advantages  over  the  one 
and  two-inlet  geometries  such  as  steer-to-turn 

capability  and  the  tendency  to  passively  divert 
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Figure  1:  A  ducted  rocket  engine  schematic 


necessary  inlet  flow  into  the  dome  region  of  the 
combustors  for  flameholding  requirements  without  the 
use  of  diverter  devices  commonly  needed  on  other 
combustors.  The  effects  of  inlet  flow  angle,  dump 
station,  air/fuel  mass  ratios,  air/fuel  momentum  ratios, 
dome  height,  and  injector  performance  were 
investigated  over  a  wide  range  with  a  water  flow 
visualization  facility  and  resulted  in  a  60  degree  same- 
station  combustor  geometry  to  be  chosen  for  the 
combustion  diagnostics  tests. 

The  reacting  flow  experiments  will  duplicate  the 
momentum  and  mass  flow  range  investigated  in  the 
flow  visualization  study  in  order  to  compare  the  flow 
fields  observed  for  each  study.  Planar  Laser  Induced 
Fluorescence  and  Raman  spectroscopy  will  be  used  to 
view  the  flame  front  (CH)  and  fuel  (CR,)  distributions 
respectively.  Additional  information  will  be  obtained 
pertaining  to  the  flow  Structure  existing  near  the  flame- 
out  condition  for  four-inlet  same  station  side  dump 
combustors. 


Experimental  Setup 

Flow  Visualization  Facility 

The  flow  visualization  facility  (Figure  2)  has 
overall  dimensions  of  6.1  meters  by  4.6  meters.  Not 
shown  in  the  figure  are  the  38(X)  liter  reservoir  tank 
and  the  centrifugal  pump  located  below  the  deck.  The 
power  plant  fw  the  facility  is  a  22.4  kW  variable  RPM 
centrifugal  pump  capable  of  dehvering  3(X)0  liters  per 
minute  at  207  kPa  to  the  combustor  model.  This  pump 
was  dedicated  to  delivering  high  volume  at  low 
pressure  and  was  used  to  supply  water  as  the  "air" 
source  for  the  combustor.  A  supplemental  517  kPa 
water  supply  arm  was  used  to  deliver  the  water  for  the 
"fuel"  source  at  up  to  132  liters  per  minute.  The 
amniint  of  water  delivered  to  the  model  was  controlled 
both  by  the  pump  RPM  and  by  throttling  the  supply 
arms  with  globe  valves.  A  PC  monitored  Blancett  401 
flowmeters  used  on  each  arm,  two  pressure 
transducers,  and  provided  remote  pump  control  if 
desired.  The  combustor  models  were  placed  inside  a 
0.92m  X  l.CXhn  x  0.60m  viewing  tank  which  contained 
water  to  minimize  refraction  effects. 


The  fuel  injection  geometry  (Figure  3)  is  a  concave 
done  head  which  contains  the  fuel  injector  and  allows 
positioning  at  any  axial  location  in  the  front  end  of  the 
combustors.  The  dome  height  is  defined  as  the  axial 
distance  between  the  beginning  of  the  inlet  flow  dump 
station  and  the  perimeter  of  the  dome  head.  A  five 
orifice  injector  was  used,  one  axial  orifice  and  four  off- 
axis  orifices  angle  at  45  degrees  to  the  combustor  axis. 

Visualization  of  the  flow  field  was  achieved  by 
seeding  the  reservoir  tank  with  0.0254  cm  (0.010") 
polystyrene  spheres.  The  fuel  jet  was  left  unseeded. 
“Outgassing”  of  air  in  the  “fuel”  supply  after 
experiencing  the  pressure  drop  across  the  injector 
resulted  in  a  mist  of  tiny  inlet  flow  bubbles  which 
provided  adequate  visualization  of  the  fuel  flow.  A  5 
watt  Argon  Ion  laser  was  used  to  illuminate  the  flow 
field.  The  laser  output  was  first  focused  and  then 
expanded  through  a  cylindrical  lens  to  produce  a  light 
sheet  ^jproximately  1.0  mm  thick.  The  sheet  was 
capable  of  being  oriented  either  vertically  or 
horizontally.  The  images  observed  were  captured  on 
Ifi8  video  and  by  conventional  photogr^hy. 

Preliminary  Hot-Rre  Setup 

A  101.4  mm  internal  diameter  optically-accessable 
stainless  steel  combustor  was  constructed  based  on 
conclusions  from  the  flow  visualization  woiic.  The 
combustor  has  four  same-station  60  degree  inlet  arms 
and  a  dome  height  of  1.27  cm.  A  30.0  mm  x  0.75  mm 
slot  exists  in  the  dome  head  of  the  combustor  to  allow  a 
laser  sheet  to  enter  the  combustor.  Viewing  of  the 
illuminated  area  was  performed  by  employing  a 
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pinbole-iike  camera  design  (Figure  4).  This  design 
allows  for  a  constant  nitrogen  flow  over  the  primary 
collection  optic  in  order  to  keep  it  clean  during  future 
direct-connect  tests.  The  amount  of  nitrogen  purged 
through  the  aperture  will  be  typically  less  than  0.30% 
of  the  total  flow  rate  of  the  combustor. 

The  preliminary  combustor  tests  have  been 
performed  in  a  lab  on  campus  in  order  to  calibrate  the 
collection  optics  and  also  to  become  familiar  with  the 
combustor  and  the  methods  used.  The  combustor  was 
located  on  an  optical  table  and  burned  methane 
through  the  fuel  injector  to  generate  both  laminar  and 
turbulent  diffusion  flames.  The  simple  flame 
geometries  allow  calibration  and  evaluation  of  the 
optical  techniques  to  be  used.  The  physical  layout  is 
shown  in  figure  5. 

The  third  harmonic  of  a  Continuum  SL-10 
Nd:YAG  was  used  to  optically  pump  a  Continuum 
ND60  dye  laser  containing  a  Coumarin  120  dye. 
Typical  output  energies  were  6  mJ  at  431.5  nm  over  a 
5  ns  pulse  length.  The  output  from  the  dye  laser  was 
vertical  polarized  and  possesed  a  bandwidth  of 
approximately  .05  cm*‘.  The  beam  first  passed  through 
a  1.0  m  focal  length  plano-convex  lens  and  then  two 
cylindrical  lenses  to  create  a  laser  sheet  2.0  cm  wide 
and  0.3  mm  thick  which  was  then  directed  into  the 


105mm  focal  length  lens  at  a  f-stop  of  2.8.  The 
imaged  area  was  approxiamateiy  3.0  cm  by  1.4  cm. 

Results  and  Discussion 

The  flow  visualization  studies  investigated  two 
types  of  combustors,  same-station  and  alternating 
station  combustors.  All  were  six-inch  I.D.  models  and 
consisted  of  45, 60,  and  90  degree  inlet  flow  inlet  arm 
dump  angles  each  with  an  internal  diameter  of  3.81 
cm.  All  were  investigated  up  to  a  Reynolds  number  of 

5.5  X  10 ,  but  no  Reynolds  number  effects  were 
observed  above  Reynolds  numbers  of  1.5  x  10'^.  The 
mass  flow  ratio  (ihg/lilf)  and  momentum  ratio 

(ihjVg  /rtlfVf )  were  varied  between  10  to  60  and 

0.43  to  2.58  respectively. 

Dome  Height 

Dome  height  was  the  furst  parameter  investigated. 
The  first  combustor  tested  was  the  90  degree  same- 
station  combustor.  A  dome  height  of  0.25D  (3.81  cm) 
resulted  in  symmettic,  stable  recirculation  regions  for 
the  90°  same-station  combustw  over  a  wide  range  of 
momentum  ratios  and  appeared  to  be  the  preferred 
height  for  that  combustor  and  the  fuel  injector  used 
based  on  the  size  and  strength  of  the  recirculation 
regimes.  The  preferred  height  was  determined  after 


combustor  through  the  dome  head. 

The  collected  image  was  collimated,  filtered  by  a  9 
nm  bandpass  filter  centered  at  490  nm,  and  then 
refocused.  A  Princeton  Instruments  intensified  CQ> 
camera  collected  the  image  through  a 


testing  each  dome  height  setting  over  the  range  of 
momenuim  ratios  Specified  earlier.  Figure  6  shows  the 
operating  zones  for  momentum  ratio  versus  dome 


Figure  3:  “Fuel”  injeaor  geometry 


Figure  5:  Preliminary  Hot-Fire  Setup 
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height  for  the  90  degree  combustor.  The  preferred 
dome  height  was  then  chosen  as  the  height  which 
minimized  the  fluctuating  behavior  between  the 
impinging  and  passive  recirculation  regimes  while  also 
avoiding  wall  impingment  of  the  fuel  jet.  Dome 
heights  greater  than  0.27D  resulted  in  the  fuel  jets 
impinged  on  the  combustor  wall  and  not  the  incoming 
inlet  flow  columns.  The  preferred  dome  height  for  the 
60®  same-station  combustor  was  determined  to  be 
=0.083D  (1.27  cm)  by  the  same  manner.  An 
acceptable  dome  height  for  the  45®  combustor  was  not 
found,  so  a  dome  height  of  0  cm  was  used  in  order  to 
provide  the  strongest  recirculation  possible.  The 
alternating-station  combustors  had  problems  with  two 
of  the  fuel  jets  impinging  on  the  wall  because  the  same 
fuel  injector  was  used  on  all  models. 


for  Reynolds  numbers  above  5.5  x  10\  but  fluctuations 
in  size  and  shape  were  observed. 

The  low  “fuel”  flow  rate  resulted  in  a  momentum 
ratio  of  2.39.  For  the  high  momentum  ratios,  the  fuel 
jet  was  seen  to  “bend  back”  and  follow  the  governing 
“inlet  flow”  flow.  This  was  termed  a  “passive” 
flameholding  scenario.  There  was  no  significant 
change  in  this  flow  field  until  a  momentum  ratio 
between  1.90  and  1.61  was  reached.  In  this  range,  the 
fuel  jet(s)  fluctuated  between  the  “passive”  scaiario 
and  a  straight  fuel  jet,  termed  “impinging”,  which 
impinged  on  the  incoming  “inlet  flow”  column.  Below 
a  momentum  ratio  of  1.61,  the  fuel  jet  would  partially 
peneuate  the  inlet  flow  column  and  then  bifurcate  into 
two  branches.  One  would  be  deflected  towards  the 
wall  and  the  other  towards  the  stagnation  point  of  the 
four  jets.  The  result  was  two  recirculation  zones  split 


■Samp.-.Station  Combustors 

The  90  degree  same-station  combustor  was 
investigated  first  with  a  low  “fuel”  flow  rate.  The  flow 
field  contained  symmetric  counter-rotating  zones  on 
both  sides  of  the  inlet  flow  inlet  station  for  plane  A-A. 
Figure  7  shows  the  orientation  of  the  sheet  lighting 
relative  to  the  combustor.  The  recirculation  zones 
increased  in  strength  as  the  “inlet  flow”  flow  rate 
increase,  but  the  downstream  zone  appeared  to 
fluctuate  in  size  and  shape.  Figure  8  shows  the  pairing 
of  two  counter-rotating  regions  found  between  each  of 
the  inlet  arms  for  plane  B-B.  These  regions  appeared 
to  be  tube-like  and  were  seen  to  contain  substantial 
amounts  of  “fuel”.  The  endpoints  of  the  tube  could  not 
be  isolated,  but  appeared  to  begin  in  the  dome  head 
region  and  terminate  approximately  1.0  inlet  diameter 
downstream  from  the  inlet  flow  station.  The  structures 
seen  in  the  flow  field  did  not  appear  to  develop  further 


Figure  6:  Recirculation  regimes  for  90  degree 


and  driven  by  the  fuel  jet. 

The  60  degree  combustor  was  investigated  at  a 
dcrnie  height  setting  of  0.083D  and  over  the  same 
momentum  ratio  range.  The  recirculation  regimes 
observed  near  the  dome  head  on  the  90  degree 
combustor  also  occurred  in  the  60  degree  combustor 
but  were  seen  to  tiansifion  at  a  higher  momentum 
ratio.  The  “passive”  regime  (Figure  9)  existed  down  to 
a  momentum  ratio  of  2.05,  fluctuated  in  the  range 
between  2.05  and  1.72,  and  became  “impinging”  below 
1.72  (Figure  10).  The  critical  momentum  ratio  values 
are  defined  as  the  value  at  which  the  fuel  jet  transitions 
from  passive  to  impinging.  This  was  taken  as  the 
center  value  in  the  fluctuating  region.  The  critical 
momentum  ratios  for  the  60  and  90  degree  combustors 
are  shown  to  be  similar  when  the  only  lateral 
momentum  component  of  the  “inlet  flows”  were 
compared. 

The  downstream  recirculation  regions  and 


reattachment  lengths  remained  relatively  constant  over 
the  momentum  ratio  range  investigated,  but  did  reduce 
in  size  once  the  center  fuel  jet  penetrated  the 
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Figure  8:  View  B-B  of  90®  same-station  combustor 


inlet  arms  t^tpeared  to  be  weaker  and  shorter  in  length 
than  observed  in  the  90  degree  combustor.  No  fuel 
accumulation  in  the  dmne  head  region  was  observed 
for  this  combustor,  indicating  that  the  transport 
mechanisms  were  indeed  functioning. 

Tbe  45  d^iee  same-station  combustor  exhibited 
the  weakest  recirculation  zones.  This  was  likely  due 
to  the  reduced  amount  of  inlet  flow  momentum  normal 
to  flie  ounbustor  axis  which  was  seen  to  produce  much 
of  the  mixing  and  transport  that  occurs  between  the 
"fuel"  and  "air  inlet  flow”  jets.  Ibe  center  fuel  jet  was 
sear  to  penetrate  the  "inlet  flow"  jets  easily  while  the 
fuel  from  the  four  off-axis  jets  seemed  to  accumulate  in 
the  dmne  head  region.  The  flow  field  did  not  change 
significantly  over  the  entire  momentum  ratio  range 
The  transport  mechanisms  existing  in  between 
the  inlet  arms  fnr  the  60  and  90  d^ree  combustors 
could  not  be  located  for  the  45  degree  combustor. 

Alternating-Station  Combustors 

The  alternating-station  combustors  did  not  produce 
favorable  flowfields.  For  this  reason,  they  will  not  be 
further  addressed  in  this  paper.  For  a  mote  complete 
review  of  the  flow  visualization  work  and  the 
alternating  station  combustor  bdiavior,  the  reader  is 
directed  to  reference  10  and  11. 

Preliminary  Reacting  Flow  Results 

Initial  combustion  diagnostics  were  performed  on 
an  optics  table  at  campus.  By  viewing  a  simple 
methane  diffusion  flame  through  the  collection  optics, 
optical  access  was  verified.  The  low  f-number  primary 
collection  optic  resulted  in  the  imaged  being  warped  in 
a  convex  manner.  A  conformal  mapping  code  is  being 


Figure  10:  View  A-A  of  60°  same-station  combustor 


written  to  correct  for  the  deformation  by  calibrating  off 
of  a  target  slide.  Figure  11  shows  the  geometry  to  the 
calibration  study.  A  45  degree  methane  diffusion 
flame  was  used,  and  a  laser  sheet  penetrated  through 
the  center  of  the  jet  firom  below. 

The  goals  of  the  combustion  diagnostics  tests  are  to 
be  able  to  view  the  fuel  (CH4)  and  flame  front  (CH) 
bdiavior  simultaneously  in  the  combustor.  To  do  this, 
a  minimum  of  two  cameras  will  eventually  be  needed. 
At  this  point  in  time,  only  one  camera  is  available  and 
being  used  to  evaluate  various  techniques  that  could  be 
used  for  these  goals. 

The  first  method  investigated  in  this  study  in  one  in 
which  the  (0,0)  A^A-X^fl  band  of  (TH  is  excited  at 
431.5nm  and  the  fluorescence  signal  is  observed 
through  the  (0,1)  band  at  489nm.  The  problem  with 
this  method  is  that  the  Av^  are  only  qiproximately 
2%  of  the  diagonal  transitions  and  therefore  are 
typically  weak‘^.  Another  potential  problem  with  this 
method  is  that  there  is  a  Stokes  Raman  signal  existing 
for  methane  at  493.6  nm  for  the  same  excitation. 
Although  the  Raman  signal  was  seen  to  be  extranely 
weak  at  atmospheric  pressures  and  the  energies  used  in 
this  study,  the  signal  will  become  significant  at  higher 
combustor  pressures.  Reference  10  presents  methods 
to  separate  these  signals. 
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Figure  11:  Calibration  Flame  Geometry  combustor 


The  dye  laser  was  tuned  to  43 1 .5  nm  and  the  output 
energy  was  typically  around  6  mJ  over  5  ns.  The  sheet 
was  directed  into  the  combustor  and  the  signal  imaged 
by  the  collection  optics.  A  9  nm  bandwidth  filter 
centered  at  490  nm  was  used  to  filter  the  unwanted 
light,  and  the  Princeton  Instruments  intensified  camera 
was  gated  for  a  50  ns  shutter  speed. 

To  estimate  the  amount  of  background  signal 
collected,  the  camera  was  exposed  to  the  imaged  area 
while  the  laser  pulse  was  blocked.  No  image  was  seen. 
This  implied  that  the  images  collected  were  laser 
induced.  The  first  image  taken  was  that  of  a  laminar 
methane  flame  shown  in  figure  12.  The  Iowct  portion 
of  the  flame  is  believed  to  be  CH  fluorescence,  while 
the  upper  portion  of  the  flame  appears  to  be 
contaminated  with  what  is  believed  to  be  laser  induced 
incandescence.  This  was  determined  by  detuning  the 
laser  off  of  431.5  nm  and  viewing  the  altered  image. 
Further  studies  will  be  run  to  investigate  the  possibility 
of  separating  the  two  signals,  and  to  determine  exactly 
what  they  are. 

Figure  13  shows  a  turbulent  diffusion  flame 
(attached)  which  clearly  shows  the  fuel  core  and  the 
neighboring  reaction  zone.  Figure  14,  shows  a  lifted 
turbulent  diffusion  flame.  Definite  flow  structure  can 
be  seen  for  this  flame. 


Figure  12:Laminar  diffusion  flame 


Figure  11:  Turbulent  diffusion  flame 
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Figure  14:Turbuient  lifted  diffusion  flame 


Summary 

The  flow  fields  observed  in  the  four-inlet  geometries 
have  the  intrinsic  characteristic  of  diverting  a 
substantial  amount  of  the  inlet  flow  into  the  dome 
region  of  the  combustors  resulting  in  two  recirculation 
regions.  Giaracteristics  of  the  recirculation  patterns 
and  geometrical  limitations  resulted  in  preferred  dome 
heights  of  0.25D  and  0.083D  fm  the  90  and  60  degree 
combustors  respectively.  A  preferred  dtxne  height  was 
not  found  for  the  45  degree  combustor. 

The  90”  same-station  combustor  contained  strong 
and  stable  recirculation  zones  both  upstream  and 
downstream  of  the  inlet  flow  inlet  locations.  The 
upstream  recirculation  zone  interacted  with  the  off-axis 
fuel  jets  in  two  distinct  patterns.  It  was  seen  that  above 
a  critical  momentum  ratio,  the  fuel  jet  “bent 
backwards”  and  foOowed  a  pattern  dictated  by  the 
“inlet  flow”  flow  resulting  in  one  primary  recirculation 
zone  per  jet  Once  the  momentum  ratio  was  lowered 
below  a  critical  momentum  ratio,  the  same  fuel  jet 
impinged  and  penetrated  the  incoming  inlet  flow 
column  resulting  in  two  counter-rotating  recirculation 
zones  located  on  either  side  of  the  fuel  jets.  In  between 
each  of  the  “inlet  flow”  inlet  columns  existed  two 
counter-rotating  tubes  aligned  with  the  combustor  axis. 
These  tubes  were  believed  to  be  the  primary  transport 
mechanisms  responsible  for  transferring  the  dome 
region  products  to  the  aft  combustor  chamber. 

The  60  degree  same-station  combustm  contained 
similar  flow  structures  as  seen  in  the  90  degree  same- 


station  combustor,  but  they  were  observed  to  be  less 
organized  and  weaker  in  strength. 

None  of  the  45  degree  same-station  combustors 
generated  a  favorable  flow  field.  The  shallow  inlet 
flow-inlet  angle  resulted  in  directing  most  of  the  “inlet 
flow”  momentum  downstream,  thereby  greatly 
reducing  the  amount  of  mixing  which  occurred.  This 
produced  a  fuel-rich  dmne  head  region.  It  is 
recommended  that  only  inlet  flow  inlet  angles  of  60 
degrees  or  greater  be  considered  for  application. 

The  preliminary  combustion  tests  have  begun  to 
explore  various  diagnostic  techniques  fcM*  view  the  fuel 
and  flame  front  distribution.  Early  results  are 
promising,  but  higher  signal  levels  will  be  needed  for 
actual  direct-connect  testing.  This  requires  that 
additional  transition  bands  be  explored  in  order  that 
the  best  method  may  be  employed.  The  possible 
detection  of  laser  induced  incandescence  introduces  a 
problem  that  must  be  addressed  before  further  tests  are 
made.  Ideas  for  solving  some  of  these  problems  are 
presented  in  the  next  section. 

Future  Work 

Future  work  will  include  CH  imaging  by  using  a 
stronger  transition  band.  Allen,  M.  et  al.^^  concluded 
fhat  the  maximiim  fluoiescence  signal  could  be 
obtained  by  exciting  at  428  nm  and  observing  the 
fluorescence  at  431  nm.  The  only  problem  with  this 
mmhod  is  that  a  very  narrow  bandpass  filter  would  be 
needed,  and  care  must  given  to  off-axis  rays. 

The  camera  gating  sequence  will  also  be 
reevaluated  in  order  to  guarantee  the  ieast  amount  of 
noise  is  introduced  into  the  signal. 

Finally,  once  a  satisfactory  experimental  method  is 
chosen,  the  optical  combustor  will  be  relocated  to  the 
Propulsion  Directorate  of  the  U.S.  Army  Missile 
Command  for  direct-connect  testing. 
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an^e,  air  inlet  location,  air/foel  n^s  t^os,  ^ 

5  Sg 

^Lrational  range.  Information  gamed  fr^  the 
operational  nmge.  visualizanon 

qualitative  woik  pafoimed  m  the  onticallv 

facility  will  be  used  in  the  design  of 
'^Me  sub-acal.  comhuam.  “f  f"' 

connect  testing  with  methane  as  a  fuel  smmlanL 

P„r.ility  Description 
The  water  flow  visualizafion 
dimensions  of  6  meters  by  4.6 
is  located  on  the  UAH  campus.  A  t^  view  ot  m 
facility  showing  the  overall  I^yo^ 

2  Not  shown  in  the  figure  are  the  3780  li^  0^ 
gallon)  reservoir  tank  and  the  cenu^Sf 
Scared  below  the  deck.  The  power  pl^t  to  Ae 
fadlity  is  the  30  hp  variable  RPM  centnfug  P  P 
^Te  of  delivering  3080  Ut«s  ^ 

5^)  a.  206  H-a  (30  psi)  . 

This  pump  is  designed  to  dehver  a  voi 
low  pressure  and  is  used  to  supply  w 
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Fuel  Rich  Propellant 

\ 

Air  Intake 
/ 

Nozzle 

- - - — — 

7^  ^ 

Gas  Generator  Variable 

Fuel  Valve 

Ram  Burner 

Figure  1 :  A  ducted  rocket  engine  scbemadc 


the  "air"  source  for  the  combustor.  A  supplranental 
517  kPa  (75  psi)  water  supply  arm  is  used  to  deliver 
the  water  for  the  "fuel''  source  at  iq)  to  132  litas  per 
miniitp.  (35  gpm).  The  amount  of  water  delivaed  to 
the  model  is  controlled  both  by  the  pump  RPM  and 
by  throttling  the  supply  arms  with  ^obe  valves.  A 
PC  is  used  to  monitor  the  Blancett  W1120  flow 
meters  used  on  each  inlet,  two  pressure  transducers, 
and  provide  extanal  pump  control  if  desired. 

The  combustor  models  are  placed  made  a  91cm 
X  X  61cm  viewing  tank  which  contains  water 

to  minimize  refraction  effects.  The  inlet  arms  of 
model  interface  with  the  siqqily  arms  by  PVC 
unions  which  allow  fa  quid:  connect/disconnecL 
The  fuel  atm  is  connected  with  the  rear  of  the  dome 
head  with  a  section  of  PVC  and  two  Fanco 
couplings.  The  dome  head,  shown  in  figure  3,  is  a 
hemiqjherical  piston  head  with  0-rings  that  contains 
the  fuel  injecta  and  is  permitted  at  any  axial 
location  in  the  front  end  of  each  combusta.  The 
dnnnft  height  is  defined  as  the  axial  distance  betweai 
the  air  inlet  dump  station  and  the  lip  of  the  dome 
hftad  It  may  be  adjusted  by  varying  the  PVC  length 
between  the  Femco  couplings.  IMfferent  fuel 
injecta  designs  may  be  installed  into  the  dome  head 
by  simply  loosening  two  alien  saews  and 
substituting  the  new  core  into  a  5  cm  diameter 
receptacle  and  retightening  the  saews. 

Visualization  of  the  flow  field  is  achieved  by 
eitha  seeding  the  fuel  jet  a  the  air  jets  with  0.0254 
cm  (0.010")  polystyroie  spheres.  Seeding  of  the  fuel 
jet  is  acconq)lished  by  pressuring  a  highly  seeded 
hffTh  of  the  sphaes  and  bleeding  them  into  the  fiiel 
arm  through  a  gate  valve.  This  method  produces  die 
best  tracking  of  the  fuel  distribudon.  The  air  jets  can 
be  seeded  by  simply  placing  the  beads  into  the  large 
reservoir  up  to  a  lO^pardcles/m^  density.  The  fuel 
jet  is  then  left  unseeded,  but  does  illuminate  well  due 
to  cavitadon  of  the  fuel  stream  after  experiencang  the 
pressae  drop  of  the  injector.  The  latter  methcxl 
appears  to  work  best  fa  long  daadon  testing.  A  5 
watt  Argon  Ion  laser  is  used  to  illuminate  the 
polystyrene  sphaes  and  air  bubbles.  The  lasa 


Figure  3:  The  dome  head  and  fud  injeda  installed 
in  the  90“  same-stadon  combustor,w/mjeaor  details 


output  is  first  focused  and  then  expanded  through  a 
cylindrical  lens  to  produce  a  light  sheet 
apprriTimatdy  1.5  mm  thick.  The  sheet  may  be 
oriented  eitha  vadcdly  a  horizontally  and  its 
locadon  controlled  by  moving  the  sheet  lighting 
assembly  with  a  stqppa  motor.  The  images  are  then 
recaded  on  both  Hi8  video  t^  and  by  convendottal 
photognqthy. 


Initial  Results 

The  fatality  has  evaluated  3  of  the  6  ckosen 
ducted  rocket  engine  ctanbustors,  each  having  an 
inner  diameta  of  15.24  cm  (6”).  The  combustors 
evaluated  to  date  are  the  45, 60,  and  90  degree  same- 
stadon  inlet  combustas  shown  below  in  figure  4. 
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Parameter 

Range  1 

Mass  flow  ratios 

s 

VI 

VI 

o 

Momentum  ratios 

0. 05  <20.0 

ihfVf 

Domehdght 

0<d^  <20.0cni 

Air  inlet  angles 

45®,  60®,  and  90* 

Fuel  iniectors 

Limited  to  a  5  cm  diameter 

All  have  been  investigated  up  to  a  Reynolds  number 
of  5.5  X  10^  and  over  the  range  of  mass  and 
momentum  ratios  described  in  the  facility 
description  section. 

The  flow  field  for  the  90  degree  same-station 
combustor  produced  an  apparoitly  fuel-rich  region 
between  the  dome  head  and  the  inlet  dump  station 
for  most  of  the  flow  rates.  Although  this  region  did 
have  recirculation,  the  apparent  lack  of  "air" 
indicates  that  this  inlet  dump  configuration  is  not 
ideal.  Figure  5  shows  the  flow  field  behavior  for  the 
90  degree  same-station  dump  combustor  for  three 
different  flow  rate  conditions.  The  photogr^hs  are 
of  sheet  lighting  along  the  combustor  axis  at  mid 
plane.  The  bottom  picture  in  figure  5  is  for  a 
Reynolds  number  of  2.2  x  10^,  a  mass  flow  ratio  of 
25.0,  and  a  dome  height  of  7.00  cm.  There  is 
evidence  that  the  fuel  jet  is  penetrating  the 
stagnation  point  of  the  four  "air"  jets  and  generating 
significant  mixing  after  the  dump  station.  The 
problem  still  remains  of  the  apparently  fuel  rich 
region  in  the  front  end  of  the  combustor  which 
would  likely  lead  to  flame  extinction  once  the  rich 
limit  is  reached. 

The  60  degree  combustor  generated  a  pair  of 
stable  counter-rotating  regions  on  each  side  of  the 
angled  fuel  jets.  The  recirculation  regions  existed 


than  this  value  produced  regions  which  became 
unstable  and  periodic.  The  evidence  of  increased 
"air"  in  the  dome  head  region  implies  a  more 
suitable  combustion  xnixture  thereby  improving 
flame  holding  and  auto-ignidon  characterisdcs. 
Figure  6  shows  the  60  degree  combustor  with  three 
mass  flow  ratios.  Nodce  that  the  center  fuel  jet 
appears  to  penetrate  the  stagnadon  point  of  the  “air” 
idlets  for  the  lower  mass  flow  rado. 

The  45  degree  same-stadon  combustor  speared 
to  have  similar  but  weaker  recirculadon  zones  than 
the  60  degree  same-stadon  combustor.  This  is  most 
likely  due  to  the  amount  of  momentum  introduced 
normal  to  the  combustor  axis  by  the  "air”  inlets.  The 
reduced  amount  of  momentum  normal  to  the 
combustor  axis  obviously  reduces  the  amount  of 
mixing  that  occurs  between  the  "fuel”  and  "air”  jets, 
therefore  increasing  the  necessary  combustor  length 
for  conq)lete  combustion.  The  center  fuel  jet  appears 
to  penetrate  the  "air”  jets  more  easily  than  the  other 
combustors,  which  is  expected.  There  also  appears 
to  be  a  larger  "slow"  zone  around  a  fast  moving  core 
downstream  of  the  air  inlet  station  for  this 
combustor,  emphasizing  the  reduced  amount  of 
mixing  for  this  configuration. 


^  “ 


36.0' 

i 


<b) 
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Figure  4:  Ducted  rocket  engine  model  combustors  with  various  air  injection  geometries 
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U=  1.139  m/s 
ReD=2.17xl05 
iha  /  ihf  =  37.5 
niaVa  /  ihfVf  =  16.32 
h=  7.00  cm 


U=  1.146  m/s 
Re£)=  2.18  X 10^ 
rila  /  lilf  =  30.0 
ihaVa  /  ihfVf  =  10.44 
h»  7.00  cm 


U=  1.154  m/s 
ReD=2.20xl05 
riia  /  ihf  =  25.0 

h=  7.00  can 


Figure  5: 90°  same-station  combnstor  flow  fidd  at  different  fuel  flow  rates 
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0.S  cm  upstream  of  air  inlets 


0.5  cm  off  injector  face 


Overhead  view 


Figure  8:  Other  combustor  viewpoints 
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The  pafonnance  and  initial  results  from  the 
ducted  rocket  engine  flow  visualization  facility 
indicate  that  it  will  be  extremely  valuable  in 
contributing  to  the  DRE  research  program  underway 
at  UAH.  The  information  provided  by  the  various 
seeding  and  laser  lighting  techniques  yidded  insight 
into  possible  combustor  characteristics  such  as  flame 
holding  and  auto  ignition  characteristics. 

The  seeding  of  the  "air"  inlets  with  polystyrene 
spheres  and  allowing  the  "fuel"  injector  to  cavitate 
produced  an  easily  visualized  flow  field.  The 
viewing  tank  constructed  for  the  models  worked 
excq)tionally  well  at  minimizing  refraction  effects, 
therefore  allowing  acceptable  viewing, 
photographing,  and  video  recording  of  the  combustor 
flow  Adds. 

The  conclusions  reached  from  work  performed 
with  the  facility  wiU  be  used  in  choosing  the 
geometry  of  the  optically  accessible  combustor  to  be 
constructed  for  direct-connect  testing.  Results  from 
the  direct-cormeci  test  will  provide  valuable 
information  about  the  mixing,  flame  holding,  and 
auto-ignition  properties  of  ducted  rocket  engines 
over  a  wide  operating  range  while  revealing  the 
validity  of  using  flow  visualization  to  gain  insight 
into  four-inlet  DRE  combustor  flow  flelds. 
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The  MICOM  Ducted  Rocket  Engine  Research  Program  is  divided  into 
the  following  task  with  the  assigned  personnel  as  shown: 

Task  1.0  Ballistic  Properties  of  Booster  Propellant 
Candidates  with  Dr.  Robert  A.  Frederick,  Jr.  as  P.I.  assisted  by  Mr. 
Brian  Greiner,  GRA. 

Task  2.0  Benchtop  Erosion  Simulator  has  been  delayed  due  to 
reduced  funding. 

Task  3.0  Solid  Fuel-Gas  Generator  Fuel  Candidates  Ignition 
and  Combustion  Characteristics  with  Dr.  Douglas  Feikema  as  P.I. 
assisted  by  Mr.  John  Dempsey,  GRA. 

Task  4.0  Ducted  Rocket  Engine  Critical  Test  Parameters 
Identification  with  Dr.  Hugh  Coleman  as  P.I.  assisted  by  Mr.  John 
Blevins,  GRA. 

Task  4.2  Verification  of  Critical  Parameters  Through  Flow 
Visualization  with  Dr.  Clark  Hawk  as  P.I.  assisted  by  Mr.  Chris 
Brophy,  GRA. 

Task  5.0  Program  Management  and  Support  with  Dr.  Clark 
Hawk  as  Program  Director  assisted  by  Ms.  Linda  Marion, 
Administrative  Assistant 

This  fifth  quarterly  report  contains  a  overall  program  management 
summary  followed  by  detailed  discussions  of  accomplishments  and 
plans  for  the  technical  tasks  in  numerical  order. 
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Labscale  Hybrid  Uncertainty  Analysis 

B.  Greiner*  and  R.  A.  Frederick,  Jr.'f 
Propulsion  Research  Center,  University  of  Alabama  in  Huntsville 
Huntsville,  Alabama  35899 


Abstract 

A  past  investigation  by  The  University  of  Alabama 
in  Huntsville  Propulsion  Research  Center  into  the 
mechanisms  of  hybrid  rocket  instability  led  to  a 
detailed  experimental  uncertainty  analysis  of  a 
labscale  hybrid  motor  used  in  the  study.  This 
uncertainty  analysis  was  carried  out  to  determine  the 
quality  of  the  experimental  data  obtained  from  the 
labscale  device  and  the  applicability  of  labscale  hybrid 
testing  in  general.  In  particular,  calculations  of  the 
uncertainties  in  the  measurements  of  fuel  regression 
rate,  oxidizer  flux,  motor  characteristic  velocity,  and 
the  oxidizer-to-fuel  mass  ratio  were  made.  These 
calculations  were  performed  using  the  accepted 
methodology  of  Coleman  and  Steele. 

The  results  showed  that  the  typical  uncertainties  in 
the  measurements  of  the  fuel  regression  rate  (±8.7%), 
oxidizer  flux  (±10.8%),  motor  characteristic  velocity 
(±10.5%),  and  the  oxidizer-to-fiiel  mass  ratio  (±9.2%) 
were  relatively  high,  thus  limiting  the  use  of  the 
labscale  hybrid  motor  for  quantitative  investigations. 
The  uncertainty  levels  were  dominated  by  1)  high 
measurement  uncertainties  in  the  diameters  of  the 
initial  fuel  grain  initial  port  and  the  nozzle  throat, 

2)  a  large  conceptual  bias  in  the  burn  duration, 

3)  a  conceptual  bias  in  the  average  chamber  pressure, 
and  4)  a  high  precision  uncertainty  showing  low 
repeatability  of  some  results.  These  high  levels  of 
uncertainty  lead  to  limited  utility  of  labscale  data  in 
the  generation  regression  rate/oxidizer  flux  curve  fits 
and  to  low  resolution  in  the  determination  of 
characteristic  velocity  efficiency. 

Methods  that  were  identified  to  reduce  the  level  of 
uncertainty  in  the  data  reduction  results  were  to: 
1)  increase  the  scale  of  the  test  motor,  reducing  the 
uncertainties  in  the  measurements  of  the  diameters  of 
the  fuel  grain  port  and  nozzle  throat;  2)  increase  the 
bum  duration,  which  would  reduce  the  domination  of 
the  bum  duration  conceptual  bias;  and  3)  reduce  the 
uncertainty  in  the  average  chamber  pressure 
measurement  through  either  an  increase  in  the  test 
pressure  or  in  the  accuracy  of  the  calibration 
methodology,  both  of  which  were  believed  to  have 
questionable  or  limited  effectiveness.  A  parametric 
study  of  these  changes  is  recommended  to  determine 
the  magnitude  of  expected  improvements  in  the 
accuracy. 


*  Graduate  Research  Assistant,  AIAA  Student  Member 
t  Assistant  Professor.  Mechanical  and  Aerospace  Engineering,  Senior 
Member  AIAA 


Nomenclature 

Ba  =  bias  (fixed)  uncertainty  in  parameter  a 
C*  =  characteristic  velocity 

Cd  =  metering  venturi  discharge  coefficient 

Dj  =  initial  fuel  grain  port  diameter 
Dth  =  throat  diameter  of  motor  nozzle 
Invent  =  throat  diameter  of  metering  venturi 
Gox  =  gaseous  oxygen  mass  flux 
i  =  Aft  mixing  section  length 

L  =  fuel  grain  total  length 

Am  =  change  in  fuel  grain  mass 

mj  -  initial  mass  of  fuel  grain 

n\f  -  final  mass  of  fuel  grain 

%  =  mass  ratio  of  oxidizer  to  fuel 

=  precision  (random)  uncertainty  in  a 
Pc  =  average  chamber  pressure 

Pup  =  average  pressure  upstream  metering  venturi 
At  =  change  in  fuel  grain  port  radius 

f  =  average  fuel  regression  rate 

R  =  gas  constant  for  oxygen 

At  =  change  in  time 

t,  =  ignition  time 

tf  =  shutdown  time 

Tup  =  average  temperature  upstream  of  venturi 

Ua  =  total  uncertainty  in  component  a 
Xi  =  test  input  parameter 

Greek 

y  =  specific  heat  ratio  for  oxygen 

\i  =  mean  value 

It  =  pi,  3.141592654... 

Pf  =  density  of  solid  fuel  component 

0  =  standard  deviation 


Introduction 

In  order  to  reduce  the  cost  and  complexity  of 
experimental  programs  quite  often  experimenters 
choose  to  use  subscale  or  labscale  test  rigs  in  place  of 
full-scale  testing.  From  these  labscale  tests,  data  are 
retrieved  in  order  to  gain  more  information  about  the 
basic  operational  phenomenea  and  physics.  This 
methodology  has  been  widely  applied  in  the  area  of 
propulsion  testing  for  the  past  fifty  years.  However, 
until  recently,  little  more  than  afterthought  has  been 
given  to  the  experimental  uncertainty  of  the  data  from 
which  any  empirical  correlations  are  made.  This  is 
particularly  the  case  with  labscale  testing  surrounding 
hybrid  rocket  motors.* 

A  test  program  conducted  by  The  University  of 
Alabama  in  Huntsville  Propulsion  Research  Center 
into  the  mechanisms  of  hybrid  rocket  instability  led 


Copyright  ©  1995  by  the  American  Institute  of  Aeronautics  and  Astronautics,  Inc.  All  rights  reserved. 


1 

American  Institute  of  Aeronautics  and  Astronautics 


to  a  detailed  experimental  uncertainty  analysis  of  a 
labscale  hybrid  motor  used  in  the  investigation The 
primary  purpose  of  the  uncertainty  analysis  was  to 
determine  the  quality  of  the  experimental  data 
obtained  from  the  labscale  device.  However,  the 
results  of  the  analysis  were  also  used  to  make  a 
determination  of  the  overall  applicability  of  labscale 
hybrid  testing.  In  the  process,  several  key  parameters 
were  identified  which  dominated  the  uncertainties  in 
the  data  reduction  equations.  This  paper  details  the 
experimental  uncertainty  analysis  that  was  performed 
on  the  Labscale  Hybrid  Motor  (LHM)  by  describing: 
1 )  the  experimental  motor  hardware  and  data 
acquisition  system,  2)  the  development  of  the  data 
reduction  and  uncertainty  equations,  and  3)  the 
analysis  of  the  results  from  these  equations.  This 
will  suggest  parametric  study  of  the  key  data 
reduction  parameters  to  identify  the  optimum  test 
conditions  for  the  highest  quality  hybrid  data. 

Experimental  Setup 

Motor  Hardware 

The  LHM  hardware  consisted  of  a  forward  closure,  a 
grain  holder,  and  an  aft  closure,  as  shown  in  Figure  1 . 
A  gaseous  oxygen  feed  system  was  connected  to  the 
motor  at  the  forward  closure.  The  grain  holder 
surrounds  cartridge  loaded  fuel  grains  which  are  held 
in  place  by  a  retaining  washer  at  each  end  of  the  grain 
holder.  The  aft  mixing  section  consisted  of  an  aft 
closure  and  a  carbon  insert  surrounding  a  tungsten 
nozzle.  The  LHM  offered  flexibility  in  selecting 
various  lengths  of  aft  mixing  section,  £  in  Figure  1 , 
different  nozzle  diameters,  and  different  grain  lengths. 

During  operation  of  the  LHM,  gaseous  oxygen  was 
injected  in  the  forward  end  of  the  motor  and  diffused 
by  the  forward  closure  diffuser  screen.  The  gaseous 
oxygen  flow  rate  was  controlled  by  a  dome-loaded 
pressure  regulator  upstream  of  a  metering  venturi. 
After  a  10-second  oxygen  pre-flow,  a  current  applied 
to  a  squib  in  the  forward  closure  ignited  the  fuel. 
After  a  bum  time  of  approximately  two  seconds,  the 
oxygen  flow  was  terminated  and  a  nitrogen  purge- 
flow  was  initiated  to  ensure  complete  extinguishment 
of  the  fuel  surface.  Short  bum  times  were  used  due  to 
the  thermal  limits  of  the  uncooled  tungsten  nozzles 
used  for  these  experiments. 


The  primary  type  of  fuel  were  used  in  the  test 
firings  was  hydroxyl-terminated  polybutadiene 
(HTPB).  Standard  operating  procedure  called  for  four 
grains  which  were  cartridge  loaded  in  the  fuel  grain 
holder.  Each  grain  was  nominally  2.5  inches  long, 
had  a  0.820-inch  diameter  port  and  a  1.375-inch  outer 
diameter.  The  HTPB (R-45)/N- 100  fuel  grains 
(C6.939H9.855N0.178O0.264)  were  mixed  and  cast  under 
a  vacuum  in  phenolic  sleeves.  The  castings  were 
then  cured  for  24  hours  at  140  °F. 

Data  Acouisition  System 

The  instrumentation  used  for  the  LHM  experiments 
supplied  data  on  the  motor  pressures,  temperatures, 
and  oxygen  flow  rate.  Forward  and  aft  pressure 
transducers  were  used  to  measure  the  pressures  inside 
the  motor.  The  pressure  and  temperature  of  the 
incoming  oxygen  were  measured  upstream  of  a 
choked  metering  venturi  to  determine  the  oxidizer 
mass  flow  rate.  The  various  measurements  were 
recorded  on  a  digital  system  (0-40  Hz)  and  on  an 
analog  FM  tape  system. 

The  motor  chamber  pressure  data  were  obtained 
from  the  test  analog  data  recordings.  The  analog  data 
was  passed  through  a  5  kHz  low-pass  anti-aliasing 
filter  and  then  digitized  at  10  kHz.  This  high- 
frequency  data  was  calibrated  using  the  low-frequency 
data  since  the  calibration  techniques  used  on  the 
NASA/MSFC  analog  system  were  found  to  be 
unacceptable.  Instead,  a  technique  of  calibration  was 
selected  which  consisted  of  selecting  two  pairs  of 
corresponding  points  from  the  high  frequency  data 
(voltage-time)  and  the  low  frequency  data  (pressure- 
time)  from  which  a  linear  calibration  curve  was 
constructed.  The  technique  of  calibrating  using  the 
low  frequency  data  introduced  a  ±10psia  bias 
uncertainty  into  the  high  frequency  data  due  to  the 
scatter  of  the  high  frequency  data  at  the  calibration 
points.  This  data  was  used  for  qualitative  visual 
analysis  as  well  as  for  digital  analysis  of  the 
frequency  content  through  the  use  of  Fast  Fourier 
Transform  (FFT)  techniques.  The  motor  operating 
times  were  also  determined  from  the  high-frequency 
digital  data.  A  precision  balance  was  used  to  weigh 
each  fuel  grain  pre-  and  post-test  in  order  to  determine 
the  average  regression  rate. 


Figure  1.  Sketch  of  Thiokol/UAH/MSFC  LHM  Layout. 
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Analysis  of  Labscale  Uncertainty 
Using  methods  established  by  Coleman  and 
Steele,^  an  experimental  uncertainty  analysis  was 
performed  to  determine  the  uncertainty  in  the 
experimental  results  and  to  determine  the  important 
data  reduction  parameters  in  LHM  testing.  This  also 
enabled  an  evaluation  of  the  overall  applicability  of 
LHM  testing.  The  basic  uncertainty  technique 
consists  of  determining  the  sensitivity  of  the  data 
reduction  equations  to  each  input  parameter, 
multiplying  by  the  parameter  bias  limit,  calculating 
the  precision  limit  from  the  standard  deviation  of  the 
test  data,  and  combining  these  values  in  a  root-sum- 
square  form  to  determine  overall  uncertainty.  In  order 
to  carry  out  this  analysis,  the  equation  for  the 
experimental  value  of  interest  must  be  written 
explicitly  in  terms  of  measurable  or  input  parameters. 
The  development  and  results  of  these  calculations  are 
shown  in  the  following  sections  for  all  the  data 
reduction  equations. 


Reduction  of  Data 

Test  data  were  reduced  to  determine  the  actual 
oxygen  mass  flux  (Gox)»  average  fuel  burning  rate 
(f ),  characteristic  velocity  (C*),  and  oxidizer-to-fiiel 
ratio  (%)  during  each  test.  The  following  equations 
express  each  of  these  results  in  terms  of  data  measured 
from  each  test.  The  full  development  of  the  data 
reduction  equations  used  in  this  paper  are  detailed  in 
Reference  4. 

The  oxidizer  mass  flux  Gox  was  determined  from 
measurements  of  the  average  temperature  (Typ)  and 
pressure  (Pup)  upstream  of  the  oxygen  venturi,  and 
the  initial  and  final  propellant  mass  (mj  and  mf): 


^CdpQpDvcntY*^ 

r  - \ 

+ 

1 _ 1 

1  JtpfL  J 

where  L  is  the  length  of  the  fuel  grain,  Dj  is  the 
initial  grain  port  diameter,  pf  is  the  solid  fuel 
density,  Q  is  the  venturi  discharge  coefficient,  and 
Dvent  is  the  venturi  throat  diameter.  The  ideal  gas 
constant  (R)  and  ratio  of  specific  heats  (y)  are  for 
oxygen. 

For  the  calculation  of  regression  rate,  the  mass 
consumption  was  assumed  to  occur  evenly  only  over 
the  fuel  grain  bore  surface  during  the  firing.  This 
assumption  is  reasonable  because  of  the  short  bum 
times  used  in  this  investigation  but  under  the 
additional  assumption  that  the  ignition  and  shutdown 
transients  were  small  compared  to  the  bum  times. 


This  leads  to  the  following  equation  for  the  average 
regression  rate: 


r  s= 


^  =  1 


4(mi  —  mf) 


nPfL 


+  Di^-D: 


At 


(2) 


where  t^  and  tf  are  the  ignition  and  shutdown  times. 

In  studies  of  instabilities  it  is  sometimes  useful  to 
examine  combustion  efficiency.  This  may  be 
accomplished  by  dividing  the  actual  motor 
characteristic  velocity  C*  by  the  theoretical  C*  at  the 
same  actual  %  ratio.  The  overall  % 
investigation  was  found  by  taking  the  ratio  of  the 
oxidizer  and  fuel  mass  flow  rates.  This  leads  to  the 
equation  for  the  %  ratio  of: 


oxidizer  mass  flow  rale 


The  expression  for  the  motor  characteristic  velocity 
was  developed  from  the  definition  of  C*  which  is  the 
product  of  the  average  chamber  pressure  (Pc)  and  the 
nozzle  throat  area,  divided  by  the  total  mass  flow  rate. 
This  leads  to  the  final  data  reduction  equation: 


C*  =  - 


throatarea 
c  4  th 


C^P„pJtDL 


(4) 


Y+l 


oxidizer  mass  flow  rate 


mj  —  mf 

fuel  mass  flow  rate 


Generalized  Uncertainty  Analysis 

The  first  phase  of  an  uncertainty  analysis  is  the 
development  of  expressions  for  the  values  of  interest 
from  the  data  reduction  equations.  For  brevity,  the 
following  expressions  for  the  experimental 
uncertainty  of  the  LHM  are  detailed  using  the  gaseous 
oxidizer  frux.  Equation  (1),  as  an  example.  These 
uncertainty  expressions  can  then  be  examined  to 
determine  the  critical  or  dominant  test  parameters. 

First,  the  basic  equation  for  the  uncertainty  must  be 
considered.  For  Gox»  the  generalized  uncertainty 
equation  is  written: 


Ug„  =  (5) 
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where  Bq  represents  the  total  bias  or  fixed 

uncertainty  limit  in  Gox  which  includes  the 
contributions  to  the  overall  uncertainty  that  are 
constant  from  test  to  lest.  Pg„  represents  the 

precision  or  random  uncertainty  limit  for  Gox  which 
is  the  contribution  due  to  the  statistically  random 
"scatter"  in  the  data.  The  concepts  of  and 
are  best  demonstrated  by  the  illustration  shown  in 
Figure  2. 


Figure  2.  Illustration  of  Bias  and  Precision 
Uncertainty  Limits. 

Figure  2  shows  a  Gaussian  distribution  of  flux 
readings  centered  about  a  point  offset  from  the  Uiie 
value  of  Gox  by  an  amount  equal  to  the  bias  limit 
Bg  .  The  methodology  used  to  determine  the  bias 

limit  is  discussed  in  the  following  paragraphs.  The 
precision  limit  Pq^^  is  found  by  determining  the 

standard  deviation,  or  o,  of  the  data  about  the  mean 
value  p  at  a  single  set  point  and  then  doubling  it  to 
obtain  the  95%  confidence  or  ”2c"  interval.  This  is 
the  limit  within  which  one  can  be  confident  that  95% 
of  the  readings  will  lie.  For  this  investigation,  the 
scatter  in  results  from  test  to  test  was  considered  to  be 
due  to  slight  uncontrollable  differences  in  test 
parameters  such  as  the  venturi  upstream  pressure  ^d 
temperature.  This  scatter  is  trapped  into  the  precision 
limit  estimate  by  calculating  the  value  of  o  for  the 
test  data. 

In  order  to  calculate  the  bias  limit  for  Gox»  one 
must  determine  how  the  uncertainty  of  each  input 
parameter  propagates  into  the  bias  limit.  As  stated 
earlier,  this  is  accomplished  by  first  writing  the  data 
reduction  equation  for  the  result  of  interest  explicitly 
in  terms  of  measurable  input  parameters.  The  data 
reduction  equation  for  Gox  is  repeated  here  for 
convenience: 


2  C  ^  P  upD  y ' 


G.  = 


r%> 


Ly+1 


«PfL  y 


(1) 


up 


In  examining  Equation  (1),  it  is  seen  that  the  input 
parameters  are:  C<i,  Pup*  Dvent*  ^f,  Pf’ 

and  T„p.  From  Reference  3,  the  expression  for  the 
propagation  of  the  bias  uncertainties  of  the  data 
reduction  equation  input  parameters  Xj  into  Gox  is 
written: 


where  the  partial  derivative  term  is  simply  the 
sensitivity  of  Gox  to  parameter  xi  and  is  determined 
for  each  parameter  from  the  basic  data  reduction 
equation  for  Gox*  Equation  (6)  has  the  additional 
assumption  that  correlated  biases  between  parameters, 
if  they  exist,  are  negligible.  Expanding  Equation  (6) 
using  Equation  (1): 


Note  that  for  the  fuel  grain  mass  terms  a  summation 
is  included  to  account  for  the  bias  limits  in  the 
individual  fuel  grains. 

An  indication  of  the  relative  importance  of  each 
input  parameter  can  be  determined  by  conducting  a 
"zero-order"  uncertainty  analysis.  In  this  analysis  ^e 
precision  uncertainties  are  neglected  and  the  bias 
uncertainties  in  the  individual  parameters  are  assumed 
to  be  1%  of  the  nominal  parameter  value.  The  partial 
derivatives  were  calculated  to  determine  the  parameter 
sensitivities.  Due  to  the  relative  complexity  of  the 
data  reduction  equations,  a  finite  differencing 
technique  was  used  to  approximate  these  derivatives. 
These  calculations  were  carried  out  for  each  test  and 
each  data  reduction  equation  using  a  spreadsheet 
utility.  Table  1  shows  the  results  of  these 
calculations  for  the  individual  Gox  parameters  on  a 
nominal  HTPB  test. 
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Table  1.  Zero*Order  Uncertainty  Analysis 
Results  for  Gaseous  Oxidizer  Flux  (Gpjj)  using 
1%  Nominal  Uncertainties. 


Xi 

dG 

o» 

8x, 

1%0f 

Nominal 

dG 

— s*-B 

3x, 

Cd 

0.327 

in^  sec 

0.02 

0.0031 

in^  sec 

^up 

fbty 

0.0002 

psi 

11  psi 

0.0024 

in^  sec 

'-£151' -3^^,  sec  j 

:W:;;  ^  '  I f  ^ 

'"S' 

/A. 

.a6o6i5'4^^ 

r 

0.0763 

.  2 

in  sec 

0.0139 

0.001 1 

.  2 

in  sec 

mj 

\bty 

-0.0016  /'”^sec 
gm 

0.50  gm 

_  ___ _  Ibm 

-0.0008  - 

in^  sec 

m, 

tony' 

0.0016  '' 

gm 

0.41  gm 

0.0006 

.  2 

in  sec 

Pt 

1.673  ' 

‘/o 

0.0003  — 
in® 

0.0006 

in^  sec 

L 

tony' 

0.0056 

in 

0.1  in 

AM--  ibm 
0.0006 - 

in*  sec 

tW 

'  V  '  "  ,  .iirt  '“w 

....  '  ' 

.0:0382  lr», 

S'*-- 

V 

-0.0003 

R 

5.25  R 

-  ^  Ibm 

-0.0015 - 

in^  sec 

By  combining  the  terms  in  the  last  column  of 
Table  1  according  to  Equation  (7),  the  overall  bias  for 
Gqx  using  the  zero-order  analysis  is  ±0.0093 
Ibm/in^'sec  or  ±3.0%  of  nominal.  Comparing  this 
to  the  results  in  the  last  column  of  Table  1 ,  one  can 
see  that  the  uncertainty  in  is  strongly  dominated 
by  the  propagation  of  uncertainties  in  Dvcm  and  Dj. 
Uncertainties  in  C^,  P„p,  and  y  are  minor 
contributors  while  the  remaining  input  parameters 


(mj,  mf,  pf,  L,  and  Typ)  have  only  a  weak  influence. 
This  analysis  indicates  that  by  minimizing 
uncertainties  in  Dv^nt  and  Dj,  and  to  some  extent  C^, 
Pup*  and  y,  the  uncertainty  in  Gox  is  minimized.  The 
results  also  show  that  the  other  parameters  can  be 
neglected  so  long  as  the  magnitudes  of  their 
uncertainties  are  equal  to  or  less  then  those  of  the 
dominant  parameters  on  a  percent  of  nominal  basis. 
Table  2  at  the  bottom  of  this  page  shows  the  results 
of  the  zero-order  analysis  on  all  of  the  data  reduction 
equations.  The  double  checks  represent  the  dominant 
input  parameters,  single  checks  are  the  minor 
contributors,  blanks  show  a  weak  influence,  and  the 
dash  indicates  that  the  particular  data  reduction 
equation  is  independent  of  that  input  parameter. 

Detailed  LHM  Uncertainty  Analysis 
The  next  step  in  the  examination  of  the  LHM 
uncei^nty  was  to  perform  the  detailed  level  analysis. 
This  is  a  refining  type  analysis  that  is  performed  once 
test  data  is  taken  which  permits  the  calculation  of  the 
precision  index,  or  2a,  and  a  better  estimate  of  the 
system  biases.  Table  3  below  shows  representative 
values  for  the  precision  limits  for  the  different  data 
reduction  results  as  well  as  the  precision  uncertainty 
as  a  percent  of  the  nominal  resultant  value. 


Table  3.  Representative  Precision  Limits  in 
_ the  Experimental  Results. _ _ 


Experimental 

Result 

Precision 

Limit,  Pa 

%of 

Nominal 

GoxO) 

±0.00489 

in*  sec 

±1 .6 

r{2) 

±0.001158  — 

sec 

±2.9 

%(3) 

±0.296 

±5.8 

C*(4) 

±317.2 

sec 

±6.8 

Table  4  shows  the  estimates  of  the  input  parameter 
bias  limits  in  absolute  magnitude  as  well  as  a 
percentage  of  actual  value  for  a  representative  HTPB 
test.  The  estimates  of  the  individual  biases  were 
made  at  a  95%  confidence  level.  The  biases  for  Dvcm 
and  Cd  were  based  on  the  venturi  manufacturer’s 
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Table  4.  Estimate  of  Input  Parameter  Bias 
Limits  for  LHM  Petailed  Analysis  of 


Input 

Parameter 

Bias 

Limit,  Ba 

%  of  Actual 

Cd 

±0.02 

±2.1 

Dyent 

±0.001  in 

±1.0 

P  UD 

±5  psi 

±0.5 

Tud 

±10  R 

±1.9 

y 

±0.05 

±3.6 

m 

±0.0005  qm 

±0.0001 

Pf 

±0.004  lbm/in3 

±11.9 

L 

±0.05  in 

±0.5 

Di 

±0.05  in  ^ 

±6.1 

D,h 

±0.005  in 

±1.6 

Pc 

±10  psi 

±2.1 

ti 

±0.035  sec 

±1.8 

tf 

±0.115  sec 

±5.8 

specifications.  For  Pup  and  Tup,  biases  were  b^ed  on 
the  calibration  specifications  and  on  fluctuations  of 
conditions  during  tests.  The  bias  in  y  was  based  on 
scatter  in  published  values  for  oxygen.  The  bias  for 
the  mass  measurements  was  based  on  the  accuracy  of 
the  digital  analytical  scale  used.  The  fuel  grain  length 
bias  was  taken  from  the  accuracy  of  the  calipers  used 
to  measure  the  grains.  For  the  fuel  density,  the  bias 
was  estimated  based  on  scatter  in  accepted  density 
values.  The  bias  in  the  initial  port  diameter  was 
estimated  based  on  the  casting  mandrel  tolerances. 
The  bias  in  Dth  was  based  on  drawing  tolerances. 

As  stated  earlier,  a  bias  of  ±10  psi  was  introduced 
in  the  measurement  of  chamber  pressure  due  to  the 
methodology  used  for  calibration  of  the  test  data. 
However,  an  additional,  more  dominant  bias  due  to  a 
conceptual  uncertainty  in  the  definition  of  chamber 
pressure  was  present  in  the  testing  as  well.  This  is 
illustrated  in  Figure  3  below.  Uncertainty  in  the 
chamber  pressure  measurement  (which  was  entered  in 
the  data  reduction  equations  as  the  mean  pressure) 


Figure  3.  Illustration  of  Chamber  Pressure 
Conceptual  Bias. 


were  primarily  due  to  oscillations  and  drops  or  rises 

in  pressure  during  the  bum  duration.  On  average,  this 

conceptual  bias  was  ±30  psi.  Using  the  root-sum- 
square  method  of  combining  the  conceptual  and 
calibration  biases  gave  a  total  bias  in  chamber 
pressure  of  ±3 1.6  psi.  _ 

The  primary  contribution  to  the  bias  in  the  i^ition 
and  shutdown  times  was  a  conceptual  uncertainty  in 
the  definition  of  bum  time.  This  concept  is 
illustrated  in  Figure  4  which  shows  the  pressure-time 
plot  for  a  typical  HTPB  test  conducted  using  an  aft 
mixing  section  with  i  =  4.0"  (see  Figure  1). 
Detailed  in  Figure  4  are  the  two  re^ons  within  which 
it  may  safely  be  assumed  the  ignition  and  shutdown 
points  lie.  These  two  regions  were  from  test  to  test, 
on  average,  0.07  and  0.23  seconds  long, 
consecutively.  According  to  convention,  the  bias 
limit  for  each  region  is  therefore  one-half  the  span,  or 
±0.035  and  ±0.1 15  seconds,  as  shown  in  Table  4. 
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Figure  4.  Illustration  of  Bum  Duration 
Conceptual  Bias. 

Taking  the  values  from  Table  4  and  substituting 
them  into  Equation  (7)  gives  a  bias  limit  for  a  typical 
Gox  test  of  ±0.032  Ibm/in^-sec  (for  this  particular 
test,  Gox  =  0.319  Ibm/in^.sec),  or  approximately 
±10.7%.  Next,  substituting  this  bias  value  and 
precision  limit  for  Gox  form  Table  3  into  Equation 
(5)  gives  an  overall  uncertainty  limit  for  Gox  of 
±0.0324  Ibm/in^-sec,  or  ±10.8%.  The  uncertainty 
values  vary  slightly  from  test  to  test  due  to  the 
dependence  of  the  sensitivity  parameters  to  the  actual 
test  conditions. 

Table  5  shows  the  results  of  similar  calculations 
carried  out  on  Equations  (2),  (3),  and  (4)  for  the  same 
test.  Shown  in  the  last  two  columns  are  the  bias  and 
total  uncertainty  limits  in  terms  of  percent  of  the  data 
reduction  result.  By  examining  plots  of  the  results 
generated  from  the  reduction  of  ^e  LHM  data,  with 
’’error  bars"  included  to  represent  the  uncertainty  in 
the  results,  and  by  examining  and  comparing  Tables  1 
through  5,  several  items  can  be  noted  which  are 
indicative  of  the  behavior  of  hybrid  testing. 
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Table  5.  Detailed  LHM  Uncertainty  Analysis  Results  for  a  Typical  Test, 


INPUT  PARAMETERS 

% 

Bias 

% 

Uncert. 

C, 

119 

Tup 

pvent 

n 

EB 

Pf 

L 

L£lj 

mm 

mm 

Dm 

■a 

Gox(1) 

n 

n 

- 

■ 

- 

- 

±10.7 

±10.8 

r(2) 

* 

- 

- 

- 

- 

aa 

n 

- 

- 

±8.2 

±8.7 

%(3) 

■■ 

>/ 

IH 

- 

BE 

- 

- 

±7.1 

±9.2 

C*(4) 

n 

n 
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Figure  5  shows  the  typical  hybrid  plot  of  f  versus 
Gox  with  the  uncertainty  of  the  results  indicated  by 
the  error  bars.  Also  shown  in  Figure  4  is  a  curve  fit 
taken  from  Sutton  for  HTPB.^  The  dashed  lines 
indicate  5%  error  bars  on  the  curve  fit  which  provided 
reasonable  coverage  of  the  original  data  from  which 
the  fit  was  generated.  One  of  the  indications  of 
Figure  5  is  that  the  majority  of  the  data  generated  by 
the  LHM,  within  the  uncertainty  of  the  results,  falls 
within  the  5%  bands  of  the  curve  fit.  However, 
equally  evident  is  the  fact  that  several  other  fits  could 
be  generated  that  would  equally  suffice  for  the  LHM 
data  due  to  the  high  level  of  uncertainty.  This  is  a 
demonstration  of  the  inaccuracy  in  using  labscale 
hybrid  data  to  generate  regression  rate  curve  fits  such 
as  the  one  in  Figure  5.  Further,  once  a  curve  fit  is 
generated,  the  level  of  uncertainty  associated  with  the 
data  is  lost  without  knowledge  of  the  original  data. 
Current  research  by  Brown,  Coleman,  and  Steele^  at 
The  University  of  Alabama  in  Huntsville  is 
examining  the  methodology  necessary  to  determine 
the  uncertainty  of  such  curve  fits  generated  from  test 
data  of  this  type. 

By  comparing  the  precision  and  bias  limits 
(Tables  3  and  5),  it  can  be  seen  that  the  overall 
uncertainties  of  the  regression  rate  and  flux 
measurements  are  dominated  by  the  bias  limits  which 
indicates  good  repeatability.  In  a  closer  examination 
of  Table  5,  the  primary  contributors  to  the 
uncertainty  in  the  data  points  in  Figure  5  are  the 
measurement  of  the  initial  port  diameter  and  the 
conceptual  bias  of  the  shutdown  transient  time.  This 


Flux  with  Uncertainty  Bars. 


points  to  two  methods  to  reduce  the  uncertainty  in 
this  type  of  plot.  The  first  is  to  reduce  the 
uncertainty  in  the  measurement  of  the  port  diameter 
by  either  increasing  the  size  of  the  fuel  grain  or 
increasing  the  measurement  accuracy.  The  second  is 
to  reduce  the  uncertainty  in  the  bum  times.  This  can 
only  be  accomplished  by  increasing  the  bum  duration 
since  the  bias  itself  is  conceptual  rather  than  due  to 
any  measurement  inaccuracy.  Increasing  the  duration 
reduces  the  sensitivity  of  the  regression  rate  equation 
to  the  conceptual  bias  even  though  the  bias  itself 
remains  the  same. 

As  stated  previously,  it  is  sometimes  useful  to 
compare  the  efficiency  of  combustion  in  the  form  of 
the  characteristic  velocity.  Figure  6  shows  a  plot  of 
C*  versus  %  for  several  HTPB  tests  as  compared  to 
theoretical  C*  values  generated  using  the  NASA  SP- 
273  thermochemistry  code.  The  stochiometric  % 
ratio  for  HTPB/N-100  is  approximately  3.0.  Shown 
in  Figure  6  are  the  results  from  tests  using  1.0”, 
4.0",  and  0.0"  or  butted  aft  closure  (see  Figure  7) 
configurations.  Data  points  that  lie  closer  to  the 
theoretical  curve  would  be  considered  to  have  higher 
C*  efficiencies.  Thus,  on  initial  inspection  of 
Figure  6  the  conclusion  may  be  drawn  that  tests 
conducted  using  the  4.0"  closure  have  higher  C* 
efficiencies  when  compared  to  those  using  the  1.0"  or 
butted  aft  closure.  However,  the  uncertainty  bars  on 
the  C*  results  indicate  that  the  level  of  uncertainty 
does  not  permit  the  resolution  required  to  measure  a 
test-to-test  difference  in  C*  for  these  tests.  Thus,  no 


Figure  6.  Plot  of  C*  versus  %  Ratio  for 
Selected  HTPB  Tests  with  Uncertainty. 
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Figure  7.  LHM  Butted  End  Grain 
Configuration. 

conclusions  may  be  drawn  about 
instability  on  the  C*  efficiency  in  the  LHM  ^ts.  In 
fact,  it  is  quite  possible  that  efficiencies  of  greater 
than  100%  could  be  generated  by  the  LHM  fording 
to  the  uncertainty  bars.  Finally,  the  plot  shows  that 
for  this  type  of  analysis  to  be  worthwhile, 
uncertainties  in  the  results  must  be  very  small,  on  the 
Older  of  ±1  %)  to  gain  the  necessary  resolution. 

Referring  again  to  Tables  3  and  5,  by  compmng 
the  bias  and  precision  limits  for  C*  and  %,  it  is 
evident  that  the  precision  limits  contribute  to  a 
sizable  portion  of  the  overall  uncertainties  in  these 
two  results.  This  is  indicative  of  the  difficulty  of 
repeating  particular  test  conditions  for  a  compan^n 
of  this  type.  Also,  this  comparison  indicates  that 
reduction  of  biases  contributions  in  the  input 
parameters  are  only  partially  effective.  Nevertheless, 
in  examining  Table  5  it  is  seen  that  burn  time  bm 
dominates  the  bias  in  %  and  the  that  C  bias  is 
dominated  by  biases  in  nozzle  throat  diameter  and 
chamber  pressure.  Again,  increasing  the  burn 
duration  would  reduce  the  %  bias. 

Two  methods  exist  for  reducing  the  C*  bias.  The 
first  is  to  reduce  the  bias  uncertainty  in  the 
measurement  of  the  nozzle  throat  diameter.  It  is 
doubtful  that  the  accuracy  could  be  incre^ed  due  to 
the  small  size  of  the  nozzle  throat  and  the  fact  thaUbe 
diameier  can  erode  measurably  during  a  test  which 
would  again  lead  to  a  conceptual  bias.  However,  by 
increasing  the  size  of  the  throat,  which  points  to 
increasing  the  scale  of  the  motor  as  a  whole,  would 
reduce  the  sensitivity  of  the  C*  data  reduction 
equation  to  the  throat  measurement  bias.  The  second 
method  to  reduce  the  C*  bias  is  to  select  a  more 
accurate  method  for  calibrating  the  chamber  pressure 
data.  This  will  have  only  a  limited  effect  since  the 
conceptual  bias  dominates  in  this  reading  due  to  the 
oscillations  that  exist  even  in  "stable”  tesu,  or  those 
with  a  low  level  of  pressure  oscillations.  The 
conceptual  bias  is  even  greater  for  "unstable"  t^ts  in 
which  oscillations  of  up  to  ±75  psi  (or 
approximately  ±17%  of  mean  Pc)  have  been  noted.  • 


Another  possibility  is  to  test  at  higher  chai^er 
pressures  thus  reducing  the  percent  bias  in  Pc-  Tins 
would  only  be  effective  under  the  assumption  that  the 
absolute  bias  in  chamber  pressure  remains  constant. 
There  is  no  indication  either  way  as  to  the  validity  of 
this  assumption. 


The  typical  uncertainties  in  the  measurements  of 
the  fuel  regression  rate  (±8.7%),  oxidizer  fUix 
(±10.8%),  motor  charactenstic  velocity  (±1^5^. 
and  the  motor  oxidizer-to-fuel  mass  ratio  (d5.2%) 
were  found  to  place  limits  on  the  use  of  the  labsc^c 
motor  for  quantitative  investigations.  These  limite 
were  indicated  through  an  examination  of  standard 
clots  for  determining  hybrid  performance,  rirst, 
labscale  regression  rate/oxidizer  flux  data  was  shown 
to  be  of  limited  use  in  the  generation  of  curve  tits. 
This  was  due  primarily  to  a  high  measurement 
uncertainty  in  the  grain  initial  port  diameter  and  a 
large  conceptual  bias  in  the  burn  duration.  Seco^, 
the  low  resolution  in  the  C*  versus  %  ratio  plot 
indicates  the  lack  of  utility  of  labscale  data  in 
determining  small  changes  (<10%)  in  C  efficiency. 
This  low  resolution  was  due,  in  part,  to  the  lack  oi 
repeatability,  indicated  by  the  dominant  precision 
index  for  C*  and  %  ratio;  and  in  part  due  to  the 
domination  by  the  measurement  uncer^nty  of  the 
nozzle  throat  diameter  and  conceptual  bias  in  the 
average  chamber  pressure.  a  ac 

In  general,  three  methods  were  identified  as 
approaches  to  reducing  the  uncertainties  mentioned. 
The  first  was  to  increase  the  scale  of  the  test  motor. 
This  would  reduce  the  uncertainties  in  the 
measurements  of  the  diameters  of  the  fuel  grmn  port 
and  nozzle  throat.  The  second  method  was  to  incre^e 
the  burn  duration,  which  coincidentally  would  require 
an  increase  in  the  motor  size.  The  increase  would 
reduce  the  domination  of  the  burn  duration  concepmm 
bias.  The  last  approach  identified  was  to  reduce  the 
uncertainty  in  the  average  chamber  pressure 
measurement  through  either  an  increase  test 

pressure  or  in  the  accuracy  of  the  cal^ration 
methodology.  However,  both  of  these  were  believed 
to  have  questionable  or  limited  effectiveness. 

A  parametric  study  of  the  changes  suggested  is 
recommended  to  determine  the 
improvement  in  the  accuracy  that  could  be  expe«^. 
Additionally,  this  study  should  indicate  the  optim^ 
motor  configuration  for  obtaining  high  quality  hybnd 
performance  data.  At  minimum,  it  is  cle^  that  all 
test  programs  conducting  hybrid  motor  performMce 
measurements  must  have  some  level  of  uncertmnty 
analysis  conducted  in  order  to  determine  the  quality  ot 
the  data  obtained. 
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clean  with  nonuniform  burning,  wiili  some  pronounced  dips 

and  rises.  ,  .  ^  .u 

Overall,  there  is  a  good  qualitative  correlation  between  the 
predicted  augmented  pressure  curve  and  the  static  firing  data 
of  Ref.  4.  Because  of  nozzle  throat  erosion  later  in  the  ex¬ 
perimental  firing,  the  predicted  peak  pressure  is  somewhat 
higher  than  expected  (no  throat  erosion  modeled).  Similarly, 
for  the  flight  case  with  increasing  forward  acceleration  and 
spin  as  the  firing  proceeds,  the  predicted  pressure  peaks  some¬ 
what  higher  than  expected  relative  to  the  reported  data  of  Ref. 
4.  Only  later  into  the  flight  does  the  normal  acceleration  of  the 
motor  rotation  begin  to  dominate  the  lateral  and  longitudinal 
acceleration  components  in  affecting  the  combustion  process. 
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Nomenclature 

At  =  propellant  bum  surface  area 

a  =  burning  rate  coefficient 

a'  -  burning  rate  coefficient  at  reference  pressure  and 
temperature 

c*  =  propellant  characteristic  exhaust  velocity 
K  =  bum  area  over  throat  area 
Thf  =  mass  flow  rate  of  fuel 

=  maximum  allowable  pressure  sensitivity 
n  =  pressure  exponent 

p  s=  chamber  pressure 

P'  =  chamber  pressure  over  reference  pressure 

r  =  propellant  bum  rate 

r,  =  initial  propellant  temperature 

a  =  temperature  sensitivity  coefficient 

a!  =  temperature  sensitivity  coefficient  at  reference 

pressure 

=  propellant  density 
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Subscripts 

r  =  reference  conditions 

s  =  standard  conditions 

1  =  maximum  limit 

2  =  minimum  limit 

Introduction 

HEN  formulating  gas  generator  propellants,  the  de¬ 
signer  must  achieve  a  specified  mass  flow  schedule. 
One  approach  is  to  use  propellants  with  large  pressure  sensi¬ 
tivity  and  a  variable-area  nozzle.  This  promotes  substantial 
variations  in  mass  flow  while  maintaining  a  reasonable  pres¬ 
sure  range.  The  propellant  temperature  sensitivity  and  com¬ 
bustion  efficiency  are  additional  factors  that  influence  the  de¬ 
sign,  The  objective  of  this  work  is  to  calculate  a  propellant 
solution  space  that  fulfills  the  design  requirements  for  a  vari¬ 
able-flow  gas  generator.  The  scope  includes  1)  developing  a 
design  methodology  that  incorporates  mass  flow  rate,  pressure, 
and  temperature  requirements;  2)  deriving  design  equations  for 
propellants  with  temperature-dependent  pressure  exponents; 
and  3)  calculating  regions  of  compliant  propellant  ballistics  for 
an  example  ducted  rocket  application. 

Approach 

Gas  Generator  Constraints 

The  gas  generator  design  is  assumed  to  be  bounded  by  mass 
flow  rate,  propellant  temperature,  and  chamber  pressure  re¬ 
quirements.  The  mass  flow  rate  range  of  a  ducted  rocket  motor 
can  be  determined  from  anticipated  altitude,  flight  M^h  num¬ 
ber,  and  operational  oxidizer-to-fuel  ratio  considerations: 

<  <  m/2.  The  gas  generator  must  deliver  this  entire  range. 

The  operating  environment  dictates  a  range  of  initial  propeUant 
temperatures:  Tl.i  ^  T  ^  7),2*  The  maximum  chamber  pressure 
is  defined  by  structural/weight  considerations  of  the  missiie. 
The  minimum  pressure  could  be  guided  by  either  the  choked- 
flow  or  propellant  extinguishment:  Pj  <  F  ^  ^2*  It  is  not 
required  that  this  pressure  range  be  spanned,  but  the  pressure 
limits  must  not  be  violated. 

For  steady-state  operation,  the  propellant  burning  rate  range 
is  derived  from  extreme  values  of  the  required  mass  flow  rate 
and  the  propellant  bum  surface  area  by 

^2  ”  tfly  / PpAt  nala 
Ti  —  mfndn^PpAtTOMX 

A  stable  equilibrium  chamber  pressure  also  requires  the  pres¬ 
sure  exponent  must  be  less  than  one. 

Figure  1  shows  how  these  bum  rate,  pressure,  and  pressure- 
exponent  limits  form  a  parallelogram-shaped  region  (A-B- 
C-D)  of  compliant  bum  rate/pressure  combinations.  Initial 
temperature  effects  are  illustrated  by  two  propellant  burn  rate 
curves.  The  maximum  temperature  curve  intersects  point  A, 
while  the  minimum  temperature  curve  intersects  point  C.  Ex¬ 
pressing  these  design  criteria  mathematically,  we  have 

r(F|,  ^  r,  and  r(P2.  TO  s  r2  (3) 


(1) 

(2) 


and  n  <  1. 

Malhcnuitical  Derivaflons 

Propellant  burning  rate  is  assumed  as  a  function  of  pressure 
and  temperature  with 

r  =  a,  exp(a(r,  -  (4) 

Small  variations  in  vrill  result  in  large  changes  in  burning 
rate  at  operational  pressures.  Normalizing  with  a  reference 
pressure  at  which  the  contribution  of  the  term  is  zero, 
yields, 

r  =  a',  exp[o'(r,  -  (5) 
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Fig.  1  Ballistics  limitations  diagram. 

where  a',  =  r(P„  T,J  and  P'  =  PIP,,  making  o'  =  “  +  ^ 
p,.  Normalizing  makes  a'  the  temperature  sensitivity  at  the 
reference  pressure.  Changes  in  /3  then  rotate  each  constant- 
temperature  bum  rate  line  about  the  point  that  it  intersects  the 
reference  pressure. 

Propellant  Design  Relationships 

Equations  to  define  a  burning-rate,  pressure-exponent  solu¬ 
tion  space  are  now  developed  that  satisfy  the  design  ^  con¬ 
straints.  Three  limiting  conditions  are  examined:  condition  I, 
the  minimum  exponent  propellant  sketched  in  Fig.  1;  conditicm 
II.  the  maximum  exponent  propellant  that  overlays  line  A-D 
in  Fig.  1;  and  condition  in.  the  maximum  exponent  propellant 

that  overlays  line  B-C  in  Fig.  1.  .  t 

The  range  of  acceptable  pressure  exponents  for  condition  1 
is  found  by  substituting  Eq.  (5)  into  Eqs.  (3),  yielding. 

^  in  rx  - (6) 
€nP|  inP[ 

>  in  -  inaj_  _  T,,  -  Tu  j  (7) 

“  tn  Pi  in  Pi 

The  reference  burning  rate  for  condition  I  is  found  by  equating 
Eqs.  (6)  and  (7)  and  solving 

€na'.U^A',- BW  -  (8) 

where 

,  €n  Pi  in  r,  -  in  P\  in 
£nPi-inP[ 

tnPUT,^  -  rj  -  inP\{Ta  -  TJ 
inPi-inP[ 
tnP\inPi{Ta-T,,) 

€nPi-inP{ 

Substituting  Eqs.  (8)  into  either  Eq.  (6)  or  (7)  yields  the 
minimum  exponent  as  a  function  of  a  and 

„.|,  =  Ai  + Bla'  +  O  (9) 


€n  Ti  —  fa  n 
tnPi-  in  PI 


Di 

’’  in  Pi  -in  Pi 

..  in  Pi(T,2  -  Tu)  -  in  Pi(T,, 

"  in  Pi  -  in  Pi 

Now,  substituting  into  Eqs.  (6)  and  (7)  gives  equations 
for  the  reference  burning  rate  at  conditions  II  and  HI. 

in  a',  |yn  =  inri-  in  Pi  -  (Ta  -  Tu)(a'  +  in  Pi) 

(10) 

in  ailm  =  in  n  -  n,„a«  in  Pi  -  (J1.2  ~  P/jX®'  +  P  in  Pi) 

(11) 

where  n„„„,  the  maximum  allowable  standard  exponent,  and 
Nmax,  the  overall  maximum  exponent  that  the  designer  aUows, 
are  related  by 


=  A(max  ^  10. 

if 

/3  =  0 

(12a) 

=  MwIAX 

-  Tu), 

if 

/3>0 

(12b) 

=  ^MAX  “ 

-  Tu), 

if 

^<0 

(12c) 

Equations  (8-12)  now  define  the  reference  propellant  bun^g 
rate  a',  and  its  corresponding  reference  exponent  n,  as  functions 
of  a'  and  /3. 

Results  for  an  Example  Application 
An  example  application  is  now  shown  for  a  ducted  rocket 
gas  generator.  Table  1  lists  the  selected  and  derived  parameters 
for  the  application.  An  8.0-in.-diam,  end-burning  gas  generator 
was  assumed  to  operate  over  a  Mach  number  range  of  2.5- 
3.5,  at  altitudes  from  6500  to  65,000  ft.  Other  assunqitions 
include  1)  the  reference  pressure  is  1000  psia,  2)  the  reference 
temperature  is  65^,  and  3)  JVmax  =  10.  The  effect  of  traper- 
atiiw.  sensitivity  parameters  a'  and  p  on  the  resulting  a',  vs  n, 
solution  space  will  now  be  presented  and  discussed  for  this 
application. 

Conventional  Propellants:  Effect  of  a  when  B  -  0 

The  design  relationships  are  first  applied  to  conventional 
propellants  that  have  no  pressure  exponent  sensitivity  (fi  -  0). 
Figure  2a  is  a  propellant  diagram  that  describes  the  re^on  of 
acreptahi"  bum  rate  properties  as  a  function  of  increasing  a' 
(equivalent  to  conventional  propellant  temperature  seimtiviQr 
for  this  case).  The  x  axis  represents  the  propellant  burning  rate 
at  the  reference  pressure  (P,  =  1000  psia)  and  reference  tOT- 
perature  {Tu  =  65”F).  The  y  axis  represents  the  corresponding 
bum  rate  exponent  at  the  reference  temperature 

The  0  =  0  case  shows  the  largest  area  of  ballistic  properties 
that  will  satisfy  the  gas  generator  constraints.  The  solution 
space  is  a  triangular  region  with  vertices  at  condition  I  [defined 
by  Eqs.  (8)  and  (9)],  condition  n  [defined  by  Eqs.  (10)  and 
(12a)],  and  condition  HI  [defined  by  ^s.  (1 1)  and  (12a)].  Con- 
riirinn  I  represents  the  reference  burtting  rate  for  tiie  minimum 
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Fig  2  Propellant  ballistics  diagram;  a)  a  =  0.6-0.00866,  - 

0.000  and  b)  a  =  0.00318,  ^  =  0  to  -0.00289. 


exponent  propellant  that  satisfies  the  constraints.  The  maxi¬ 
mum  exponent  Nmax  and  standard  exponent  n,  are  identical  in 
this  case. 

As  a'  is  increased,  the  triangular  region  becomes  smaUer. 
The  minimum  exponent  increases  until  the  bum  rate  curves 
span  from  line  A-C  to  fine  B-D  in  Fig.  1,  in  which  case 

«i,Ax  =  (/Vmax  -  A\)IBl  @(^  =  0)  (13) 

giving  the  value  of  0.00866  1/*F,  which  is  plotted  as  the  center 
point  of  Fig.  2a. 

Effect  of  P  at  Constant  a' 

Figure  2b  illustrates  the  effects  of  variations  in  the  pressure 
exponent  sensitivity  for  a  fixed  value  of  In  this  case,  a  is 
assumed  constant  at  0.00318  1/T,  while  the  exponent  sensitiv¬ 
ity  ranges  from  0  to  -0.00289  1/®F.  The  variation  in  exponent 
sensitivity  causes  a  reduction  of  the  triangular  solution  space 
until  it  becomes  a  point  at  the  maximum  allowable  value.  The 
minimum  exponent  increases  while  its  corresponding  reference 
burning  rate  decreases. 

Increasing  the  pressure  exponent  sensitivity  then  reduces 
and  translates  the  propellant  solution  space  to  narrower  ex¬ 
ponent  ranges  and  lower  reference  burning  rates.  The  maxi¬ 
mum  value  of  pressure  exponents  also  decreases  as  p  increases 
[Eq.  (12b)]. 

Conclusions 

Propellant  design  relationships  for  variable  exponent  pro¬ 
pellants  have  been  developed  and  applied  to  a  gas  generator 
design  problem.  The  method  defines  a  region  of  acceptable 
propellant  ballistic  properties  as  a  function  of  two  propellant 
temperature  sensitivity  parameters.  The  design  achieves  mass 
flow  rate,  pressure,  pressure  exponent,  and  temperature  range 
constraints.  The  propellant  diagram  approach  gives  the  de¬ 
signer  or  propellant  formulator  a  region  of  compliant  ballistic 
properties  for  a  given  application  instead  of  a  single  design 
point. 
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Introduction 

ICROSCOPIC  chopped  Kevlar  fibers  were  added  to  am¬ 
monium  perchlorate  (AP)  composite  propellants  to  in¬ 
vestigate  their  bum  rate  enhancing  features.  Kevlar  fibers  have 
been  used  in  the  past^  to  increase  the  strength  of  AP  composite 
propellants,  but  as  a  side  effect  it  was  observed  that  the  fibers 
increased  the  burning  rate  of  an  AP/Al  composite  as  much  as 
27%  at  3.5  MPa.  In  the  present  work,  steady  burning  rate 
measurements  and  combustion  photography  were  used  to 
quantify  the  burning  rate  enhancement  and  suggest  a  possible 
explanation  for  the  observed  increases. 

Experimental 

The  steady  burning  rate  was  measured  for  the  scries  of  AP 
composite  propellants  shown  in  Table  1.  TTie  burning  rate  of 
each  propellant  was  measured  using  the  fuse  wire  technique 
in  a  nitrogen-purged  combustion  bomb.  Strands  of  7  X  7  inm 
cross  section  and  30-60  imn  length  were  coated  lightiy  with 
vacuum  grease  as  an  inhibitor  and  were  ignited  by  nichrome 
wire.  The  initial  temperature  of  strands  was  ambient  room  tem¬ 
perature  (20-25X). 

High-speed  and  microscope  photography  were  used  to  in¬ 
vestigate  the  qualitative  differences  between  the  gas  phase 
combustion  and  surface  condition  of  the  various  propellants. 
Conventional  VHS  camcorder  movies  and  3S-inm  SLR  m^- 
rolens  photography  were  also  used  to  reveal  macroscopic  dif¬ 
ferences  between  the  propellants.  Although  only  photogr^hs 
from  the  35-mm  photography  are  presented  here,  some  obser¬ 
vations  based  on  the  other  photographic  techniques  arc  dis¬ 
cussed  when  applicable. 

Results 

Burning  rate  measurements  of  the  nonmetalfized  propellants 
arc  shown  in  Fig.  1  along  with  the  corresponding  burning  rate 
equation: 

r^aP^ 

where  burning  rate  r  is  in  mm/s  and  pressure  P  is  in  MPa.  The 
coefficient  a  and  exponent  n  were  determined  by  a  least- 
squares  fit,  and  a  linear  correlation  coefficient  of  0.98  or  better 
was  calculated  for  each  of  the  curve  fits.  Figure  la  shows  that 
the  addition  of  small  amounts  of  Kevlar  increased  the  burning 
rate  and  lowered  the  burning  rate  exponent  slightly  in  non- 
metallized  AP  systems.  Rgure  lb  compares  propellants  with 
fiber  lengths  of  2  and  5  mm  to  a  fibcrless  AP  composite  pro¬ 
pellant  and  demonstrates  increased  burning  rate  with  increased 
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Abstract 

Aninoniura-nitrate(AN)  based  propellant  are 
candidates  for  fuel  generation  in  a  gas 
hybrid  rockets (GHR).  When  used  as  a  gas 
generator  in  conjunction  with  HNO3.  the 
calculated  specific  impulse  increases  as 
much  as  25  seconds  compared  with  AN  solid 
propellant  rockets.  This  study  examines  the 
combustion  mechanisms  of  metal ized,  AN  gas 
generator  propellants  for  this  application. 
Analysis  of  the  combustion  wave  structure  of 
the  AN  gas  generator  propellants  showed  that 
Granular  Diffusion  Flame(GDF)  theory  can 
explain  the  effects  of  the  pressure  and 
temperature  on  the  burning  rate 
characteristics  of  the  AN  propellants. 
Experimental  data  and  theoretical 
calculations  showed  excellent  correlation 
within  the  pressure  and  temperature  range 
investigated.  Results  of  the  analysis  also 
showed  qualitatively  that  chemical  reaction 
in  the  gas  phase  is  a  dominant  factor  on 
determining  the  burning  rate  of  these 
propellants.  The  gas  phase  sensitivity  of 
the  combustion  wave  dominated  the  propellant 
temperature  sensitivity.  Thus,  it  is 
concluded  that  the  burning  rate 
characteristics  and  temperature  sensitivity 
of  these  AN  propellants  are  controlled  by 
the  gas  phase  chemical  reactions  despite 
existence  of  thick  melt  layer  on  the  burning 
surface. 

Introduction 

Ammonium  nitrate  (AN)  is  an  insensitive, 
affordable  and  Cl -free  oxidizer  crystal  for 
solid  propellants.  Recently,  these 
characteristics  have  become  important  for 
next  generation  solid  propellant  rockets  and 
hence.  AN  has  regained  interest  as  a 
replacement  of  aimonium  perchlorate(AP)[l- 
8].  Compared  with  AP  composite  propellants, 
the  AN  propellants  must  overcome  the 
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problems  of  low  burning  rate,  high  pressure 
exponent,  and  low  specific  impulse(Isp) . 

For  burning  rate  augmentation,  it  is 
known  that  some  chromate  compounds,  such  as 
CrjOj,  CuoCrjO;.  or  (NH4)2Cr207  act  as 
catalysts[2,9,10].  However,  these  additives 
have  not  achieved  the  burning  rate  levels 
that  are  required  for  practical  systems. 

The  high  pressure  exponents  associated 
with  AN  (n-0.8)  are  not  desirable  for  solid 
propellants.  However,  this  characteristic 
is  desirable  for  variable  flow  gas  generator 
propellants  as  some  studies  show[ll-13]. 

Metal  particles  such  as  A1  and/or  Mg  have 
been  introduced  into  AN  composite 
propellants  to  increase  the  specific 
impulse[1.3].  Unfortunately,  all  such 
attempts  were  not  completely  successful 
without  some  AP  to  achieve  the  temperature 
required  for  metal  particle  ignition. 

The  hybrid  rocket  also  has  been  regained 
interest  as  an  insensitive.  Cl -free  and 
affordable  propulsion  system.  Especially, 
the  gas  generator  type  hybrid 
rocket(GHR)[14-22]  combined  with  insensitive 
gas  generators,  which  is  shown  in  Figure  1, 
is  considered  as  a  near-term  application. 
Although  some  studies  focused  on  utilization 
of  high  energy- density  materials  such  as 
glycidyl  azide  polymer(GAP).  experimental 
work[19.20]  revealed  that  hydroxyl 
terminated  polybutadiene(HTPB)  as  a  fuel 
binder  still  produces  a  good  performance,  so 
that  a  fuel -rich,  insensitive  conventional 
gas  generator  is  a  possibility  for  the  GHR 
gas  generator.  For  this  conventional  gas 
generator,  AP  cannot  be  used  without  a  Cl 
scavenging  agent  because  it  contains  large 
amount  of  Cl. 

On  the  other  hand,  problems  of  AN 
composite  propellants  described  above  can  be 
overcome  if  GHR  is  introduced.  Thus, 
ammonium  nitrate(AN)  was  chosen  as  an 
oxidizer  in  this  study. 

It  is  known  that  there  are  several 
reaction  paths  on  the  decomposition  and 
reaction  of  AN[9, 23-26]  which  makes  detailed 
investigation  very  difficult.  Thus, 
detailed  combustion  mechanisms  of  AN 
composite  propellants  have  not  yet  been 
verified. 

In  this  study,  theoretical  performance  of 
GHR  with  AN  composite  propellants  is  briefly 
described.  The  burning  rate  and  temperature 
sensitivity  of  AN  composite  propellants  are 
investigated  in  order  to  gain  basic 


1 

American  Institute  of  Aeronautics  and  Astronautics 


understanding  of  combustion  mechanisms  of 
the  propellants. 

Theoretical  Performance  of  GHR 

GHR  performance  with  AN  composite 
propellants  is  calculated  with  using 
thermochemical  equilibrium  codes [27]. 

Results  are  shown  in  Figure  2.  The  liquid 
oxidizer  chosen  is  HNO3,  which  is  storable 
in  room  temperature.  For  comparison,  the 
Isp  of  AN  composite  propellant  which  refers 
to  solid  propellant  rockets  is  also  shown. 
Equivalence  ratio.  Ogq.  is  equal  to  1.0  for 
both  cases.  Chamber  pressure.  Pg.  is  varied 
from  5  to  20  MPa.  Values  of  Isp  are  the 
ones  at  optimum  expansion. 

As  shown,  the  Isp  of  the  GHR  ranges  from 
247  to  278  seconds,  which  is  greater  than 
that  of  solid  rockets  at  the  same  chamber 
pressure.  For  example,  the  Isp  of  the  GHR 
reaches  approximately  263  seconds  at  Pc  =  10 
MPa.  while  that  of  AN  composite  is 
approximately  238  seconds.  The  GHR  produces 
the  same  Isp  at  lower  pressure,  or  higher 
Isp  at  the  same  pressure.  The  Isp  increase 
shown  here  is  considered  to  be  enough  to 
compensate  weight  increase  caused  by 
additional  oxidizer  tank,  feed  system  and 
secondary  chamber.  Mean  molecular  weight  of 
the  exhaust  gas  is  relatively  constant 
through  chamber  pressure  range,  while 
adiabatic  flame  temperature  increases  as  the 
pressure  increases.  Thus,  the  pressure 
dependency  of  Isp  is  caused  by  the  flame 
temperature  increase. 

Figure  3  shows  relationships  between 
vacuum  specific  impulse,  Ispygc  various 
liquid  oxidizer  /fuel  gas  mass  flow  ratio. 
0/F.  HTPB  concentration  within  the  gas 
generator  propellant  is  changed  from  20  to 
60  twt.  Chamber  pressure  is  lOMPa.  No 
metal  compounds  are  included  in  this 
calculation.  Optimum  0/F  which  gives 
maximum  Isp.ygc  remained  unchanged  except 
for  the  50*  HTPB  propellant.  For  a  system 
stand  point  lower  0/F  ratio  is  suitable  to 
decrease  the  oxidizer  tank  size.  In  this 
point  of  view,  lower  HTPB  concentration  or 
higher  AN  concentration  within  the 
propellant  is  desirable.  Ultimately  the 
optimum  0/F  and  fuel  concentration  have  to 
be  determined  for  each  specific  application 
with  system  trade  studies. 

Burning  Rate  Characteristics 

Typical  burning  rate  of  AN  composite 
propellants  is  shown  in  Figure  4[3,7].  This 
composition  includes  60*  AN.  20*  HTPB  and 


20*  Al.  Initial  temperature  was  changed 
from  243  to  343  K.  Vieille’s  law  (  or  St. 
Robert’s  law)  for  the  propellant  yields 
following  expressions: 

rb  =  a  Pc"  (1) 

where  a  and  n  are  tabulated  as  Table-1.  Tq 
is  initial  temperature  of  the  propellant. 

As  is  clearly  seen,  the  burning  rate 
linearly  increases  as  pressure  increases  in 
In  Pc  vs.  In  ri,  plane.  Pressure  exponent  of 
the  burning  rate,  n,  does  not  change  as  Tq 
changes,  although  coefficient,  a,  decreases 
as  Tq  decreases. 

The  low  burning  rate  of  the  AN  composite 
propellant  is  generally  explained  by 
endothermic  phase  change  from  solid  to 
liquid  which  produces  a  thick  melt  layer  on 
the  burning  surface  of  the  propellants. 

This  melt  layer  is  believed  to  inhibit  heat 
feedback  from  the  gas  phase  to  the  burning 
surface,  thus  decreases  the  burning  rate  of 
the  propellant.  However,  an  analysis  in  the 
following  section  revealed  that  melt  layer 
does  not  play  major  role  on  determination  of 
the  burning  rates. 

Analysis  on  Combustion  Mechanisms 

There  are  many  combustion  models  on  solid 
rocket  propellants,  such  as  Thermal  Layer 
model [23],  GDF  model [28.29],  Hermance’s 
model [30],  BDP  model [31.32],  PEM 
model [33.34],  etc..  Major  effort  has  been 
taken  to  build  a  numerical  code  which  takes 
into  account  many  details  of 
phenomenological  descriptions  especially 
after  Hermance[30] .  In  order  to  understand 
and  explain  the  overall  mechanism  of 
combustion,  however,  the  GDF  model  is  still 
believed  to  be  a  good  starting  point[35]. 

The  GDF  model  first  assumes  steady- state 
combustion.  The  propellant  matrix  consists 
of  oxidizer  and  binder  portion,  and  these 
are  distributed  isotropically  in  the  matrix. 
Physical  properties  of  the  solid  phase  are 
assumed  constant  in  space  and  time. 

It  has  been  observed  that  there  is  a  melt 
layer  on  the  burning  surface  of  AN  composite 
propellants[23,24.36].  In  order  to  verify 
the  phenomena  within  the  melt  layer,  let  % 
be  residence  time  in  the  melt  layer.  Then, 

(Pm^Ps^^*-m^"b^ 

~  l-m^"b 

where  Pn,  and  pg  are  densities  of  melt  layer 
and  solid  phase,  respectively.  L,„  is 
thickness  of  the  melt  layer  and  rjj  is 
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burning  rate.  According  to  experimental 
data  of  Brewster.  et.al.[l],  L^,  and  Cb  are 
50  nm  and  0.5  mm/ sec.  respectively. 

Therefore,  the  residence  time  of  the  melt 
layer  is  in  the  order  of  10"^  sec. 

On  the  other  hand,  a  procedure  similar  to 
the  60F  model  yields 

d2/D,  <3) 

where  d  and  are  oxidizer  diameter  and 
diffusion  coefficient,  respectively,  and 
details  of  the  derivation  is  given  in  the 
following  section.  Brewster,  et.al.  gives  d 
-  6  X  10“^  m.  On,  is  in  the  order  of  10’^ 
m2/sec  [37].  Therefore.  is  in  the  range 
of  10^  sec. 

Equation. (3)  is  very  sensitive  to  d  since 
it  is  proportional  to  d^.  and  actual  d  of 
the  propellant  has  distribution  so  that  the 
residence  time  changes  both  in  time  and 
space  within  the  combustion  wave  structure. 
However,  the  result  of  this  simple  analysis 
indicates  that  the  time  required  for 
complete  diffusion  of  oxidizer  and  fuel 
components  can  be  much  greater  than  its 
residence  time  in  the  melt  layer. 
Consequently,  no  diffusion  or  mixing 
processes  were  considered  within  the  melt 
layer.  No  reactions  between  liquefied  AN 
and  HTPB  are  also  assumed  within  the  melt 
layer  since  decomposition  of  AN  takes  place 
just  above  the  burning  surface[23-25]. 

Hence,  the  burning  surface  of  the 
propellant  is  defined  just  above  the  melt 
layer  so  the  liquefaction  is  treated  as  a 
condensed  phase  reaction.  Above  the  burning 
surface,  fuel  and  oxidizer  gas  pockets 
produced  by  decomposition  of  both  AN  and 
HTPB  are  assumed  to  exist  as  shown  in  Figure 
5.  Heterogeneity  in  the  gas  phase  in  which 
exothermic  reactions  generate  heat  feedback 
to  the  burning  surface  is  taken  into 
account . 

It  is  interesting  to  note  that 
experimental  investigations  done  by 
Whittaker  &  Barham[36],  Kubota,  et.al. [4], 
Alspach[38]  and  Brewster  et.al. [1]  show  an 
apparent  tenperature  step  just  above  the 
burning  surface  as  shown  in  Figure  6. 

Though  other  investigators  suppose  this  as  a 
gas  phase  reaction  zone,  Brewster  et.al. 
explain  this  as  a  "build-up"  of  some 
residual  materials.  In  their  experiments, 
the  thickness  of  the  build-up  layer  is  the 
same  order  of  that  of  melt  layer  if  a 
microthermocouple  junction  is  placed  on  the 
AN  crystal.  The  thickness  of  the  “build-up” 


layer  reduces  to  1/4  when  the  junction  is  on 
the  "fill"  region  which  is  fuel -rich  pocket 
propellants.  If  it  is  a  carbonaceous  build¬ 
up,  then  the  layer  must  be  thicker  when  the 
junction  is  on  the  fuel -rich  region. 
Consequently,  it  is  not  likely  that  the 
build-up  is  carbonaceous  materials. 

Brill  et.al. [25]  showed  that 
decomposition  reaction  of  AN  is  very 
complicated  and  even  may  consist  of  re¬ 
combination  of  anmonia  and  nitric  acid. 

Though  the  first  stage  decomposition  is 
believed  as  decomposition  into  ammonia  and 
nitric  acid  molecules,  their  experiments 
confirmed  complication  of  decomposition  and 
reaction  processes  of  AN.  On  the  other 
hand.  Alspach’s  data  clearly  show  the 
existence  of  a  two  stage  reaction  in  the  gas 
phase,  one  just  above  the  burning  surface 
which  corresponds  to  "build-up"  temperature 
and  the  other  far  from  the  burning  surface. 
Although  data  of  Kubota  et.al [4]  are  not 
conventional  AN  propellants,  these  also  show 
the  same  tendency  in  general.  Therefore,  it 
is  concluded  that  no  matter  how  the  first 
stage  gas  phase  reaction  is  completed,  it 
nevertheless  exists.  Since  AN  consists  of 
anmonia  and  nitric  acid.  NO2  and  N2O  may  be 
considered  as  intermediate  combustion 
products  as  in  case  of  double-base 
propellant  combustion.  Unfortunately, 
detailed  gas  phase  reactions  are  difficult 
to  measure  by  experiments.  Further 
investigations  are  required  to  verify  the 
reaction  processes  of  the  combustion  wave  of 
the  propellant. 

Although  mechanisms  of  these  reaction 
zones  are  still  to  be  investigated,  it  also 
can  be  possible  to  treat  overall  mechanism 
as  a  combination  of  monopropellant  flames  of 
AN  and  diffusion  flames  of  decomposition 
gases  of  AN  and  binder.  As  a  first 
approximation,  the  AN  monopropellant  flame 
is  assumed  to  be  in  the  vicinity  of  the 
burning  surface  so  that  it  is  considered  as 
a  burning  surface  reaction. 

The  assumptions  made  above  enable  us  to 
use  GDF  model  with  collapsed  AP 
monopropellant  flames,  which  is  mainly 
concerned  AP  composite  propellant 
combustion[29].  Thus,  following  procedure 
is  essentially  the  same  as  GDF  model  itself 
although  physical  meanings  of  each  parameter 
and  coefficient  are  different. 

Total  characteristic  length  of  combustion 
wave,  L.  is  assumed  as  a  linear  summation  of 
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that  of  physical  diffusion.  L^j.  and  chemical 
reaction.  L^. 

^  *-c 

Next,  heat  flux  in  vicinity  of  the 
burning  surface  is  considered.  Radiative 
heat  transfer  from  the  gas  phase  to  the 
burning  surface  is  neglected  in  this  case. 
Then. 

(XgCl>L)/l-  =  Ppf’bCp'J^L  (5) 

where  1  is  thermal  conductivity,  p  is 
density  and  c  is  specific  heat.  Subscript  g 
denotes  gas  phase  and  p  denotes  propellant 
liquid  phase.  <I>l  '*'l  defined  as 

Ol  =  Tg  ■  Ts  (6) 

'i'L  =  Ts  -  To  ■  Qg/Cp  (7) 

where  Tg  is  gas  phase  temperature.  Tj  is 
burning  surface  temperature  and  Qg  is  heat 
generated  on  the  burning  surface. 

Therefore,  burning  rate  is  expressed  by 
Eq.(5)  once  L  is  properly  estimated. 

However,  care  has  to  be  taken  that  primary 
flame  which  corresponds  to  AP  monopropellant 
flame  must  be  collapsed  so  as  to  be 
approximated  as  a  surface  reaction  term. 

Estimation  of  L(j  and  Lq  is  obtained  as 
the  same  manner  as  GDF  model .  In  case  of 
L(j  is  dominant  in  the  combustion  wave. 

Ld  “  “g  ^d 

where  Ug  and  t(j  are  gas  velocity  and 
residence  time  in  the  gas  phase, 
respectively.  Simplified  mass  conservation 
for  steady- state  combustion  yields, 

Ug  =  Ppr|5  /  Pp  (9) 

In  order  to  estimate  tj,  mass  and  size  of 
the  decomposed  gas  pocket  is  considered  as 

Pg  d3  =  p  (10) 

where  p,  d.  R  are  mass,  scale  and  gas 
constant  of  the  gas  pocket,  respectively. 

Let  Dg  be  a  diffusion  coefficient.  Then, 
Pick's  1st  law  gives 

J  =  -Og  OC/5y)  (11) 

where  J  is  molar  number  of  diffusive  matter 

per  unit  surface  and  time.  C  is 
concentration  and  y  is  direction  of 
diffusion.  J  is  represented  with  using  p  as 

J  -  (p/M)(d2/Td)  (12) 

where  M  is  molecular  mass.  Spacial 
distribution  of  decomposed  gases  just  above 
the  burning  surface  is  approximated  as  a 


linear  function  of  distance.  Then,  (3C/3y) 
is  expressed  as 

(ac/ay)  -  -(p/M)/d4  (13) 

Substituting  Eqs.  (12)  and  (13)  into  Eq. 

(11).  one  can  obtain  an  estimation 

Td  -  d2/Dg  (14) 

Now,  Eqs.  (9).  (10).  and  (14)  are 
substituted  into  Eq.  (8). 

Ld  =(p2/3/Dg)  (Pprb/Pg3'5)  (15) 

In  order  to  eliminate  r^  in  Eq.  (15).  Eq. 

(5)  is  substituted.  Then  Ld  is  expressed  as 

Ld  =  {Xg®L/(Cp‘FL)}l/2 

X  (p/rg5/6  Dgl/2)  (16) 

Next,  characteristic  length  for  chemical 
reaction  is  estimated.  Using  residence 
time  for  chemical  reaction, 

i-c  ~  (ppr|;)/pp)'rc  (17) 

Overall  reaction  is  supposed  as  a  one-step, 
second  order  reaction  as 

Oxidizer  +  Fuel  ->  Product 

then  the  x^  becomes 

==  (l-e)2pgZgexp(-Eg/RuTg)  (18) 

where  e  is  mass  fraction  of  decomposed 
gases,  Zg  is  pre-exponential  factor.  Eg  is 
activation  energy  and  Ru  is  universal  gas 
constant.  Although  molecular  corrosion 
theory  gives  Zg  as  a  function  of  Tg^'^,  this 
term  is  treated  as  a  constant.  Substituting 
Eq.(18)  into  Eq,(17)  and  taking  limit  e->0. 

l-c  ~  (Pp^'b^Pp^ 

{PgZgeXp(-Eg/RUTg)}-l  (jg) 

Eq.  (19)  is  substituted  into  Eq.  (5). 

Lc  =  {T-g^L  ' 

/CPg{Zgexp(-Eg/RuTg)}l/2]  (20) 

To  obtain  the  GDF  expression  of  the 
burning  rate  characteristics.  Eqs.  (4).  (5). 
(16)  and  (20)  are  combined  and  yield 

l/rp  =  PpC(/g{‘l*L^  Cp4'L)}1^2 

X  {(p/Pg^'V^^^ 

+  l/[PgZgl/2 

X  exp(-Eg/2RuTg)]‘^}  (21) 

Suppose  the  change  of  Tj  as  pressure 
changes  is  sufficiently  small.  Then, 

p//6  Ogl«  -  (P,  Dgll'V" 
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Using  Eq.  (22)  and  perfect  gas  equation  of 
state,  and  approximate  all  terms  other  than 
rj)  and  pg  be  constants,  we  obtain  the  final 
form  of  the  GOF  expression  of  the  burning 
rate  as 

P^/  rtj  =  a*  +  b*  (23) 

where  a*  and  b*  are  constants  represent 
effects  of  chemical  reaction  and  physical 
diffusion  on  the  burning  rate 
characteristics.  Note  that  a*  includes 
reaction  term  which  generates  heat  in  the 
gas  phase,  while  b*  does  not  include  any 
reaction  terms.  As  a  result,  pressure 
exponent  of  the  burning  rate  does  not  become 
less  than  1/3  in  6DF  expression. 

Experimental  data  shown  in  Figure  4  are 
replotted  in  Figure  7  using  Eq.(23). 

Despite  scatter  in  the  data,  it  agrees  well 
with  the  modified  GDF  model  described  above. 
Thus,  the  overall  combustion  mechanism  of  AN 
composite  propellant  is  explained  as  a 
combination  of  physical  diffusion  and 
chemical  reaction. 

As  Eq.  (23)  shows,  burning  rate  of  the 
propellant  is  controlled  by  a*  and  b*  in  the 
GDF  model.  Figure  8  indicates  the 
characteristics  of  these  two  parameters  as 
Tq  changes.  The  order  of  a*  is  10^  through 
Tq  is  changes  from  243  to  343  K  while  that 
of  b*  maintained  10*1  to  10®.  Thus,  it  is 
clear  that  a*  is  the  dominant  parameter  on 
determining  the  overall  burning  rate 
characteristics  of  the  propellant  . 

Temperature  Sensitivity  Analysis 

Equation  (23)  derives  that  temperature 
sensitivity  of  the  propellant  can  be 
decomposed  into  the  sensitivities  of  a*  and 
that  of  b*.  These  parameters  are  defined  as 
follows: 

Op, a*  =  3a*/aTo  (24) 

=  -15.0  X  atm(nin  /  sec  K)'' 
Onb*  =5b*/STo  (25) 

»-i 

»  *0.56  X  atm  (mm  /  sec  K) 

Both  terms  have  sensitivity  to  Tq, 
although,  the  order  of  Op  a*  is  much 
greater  than  that  of  Op  b*  as  shown. 
Therefore,  the  temperature  sensitivity  of  AN 
composite  propellant  is  considered  to  be 
chemical  reaction  controlled  just  as  the 
overall  burning  rate  is  chemical  reaction 
dominant. 


If  the  mechanism  of  temperature 
sensitivity  is  chemical  reaction  dominant, 
then  it  can  be  supposed  that  the  mechanism 
is  similar  to  that  of  double-base 
propellants  in  which  flame  structure  is  pre- 
mixed.  no  diffusion  flame  is  observed  and 
chemical  reactions  determine  the  burning 
rate  of  the  propellants.  Now,  Op  is 
expressed  as [39 -42] 

CTp  =  -H  T  (26) 

=  (Q  +  ©-hM0/2 

where  Oand  SF  are  temperature  sensitivity 
of  the  gas  phase  and  the  condensed  phase, 
respectively.  Then  each  terms  in  Eq.  (27) 
is  defined  and  expressed  as 

Q=  OTg/aTo)p  {Eg/(Ru  Tg2)}  (28) 

0=  {(3Tg/aTo)  -  (5Ts/aTo)}p 


/(Ts-To-Qs/Cp) 

(29) 

4'  =  {l-(  aTs/6To)}p 

/(Ts-To-Qs/Cp) 

(30) 

If  the  burning  rate  is  expressed  as  an 
Arrhenius  type  reaction  at  the  burning 
surface  temperature,  we  obtain 

(cTs/cTo)p  =  Op  Ru  Ts2/Es  (31) 

Estimated  values  of  Tj,  Qj,  £$  ^p 
600  K[23.24].  -1180  kJ/kg[l],  169.3 
kJ/kraol[43]  and  1.90  kJ/kg[43], 
respectively.  By  using  Eqs.  (26)  and  (30), 
results  of  temperature  sensitivity  analysis 
shown  in  Figure  9  was  obtained.  It  is  clear 
that  <b  is  a  dominant  factor  on  determining 
the  overall  temperature  sensitivity  of  the 
propellant  burning  rate. 

Although  low  burning  rate  of  AN  composite 
propellants  is  caused  by  the  melt  layer 
which  acts  as  a  thermal  barrier  for  heat 
flax  from  the  gas  phase,  characteristics 
analysis  of  the  temperature  sensitivity 
described  above  revealed  that  the  gas  phase 
chemical  reaction  is  the  dominant  factor  on 
temperature  sensitivity  determination. 

Thus,  liquefaction  of  AN  does  not  play  an 
important  role  on  temperature  sensitivity  of 
AN  composite  propellants. 

Conclusions 

An  AN  composite  propellant  was  proposed 
as  a  gas  generator  propellant  for  GHR. 
Theoretical  calculations  showed  a  25  seconds 
specific  impulse  increase  when  compared  with 
AN  solid  propellant  rockets.  The  combustion 
mechanism  of  the  propellant  was 
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investigated.  Results  of  the  investigation 
showed  that  the  burning  rate  characteristics 
of  the  propellants  can  be  explained  with  a 
GDF  prediction.  It  was  shown  quantitatively 
that  the  gas  phase  chemical  reactions  are  a 
dominant  factor  on  determining  the  burning 
rate  characteristics.  Analysis  based  on  the 
GDF  model  showed  that  temperature 
sensitivity  of  the  propellant  burning  rate 
is  also  a  chemical  reaction  dominant.  On 
the  other  hand,  analysis  based  on  the 
double-base-like  reaction  controlled 
mechanism  showed  that  temperature 
sensitivity  in  the  gas  phase  is  a  dominant 
factor.  Thus,  mechanisms  of  an  AN  composite 
propellant  combustion  are  explained  as  a  gas 
phase  chemical  reaction  dominant  phenomena. 
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Table-1 

Burning  Rate  Parameters  of  the  AN  Composite  Preopellants 


To(K) 

a 

n 

243 

0.0255 

0.887 

293 

0.0315 

0.887 

_ 

0.0400 

0.887 

Figure  1  Schematic  of  Gas  Hybrid  Rockets 


Figure  2  Performance  Comparison 
of  GHR  and  Solid  Rockets 


Liquid  Oxidizer /Sol id  Fuel 
Figure  3  Specific  Impulse  Characteristics  of  GHR 


Figure  4  Burning  Rate  Characteristics  of  Typical 
AM  Composite  Propellants [3. 7] 
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Physical  Diffusion  Length.  Ld 


Figure  5  Schematic  of  AN  Composite 
Propellant  Combustion 


AN  Conposite  Propellant  Temperature  Profile  at  200  psia 
Diameter  of  the  Thermocouple  Junction  •  150  um 


Distance,  mm 


Figure  6  Typical  Temperature  Profile 

of  AN  Composite  Propel lants[38] 
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6DF  Coefficients,  a*  and  b* 


Figure  8 


Figure  7  GDF  plot  of  the  AN  composite  Propellants 


Chemical  Reaction  and  Physical 
Diffusion  Parameters 


10'^ 


Temperature  Sensitivity. 
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Figure  9 

Mechanisms  of  Temperature 

Sensitivity  of  the  Propellants 
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NUMERICAL  ANALYSIS  OF  A  NONLINEAR  BURN  RATE  EQUATION 
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In  this  paper,  a  nonlinear  solid  rocket  propellant  bum  rate  equation  is  developed.  Temperature  sensitivity 
coefficients  are  produced  based  on  the  nonlinear  equation.  Modeling  of  the  solid  propellant  bum  rate  with  the 
nonlinear  equation  is  compared  with  equations  used  to  model  the  solid  propellant  bum  rate.  The  analysis 
compares  solid  propellant  bum  rate  data  produced  by  the  Beckstead.  Derr.  Price  Multiple  Flame  Model.  This  data 
is  taken  at  nine  pressures  and  three  temperatures.  This  paper  studies  the  effects  of  the  number  of  pressures  used 
and  the  propellant  particle  size.  Also,  a  random  error  is  added  to  the  data  points.  It  is  shown  that  a  higher  order 
version  of  the  nonlinear  bum  rate  equation  with  fewer  data  points  is  a  better  model  than  the  current  models  with 
more  data  points. 


NOMENCLATURE 

a  pressure  coefficient 

c*  characteristic  velocity 
D^  coefficient  of  determination 
K  ratio  of  propellant  burning  surface  area  to  nozzle 
throat  area 

n  pressure  exponent 

P  chamber  pressure 

r  solid  propellant  bum  rate 
Sr  sum  of  the  square  of  the  residuals 
uncertainty  around  the  mean 
Sy/  standard  error  of  the  estimate 
T  temperature 

a  empirical  constant 

P  empirical  constant 

e  residual 

xc  temperature  sensitivity  of  characteristic  velocity 
temperature  sensitivity  of  pressure  at  constant 
motor  geometry 

p  g  density  of  the  combustion  gases 
pp  density  of  the  propellant 

*  Engineer  II.  Thermal  and  Fluids  Analysis  Section.  Member.  AI AA 
t  Asst  Professor,  Dept  of  Mech  and  Aero  Engr.  Member,  AIAA. 

Tliis  paper  is  declared  a  work  of  the  U.S.  Government  and  is  not  subject 
to  copyright  protection  in  the  United  States. 


ap  temperature  sensitivity  of  bum  rate  at  constant 
pressure 

INTRODUCTION 

The  Ducted  Rocket  Engine 

A  ducted  rocket  engine  (DRE)  is  an  engine  that 
uses  air  as  an  oxidizer  and  takes  its  fuel  from  a  gas 
generator.  The  air  and  the  fuel  mix  in  the  combustion 
chamber  and  are  ejected  through  a  nozzle.  A  DRE 
requires  a  flight  speed  of  about  Mach  2  to  operate 
efficiently.  To  get  a  ducted  rocket  engine  up  to  this 
operating  speed,  a  solid  propellant  boost  motor  is  used. 
The  boost  motor  casing  serves  as  the  combustion 
chamber  for  the  air  and  fuel  mixture.  Figure  1  shows  a 
DRE  with  a  solid  boost  motor  attached. 

Accurately  describing  the  bum  rate  of  the 
propellant  is  important  when  calculating  the 
performance  of  the  boost  motor.  If  the  burning  rate 
relationship  used  to  calculate  performance  varies 
significantly  from  the  actual  ballistic  perfommce  of 
the  boost  motor  the  performance  will  not  be  accurately 
predicted.  For  example,  if  the  temperature  of  the  boost 
motor  increases,  the  bum  rate  increases.  The  effect  of 
the  initial  propellant  temperature  as  well  as  chamber 
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pressure  influence  the  propellant  burning  rate,  and 
hence  the  motor  performance,  and  are  therefore, 
important  design  considerations  for  solid  boost  motors. 


Fqd  Rkk  Propdiast  Air  lataka  Noafa 


Boost  Motor 


Figure  1  A  ducted  rocket  engine  with  a  solid  rocket 
propellant  boost  motor. 

.Solid  Propellant  Bum  Rate 

The  most  basic  form  of  the  solid  propellant  bum 
rate  expression,  called  St.  Robert's  Law”  is 

r  =  aP" 

In  this  equation,  the  pressure  coefficient  and  exponent 
are  assumed  constants  so  the  bum  rate  is  a  function  of 
chamber  pressure  only.  The  pressure  coefficient  can  be 
a  function  of  the  initial  temperature  of  the  propellant 
and  can  be  approximated  as 

a  =  aoe.\p[ap(T  -  To) 

giving  the  bum  rate  expression  the  final  form  of 
r  =  aoe.xp[CTp(T-To)]P 

Since  the  pressure  exponent  is  assumed  constant, 
and  the  pressure  coefficient  is  constant  for  a  given 
temperature,  the  family  of  bum  rate  curves  are  straight 
and  parallel  for  different  initial  propellant 
temperatures  when  plotted  in  the  log(P)  log  (r)  plane. 
In  reality  the  actual  bum  rate  lines  can  take  a 
nonlinear  form.  Figure  2  demonstrates  how  the  bum 


rate  data  can  take  a  nonlinear  shape  in  a  log-log  plot  of 
bum  rate  versus  chamber  pressure.  The  dashed  line 
represents  a  curve  fit  of  the  data  if  St.  Robert  s  Law  is 
used  as  the  model  equation. 


Figure  2  Propellant  bum  rate  characteristics. 

From  the  literature'*,  the  effects  of  making  the 
bum  rate  pressure  exponent  a  function  of  temperature 
can  be  predicted.  By  making  the  pressure  exponent  a 
linear  function  of  temperature  the  bum  rate  lines  are 
no  longer  parallel.  However,  they  are  still  linear 
because  the  pressure  coefficient  and  pressure  exponent 
are  solely  a  function  of  temperature. 

The  effect  of  the  initial  propellant  temperature  on 
motor  performance  is  shown  in  Figure  3  where 
chamber  pressure  is  shown  as  a  fiinction  of  time  at  two 
different  temperatures.  At  the  higher  temperature  the 
propellant  bums  faster  which  results  in  higher  motor 
pressures.  Correspondingly  at  the  lower  temperature 
the  propellant  will  take  longer  to  bum  and  will  bum  at 
a  lower  pressure. 

Since  the  linear  expressions  do  not  correlate 
exactly  with  the  propellant  bum  rate  curves,  error  will 
naturally  be  introduced  when  using  this  expression  to 
predict  motor  chamber  pressure.  It  is  proposed  that  a 
nonlinear  solid  propellant  bum  rate  expression  be 
established  and  evaluated.  The  nonlinear  expression 
will  describe  the  propellant  bum  rate  as  a  function  of 
pressure  and  temperature.  Evaluations  will  be 
completed  to  see  if  the  bum  rate  is  more  accurately 
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portrayed  and  can  lead  to  better  prediction  of  the 
performance  of  the  ducted  rocket  boost  phase.  Care 
must  be  taken,  however,  that  the  new  e.xpressions  are 
not  so  -nexible"  that  they  fit  "noise"  rather  than  the 
true  burning  rate  curves  of  the  propellant. 


Figure  3  Effects  of  the  initial  propellant  temperature 
on  the  solid  propellant  bum  rate.“ 

First  the  new  burning  rate  e.xpression  will  be 
proposed.  Predictions  calculated  with  the  Beckstead 
Derr  Price  combustion  model”  will  be  applied  to  the 
new  e.xpression.  The  effect  of  the  order  of  the 
nonlinear  e.\pression  on  the  least-squares  residual  will 
be  calculated  as  a  function  of  bum  rate  measurement 
uncertainty,  number  of  pressure  intervals  in  the 
burning  rate  data,  and  will  be  applied  to  propellant 
bum  rate.  Finally,  the  results  will  be  presented  and 
evaluated  to  determine  if  any  improvements  have  been 
realized  in  solid  propellant  bum  rate  prediction. 

A. REVIEW  OFTEMPFRATIJRF.  .«;FN.<;mvTTV 
EXPRESSIONS  FOR  SOLID  PROPFT.r  AMTg 

Temperature  Sensitivity  E.xpressions 

There  are  two  coefficients  commonly  used  in  solid 
propellant  temperature  sensitivity  theory.  The  first  is 
the  temperature  sensitivity  of  burning  rate  at  a  constant 
chamber  pressure*  and  is  defined  as 


This  is  the  percent  change  of  bum  rate  per  degree 
change  in  temperature  and  is  a  propellant  parameter. 
Here  the  data  are  taken  at  constant  pressure. 
Propellant  strand  burner  test  data  are  used  to  determine 
the  bum  rate  to  calculate  the  temperature  sensitivity  of 
bum  rate. 

The  second  coefficient  is  the  temperature 
sensitivity  of  pressure  and  is  defined  in  Equation  5.* 


This  is  the  percent  change  of  the  log  pressure  per 
degree  change  in  initial  propellant  temperature.  A 
motor  firing  is  used  to  determine  the  temperature 
sensitivity  of  pressure  since  it  is  defined  for  a  constant 
motor  geometry  which  is  defined  as  the  ratio  of 

propellant  burning  surface  area  to  the  nozzle  throat, 
area. 

A  third  temperature  related  coefficient  has  also 
been  defined.  It  is  the  temperature  sensitivity  of  the 
characterisUc  velocity  for  a  constant  motor  geometry”. 


To  make  use  of  strand  burner  rate  data  for  motor 
pressure  predictions,  a  relationship  between  the 
temperature  sensitivity  of  bum  rate,  op,  and  the 
temperature  sensiUvity  of  pressure,  Tt^,  coefficients 
must  be  described.  This  could  allow  early  evaluation 
of  candidate  propellant  formulations  and  potentially 
reduce  testing  costs.  Several  expressions  have  been 
proposed  by  the  researchers.  A  summary  of  the  ideas 
previously  presented  is  outlined  in  Table  1. 

For  this  work,  the  Click  and  Brooks  term  will  be 
applied  in  the  development  of  the  nonlinear  equation. 

development  of  THE  NONLINEAR  BURN 
RATE  EXPRES.STON 

To  better  describe  actual  propellant  burning  rate 
data,  a  nonlinear  bum  rate  expression  will  be 
developed  and  evaluated  to  enhance  the  performance 
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prediction  of  solid  propellant  rocket  motors.  This  will 
maximize  the  use  of  strand  burner  tests  allowing  for  a 
close  approximation  of  the  actual  bum  rate  curves.  To 
obtain  a  curved  fit  through  the  given  data  a 
least-squares  multiple  polynomial  regression  will  be 
used  to  describe  the  bum  rate  as  a  function  of  both 
pressure  and  temperature. 


Table  I  Summary  of  Temperature  Sensitivity 
Coefficients 


Author 

Geckler 

and 

Sprenger' 

(d\nr] 

\  av  Jp 

Click  and 

9 

Br(x>ks 

(dint'] 

\  ffV  )p 

Gaunce 

and 

Osborne 

f  dlnr'j 

V  ax  Vp 

H\ 

^Pp”Pg) 

L 

r  an  ^ 

Cohen 

and 

Flanigan^ 

ap^+op^ 

ap 

[*72  **2^^ 

2 

h 

(  'I 

•I2  *21  *^2 

^  »12*11  p"ir"2i 

Komai** 

and 

Frederick 

a  +  pin 

1 

(a+plnP  + Jic) 

The  Nonlinear  Bum  Rate  Equation 

A  nonlinear  bum  rate  equation  is  proposed  in  the 
form  of 

ln(r)  =  zjL,  {aj[In(P)r‘}  (7) 

where  k  will  be  denoted  as  the  order  of  the  equation 
and 

aj  =ajs+aj(T-Ts  W 

The  bum  rate,  chamber  pressure  and  initial  motor 
temperature  are  assumed  to  be  given  by  way  of  strand 
burner  tests.  The  remaining  constants  (aj  g;  aj)  will  be 
determined  using  the  least-squares  multiple  polynomial 
regression. 


I  past-Sniiares  Multiple  Polynomial  Regression 

An  example  of  the  multiple  polynomial  regression 
will  be  developed  using  k  =  3.  Thus.  Equation  7  will 
take  the  form  of 

In  (r)  =  a ,  ,s[ln  (P)]^  +  a2.s[ln  (P)]  ‘  +  aj  ^[In  (P)]^ 

-Hxi  AT[ln(P)]®  +a2AT[ln(P)]‘  +a3AT[ln(P)]2  (9) 

which  reduces  to 

In  (r)  =  a  1  +  32,sln  (P)  +  a3^snn  (P)]^  +  a  i  AT 

-HX2ATln(P)  +  a3ATtln(P)]2  (10) 

To  apply  the  multiple  polynomial  regression,  the 
equation  is  rewritten  with  the  following  substitutions. 

ln(r)  =  z,  ln(P)  =  xi  andAT  =  .X2  (H)*- 

Also.  an  error  term,  epsilon,  is  added  to  the  right  hand 
side  so  that  Equation  10  will  now  describe  the 
difference  in  the  curve  fit  and  the  value  predicted  by 
the  equation. 

'1  2 

z  =  a+bxi  +cxj  +d.X2  +ex2Xi  +fx2Xj  +e  (12) 

In  this  equation,  epsilon  is  the  residual  which  is  the 
difference  between  the  known  bum  rate  and  that 
predicted  by  die  right  hand  side  of  the  nonlinear 
equation.  This  equation  is  solved  for  the  residual. 


The  sum  of  the  square  of  the  error  is  determined  by 
summing  the  right  hand  side  over  the  range  of  known 
points. 

Sr  =  2™i(Zj  -a-bxi  -cx^  -d.X2  -ex2Xi  -fx2  (14) 

The  derivative  of  Equation  14  with  respect  to  each 
unknown  coefficient  is  taken  and  the  results  are  set 
equal  to  zero.  The  solution  to  this  set  of  equations 
represents  a  minimum  sum  of  the  squares  of  the  error. 
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^  =  -2Z“iX^(Zj-a-bx,-cx^-dx2-ex2X,  +el  (xjx^)  +  fl  (xjx^  =x(zx,X2)  (16e) 

-fx2.x^)j=0  (15c)  aZ(x^X2) +bE(xjX2) +cz(xjX2) +dx(x^X2) 

^  =  -2Ej^,X2(Zj-a-bx,-cx5-dx2-ex2Xi  -k!e(xJx2) +fE(x^x^)  =e(zx2x2)  (160 


-fx2xj)j=0  (15d) 

^  =  -2  Ej^l  X2X ,  (Zj  -  a  -  bx,  -  cx^  -  dx2  -  ex2Xi 
-fx2X^)j=0  (15e) 

^  <^rn  2/  i.  2  j 

-^  =  -2Zj^l  X2X,(Zj-a-bx,  -cx,  -dx2-ex2X, 

-fx2x2)j=0  (150 


The  constant  outside  of  the  summation  can  be  divided 
out.  The  summations  are  multiplied  through  and  the 
equations  are  rearranged  to  give  the  following 
equations. 


an  +  bZx| 


+cZX|  +dEx2  +eE(xjX2) 


To  solve  these  equations  for  the  constants,  they  are 
arranged  into  a  6  x  6  matrix  with  a  column  matrix  for 
the  coefficients  and  a  column  matrix  for  the  right-hand 
side. 

N  £x|  Sx;  £(X|X2)  ^(xiX]) 

Sx,  ExJ  ExJ  Sc-xixj)  e(x{x2)  I:(xJx2) 

Exf  Ex’  ExJ  e(xJx2)  e(x’x2)  e(x}xj) 

Ex2  E(X|X2»  t:(x|X2)  Ex|  J^(x|X^)  Z(xjx^) 

E(x,X2)  i(x?X2)  Sfxjxj)  2:(x|xj)  E(x?xj)  e(xJx^) 

e(xJx2)  e(x’x2)  e(xJx2)  e(xJx2)  E(xjxj)  e(xJx^) 


+f  Z  (x^X2)  =  Z  z  ( 16a) 

aZxi  +bZxj  +cZx]  +dz(xiX2)  +ez(x^X2) 
+fz(xjX2)  =z(zxj)  (16b) 


To  solve  this  equation,  the  6  x  6  matrix  was  inverted 
and  multiplied  by  the  right-hand  side  column  matrix, 
giving  the  column  matrix  of  constants.  For  this  paper, 
a  spreadsheet  application  was  used  to  perform  these 
operations  efficiently. 

In  a  similar  way,  St.  Robert's  Law  and  the 
Komai-Frederick  equation  can  be  adapted  for  use.  For 
St.  Robert's  Law.  start  with  the  original  equation. 
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r  =  aoexp[op(T-To)]P”  (1^) 

Following  the  same  derivation  technique  used  earlier, 
take  the  natural  log  of  St.  Robert's  Law.  find  the  sum 
of  the  squares  of  the  residuals,  take  the  derivatives  with 
respect  to  the  constants,  set  the  equations  up  into 
matrix  form  and  solve  them.  The  resulting  matrix  is 
equation  is 


N  E(T-To) 

ZP 

S(T-To)  S(T-To)2 

ZP(T-To) 

X  ’ 

EP  ZP(T-To) 

ZP2 

a|=ln(ao)  82  =  n  (23a.b) 

aj  =CTp  012  =  P  (23c,d) 

x,=ln(P)  .X2  =  (T-To)  (23e.O 

Thus,  the  two  equations  are  mathematically  equal. 
This  leads  to  the  similarities  in  sum  of  the  squares  of 
the  residuals  as  well  as  other  results,  as  will  be  seen 
later  in  the  paper. 

Temperature  Sensitivity  Coefficients  Development 
When  the  nonlinear  bum  rate  expression  is  given 
as  Equation  7  the  temperature  sensitivity  coefficients 
can  be  defined  for  use  with  the  bum  rate  data  available. 


Eln(r) 

Z(T-To)ln(r) 

ZPln(r) 


(19) 


In  a  similar  manner,  the  matrix  equation  for  the 
Komai-Frederick  equation  can  be  derived  and  is  given 
to  be 


Temperature  Sensitivity  of  Bum  Rate  at  Constant 
Pressure 

The  temperature  sensitivity  of  bum  rate  at 
constant  pressure  is  defined  as 


N  I(T-To)  ZP  SP(T-To) 
2(T-To)  Z(T-To)2  SP(T-T«)  ZPCT-To) 
ZP  IP(T-To)  ZP^  ZP(T-To) 
IP(T-T„)  SP(T-To)2  ZP^(T-T„)  Zp2(T-To)^ 


When  Equation  7  is  applied  to  Equation  24  the  result 
is  simply 


ao 

CTp 

Zlnr 

Z(T-T„)lnr 

(20) 

(25) 

n 

IPlnr 

^  P 

ZP(T-To)lnr  1 

which  reduces  to 

Note  that  the  Komai-Frederick  equation  and  the  1” 
order  nonlinear  equation  produce  the  same  results. 
The  sum  of  the  squares  of  the  residuals  equation  for  the 
Komai-Frederick  equation  is  listed  in  Equation  21 
while  the  1“  order  nonlinear  equation  is  Equation  22. 


CTp=ao+ailnP  +  ...  +  aj(lnP)l  (26) 

or  in  summation  form 

<jp  =  zjLo[aj(>nP)j]  (27) 


Sr  =  Ij"  1  [zj  -  In  (ao)  -  <Jp(T  -  To)  -  n  In  (P)] 
[-p(T-To)In(P)]j  (21) 

Sr  =  Z™i  (zj-ai-a2X,  -aiX2 -a2XiX2)  .  (22) 

A  visual  inspection  of  these  two  equations  reveals  that 


Temperature  Sensitivity  of  Pressure  at  a  Constant 
Motor  Geometry 

The  temperature  sensitivity  of  pressure  at  a 
constant  motor  geometry  is  defined  as 


JiK  = 


fainP) 
V  ffT 


(28) 
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Adopting  a  formula  that  Click  and  Brooks  use  to 
calculate  the  temperature  sensitivity  of  pressure  gives 

"K  =  -ii;7(<^P+"c)  (29) 

where  the  hum  rate  pressure  exponent  is  defined  as 

This  leads  to  the  following  expression  for  the 
temperature  sensitivity  of  pressure  at  a  constant  motor 
geometry  in  terms  of  the  nonlinear  equation. 

KK  =  IjL,  {jaj[ln(P)r'}  (31) 

Summary  of  Relevant  Equations 

The  nonlinear  solid  propellant  burning  rate 
equation  and  the  temperature  sensitivity  equations 
have  been  developed  using  a  multiple  polynomial 
regression.  The  equations  are 

•  Bum  Rate,  a  combination  of  Equations  7  and  8 

ln(r)  =  Zji,  {[aj  s  +aj(T-Ts)][ln(P>i-']}  (32) 

•  Temperature  Sensitivity  of  Bum  Rate  at  Constant 
Pressure 

ap  =  ljL,[ajln(P>)-']  (27) 

•  Temperature  Sensitivity  of  Pressure  for  a  Constant 
Motor  Geometry 

’'K  =  Sj=lljajnn(P)r’}  (31) 

It  is  apparent  that  these  equations  can  be  solved  with 
sufficient  burning  rate  data  at  various  pressures  and 
initial  temperatures.  All  of  these  parameters  are  easily 
obtained  from  such  tests. 

APPROACH 

Effect  to  be  Investigated 


The  Beckstead  Derr  Price  (BDP)  multiple  flame 
model  was  used  to  predict  the  bum  rate  of  a  solid 
rocket  propellant.  The  data  consisted  of  points  at  nine 
pressures  and  three  temperatures.  The  data  points 
were  determined  for  four  different  propellant  particle 
sizes. 

The  bum  rate  data  will  have  an  uncertainty  factor 
added  to  it  to  find  out  how  well  the  nonlinear  equation 
deals  with  more  realistic  data.  The  number  of  pressure 
intervals  will  be  changed  from  nine  to  five  to  see  if 
fewer  data  points  work  as  well  as  the  nine  when 
modeling  the  BDP  data.  E.xpressions  for  ap  and 
based  on  the  nonlinear  equation  will  be  developed. 
These  tests  will  help  determine  if  a  nonlinear  bum  rate 
equation  is  more  accurate  and  just  as  practical  to  use  as 
the  other  equations  being  used  now. 

Beckstead  Den  Price  Multiple  Flame  Model  "Data" 

The  BDP  model  was  used  to  predict  the  burning. . 
characteristics  of  a  composite  solid  propellant.  A 
composite  propellant  mixes  an  oxidizer  and  a  fuel 
binder.  The  BDP  model  assumes  the  combustion  of 
the  o.xidizer  crystal  from  the  composite  propellant  is 
made  up  of  three  flames:  the  primary  flame,  the 
premixed  oxidizer  flame  and  the  final  difhision  flame. 
The  bum  rate  is  determined  by  solving  the  heat 
balances  at  the  propellant  surface. 

The  propellant  from  the  BDP  example  was  an 
AP/HTPB  propellant.  The  oxidizer  weight  percent  was 
84.4%  and  the  oxidizer  density  was  1.95  g/cm^  The 
binder  density  was  0.90  g/cm’  giving  the  propellant  an 
overall  density  of  1.65  g/cm\  The  adiabatic  flame 
temperature  was  2968  K  and  the  molecular  weight 
equaled  20.78  g/gmole. 

The  data  collected  was  for  a  series  of  chamber 
pressures  ranging  from  0. 100  MPa  to  40.0  MPa.  The 
temperatures  modeled  were  276  K,  293  K  and  3 10  K. 
Data  was  also  collected  based  on  the  size  of  the 
oxidizer  particle.  This  varied  from  l.O  micron  to 
1000.0  microns.  Table  2  defines  the  data  for  the 
example  case. 

Equation  Analysis 

The  uncertainty  around  the  mean,  sum  of  the 
squares  of  the  residuals,  the  standard  error  of  the 
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estimate  and  the  coefficient  of  determination  were 
calculated  for  the  various  equations  under  different 
conditions.  These  terms  are  explained  in  the  following 
sections. 

Uncertainty  around  the  Mean 

The  uncertainty  around  the  mean  is  a  calculation 
based  on  the  test  or  BDP  data.  The  mean  of  the  data  is 
simplv  the  sum  of  the  bum  rates  divided  b\’  the  number 
of  bum  rate  values.  The  uncertainty  around  the  mean 
is  the  square  of  the  difference  between  a  specific  bum 
rate  value  and  the  mean  bum  rate  value.  The 
important  value  when  analyzing  the  nonlinear  bum 
rate  equation  is  the  total  sum  of  the  squares  around  the 
mean.  In  equation  form  it  is 


Table  2  BDP  Bum  Rate  Data  for  a  Mean  Diameter  of 
1000.0  Microns 


310  K 

293  K 

276  K 

Pressure 

(MPa) 

Bum  Rate 
(in/aec) 

Pressure 

(MPa) 

Bum  Rate 
(in/sec) 

Pressure 

(MPa) 

Bum  Rate 
(in/scc) 

O.l 

0.023 

0.! 

0.022 

0.1 

0.022 

0.211 

0.029 

0.211 

0.029 

0.211 

0.028 

0.447 

0.036 

0.447 

0.036 

0.447 

0.035 

0.946 

0.053 

0.946 

0.052 

0.946 

0.051 

2 

0.089 

2 

0.087 

2 

0.085 

4.23! 

0.149 

4,23 

0.146 

4.23 

0.143 

8.944 

0.242 

8.944 

0.237 

8.944 

0.232 

18.915 

0.376 

18.915 

0.368 

18.915 

0.361 

40 

0.588 

40 

0.575 

40 

0.563 

Slim  of  the  Squares  of  the  Residuals 
To  get  the  "best"  values  of  the  equation 
coefficients,  the  sum  of  the  squares  of  the  residuals  was 
set  up  in  Equation  14.  After  the  coefficients  are 
determined,  they  are  substituted  back  into  Equation  14. 
The  lower  the  sum  of  the  squares  of  the  residuals  is. 
the  better  the  answer. 

Standard  Error  of  the  Estimate 
The  standard  error  of  the  estimate  is  calculated  by 
taking  the  square  root  of  the  sum  of  the  squares  of  the 
residuals  divided  by  the  difference  between  the  number 
of  data  points  and  the  number  of  coefficients  to  be 


solved  for.  The  standard  error  of  the  estimate  in 
equation  form  is 

=  J  N-(m+l 

The  subscript  y/x  denotes  an  error  predicted  for  a  y 
value  that  corresponds  to  a  specific  x  value.  The 
standard  error  of  the  estimate  determines  the  quantity 
of  the  spread  of  the  data  around  the  regression  line. 
This  differs  from  the  standard  deviation  which 
quantifies  the  spread  of  the  data  around  the  mean. 
Once  again  it  is  seen  that  the  higher  the  order  of  the 
equation,  the  better  the  results  are. 

rnefficient  of  Determination 

The  coefficient  of  determination  is  a  way  to  show 
the  improvement  of  one  form  of  the  bum  rate  equation 
over.  This  is  accomplished  by  comparing  the  sum  of 
the  squares  of  the  residual  of  the  various  equations  to^ 
the  total  sum  of  the  squares  around  the  mean.  The 
equation  for  the  coefficient  of  determination  is 

RESULTS 

St.  Robert's  Law.  the  Komai-Frederick  equation 
and  the  nonlinear  equation  at  1“*,  and  3"*  orders 
were  compared.  Areas  of  interest  were  (1)  equation 
type  effects.  (2)  bum  rate  uncertainty  effects  and  (3) 
pressure  interval  effects.  Since  the  Komai-Frederick 
equation  and  the  I"  order  nonlinear  equation  were 
shown  to  be  mathematically  equal,  only  the  1"  order 
equation  will  be  presented  in  the  results. 

Equation  Type  Effects 

In  a  previous  section,  the  derivation  of  the 
nonlinear  equation  for  use  with  the  multiple 
polynomial  regression  method  was  presented.  In  a 
similar  way.  St.  Robert's  Law  and  the  Komai-Frederick 
equation  can  be  adapted  for  use. 

Recall  that  for  the  form  of  the  equations  being 
used,  only  chamber  pressure,  initial  temperature  and 
bum  rate  are  required.  The  standard  temperature  may 
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be  any  temperature.  For  this  study,  the  standard 
temperature  was  one  of  the  temperatures  from  the  BDP 
data.  Specifically,  the  standard  temperature  was  273 
K. 

Consider  the  BDP  data  for  an  average  diameter  of 
1000  microns.  This  data  was  shown  back  in  Table  2. 
The  data  used  in  the  evaluation  process  was  made  up 
of  all  nine  pressures  and  two  different  temperatures. 
The  temperatures  were  276  K  and  310  K.  The  data 
was  evaluated  by  St.  Robert’s  Law.  the 
Komai-Frederick  equation  and  the  nonlinear  equation 
up  to  the  3"'*  order. 

The  uncertainty  around  the  mean  was  defined 
earlier.  It's  value,  which  is  based  on  the  BDP  data, 
determines  the  spread  of  the  data  around  the  mean  and 
is  used  to  measure  improvement  in  modeling  the  data. 
The  sum  of  the  squares  of  the  residuals  was  defined  in 
a  previous  section.  Recall  that  the  lower  the  value,  the 
better  the  fit.  The  standard  error  of  the  estimate,  as 
defined  earlier,  is  used  to  determine  the  quantity  of  the 
spread  of  the  data  around  the  regression  line.  The 
coefficient  of  determination  is  a  measure  of  the 
improvement  in  modeling  based  on  the  mean  of  the 
data.  Comparing  the  different  values  of  this  number 
for  the  different  equation  types  will  show  whether  one 
equation  is  better  than  another  at  modeling  the  data. 
All  of  these  values  are  summarized  in  Table  3  for  the 
nine  pressures  and  two  temperatures  used  in  this  study. 


Table  3  Summary  of  the  Equation  Analysis  for  the 
Nine  Pressures  and  Two  Temperatures 


Uncertainty 
Around  the 
Mean 

Sum  of  the 
Squares  of  tlie 
Residuals 

Standard 
Error  of  the 
Estimate 

Coefficient  of 
Determination 

0.5792 

0.2833 

0.1374 

0.5108 

0.5792 

0.2869 

0.1432 

0.5047 

2"' Oder 

0.5792! 

0.0957! 

0.0893 

0.8348 

3"*  Order 

0.5792 

0.0153 

0.0391 

0.9736 

This  shows  that  the  traditional  St.  Robert’s  Law 
accounts  for  51%  of  the  error  associated  with  the 
uncertainty  around  the  mean  while  the  3’’*  order 
nonlinear  equation  accounts  for  97%.  Based  on  this, 
the  3"*  order  nonlinear  equation  is  a  better  model  of  the 


solid  rocket  propellant  bum  rate.  The  coefficients 
acquired  during  the  evaluation  can  be  used  with  the 
data  for  a  temperature  of  293  K. 

Uncertainty  and  the  BDP  Data 

What  happens  when  random  noise  is  added  to  the 
BDP  data?  It  is  e.xpected  that  results  similar  to  those 
of  the  previous  section  would  be  obtained.  To  check 
the  equations  against  data  with  more  fluctuation  in  the 
curve,  a  ten  percent  noise  factor  was  added  to  the  BDP 
data.  The  noise  factor  could  cause  a  bum  rate  value  to 
increase  or  decrease.  The  object  was  to  make  it  more 
difficult  to  fit  the  curves.  The  equation  analysis 
summary  table  for  the  noise  factor  is  Table  4. 

It  can  be  seen  that  the  3"*  order  nonlinear  equation 
still  provided  the  best  explanation  of  the  uncertainty 
around  the  mean  error.  Note,  however,  that  only  the 
2"**  order  equation  had  a  coefficient  of  determinatioa . 
that  was  relatively  unchanged  from  the  regular  BDP 
data  set  calculations.  Recall  that  in  the  regular  BDP 
data  set  St.  Robert's  Law  provided  a  better  fit  than  the 
1**  order  nonlinear  equation.  It  is  observed,  however, 
that  the  I**  order  equation  provided  a  better  fit  for  the 
data  set  with  the  ten  percent  noise  factor. 


Table  4  Summary  of  the  Equation  Analysis  for  the 
Ten  Percent  Noise  Factor  at  Nine  Pressures  and  Two 
Temperatures  _ 


Equation 

Type 

Uncertainty 
Around  the 
Mean 

Sum  of  the 
Squares  of  the 
Residuals 

Standard 
Error  of  the 
Estimate 

Coefficient  of 
Determination 

St 

Robert's 

Law 

0.6276 

0.3666 

0.1563 

0.4159 

l-Ordcr 

0.6276! 

0.3617 

0.1607 

0.4237 

2*' Order 

0.6276 

0.1041 

0.0932 

0.8341 

S'*  Order 

0.6276 

0.0796 

0.0892 

0.8732 

Pressure  Interval  Studies 

It  is  seen  that  the  nonlinear  equation  provides  a 
good  fit  for  the  BDP  data  with  nine  pressures  and  two 
temperatures.  What  happens  if  fewer  pressures  are 
used  to  model  the  bum  rate  curves?  The  next  check  of 
the  nonlinear  equation  involves  five  pressures  and  two 
temperatures.  The  same  pressure  range  is  used  as  was 
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used  before,  only  this  time  .the  second,  fourth,  si.xth  and 
eighth  pressure  measurements  are  removed. 

The  reference  for  the  set  with  five  pressures  and 
two  temperatures  is  the  case  with  nine  pressures  and 
two  temperatures.  The  purpose  is  to  find  out  if  the 
results  from  the  five  pressure  case  are  similar  to  the 
results  from  the  nine  pressure  case,  thereby  reducing 
the  amount  of  strand  burner  data  necessary  for 
evaluation.  The  uncertainty  around  the  mean  is  based 
on  the  mean  bum  rate  value  from  the  nine  pressure 
case.  The  constants  from  the  nine  pressure  case  were 
used  with  the  predicted  bum  rate  values  from  the  five 
pressure  case  to  determine  the  sum  of  the  squares  of 
the  residuals.  The  standard  error  of  the  estimate  and 
the  coefficient  of  determination  values  are  based  on  the 
previous  two  statements. 


Table  5  Summary  of  the  Equation  Analysis  for  Five 
Pressures  and  Two  Temperatures  _ 


Equation 

Type 

Uncertainty 
Around  the 
Mean 

Sum  of  the 
Squares  of  the 
Residuals 

Standard 
Eftot  of  the 
Estimate 

St. 

Robert’s 

I^w 

0.5857 

0.2914 

0.204 

1*  Order  | 

0.5857 

0.2957 

0.222 

2“' Order 

0.5857 

0.0988 

0.1572 

S'"  Order 

CONCLUSIONS 


The  primary  objective  of  this  work  was  to  test  and 
verify  a  nonlinear  solid  propellant  bum  rate  equation 
to  be  used  in  the  performance  prediction  of  a  boost 
motor  of  a  DRE.  The  research  began  with  a  review  of 
ternperamre  sensitivity  theory.  Then  the  nonlinear 
expression  was  identified.  This  e.\pression  was 
analyzed  using  a  new  least-squares  multiple 
polynomial  regression  method.  The  conclusions  are 
stated  below. 

It  was  observed  that 

•  using  a  nonlinear  equation  improves  the  correlation 
to  the  model  data  up  to  the  third  order  in  all  cases. 

•  using  higher  order  equations  are  more  effective 
than  mnning  more  tests  to  improve  the  correlation 
coefficient,  and 


•  uncertainty  lowers  the  correlation  coefficients. 
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Abstract 


Fullerenes  molecules)  have  produced  wide  interest  in  many  fields  because  of  the  many  suggested 

potential  applications  including  high  energy  fuels,  polymers  and  biomaterials.  The  current  research  seeks  to  expand 
upon  the  ^damental  knowledge  base  concerning  the  ignition  and  combustion  of  fiillerene  particles  injected  into 
high  temperature  combustion  gases  within  an  oxygen  rich  environment.  The  thermal  and  chemical  environment  of 
the  combustion  gases  are  produced  by  a  laminar,  premixed  methane/air  flame  stabilized  on  a  porous  flat  flame 
burner  apparatus.  An  optical  measurement  system  is  used  to  non>intrusively  measure  1}  ignition  characteristics  and 
2)  the  total  burnout  time  under  varying  oxygen  concentrations  and  flame  temperatures. 


Introduction 

Air  breathing  rocket  technologies  having  high 
energy  density  can  provide  longer  range  with  lower 
propellant  volume.  These  capabilities  have  maintained 
a  continued  interest  in  long  range,  volume  limited 
propulsion  systems.*  Advances  in  missile  propulsion 
system  technology  have  provided  large  increases  in 
missile  performance.  In  the  ducted  rocket  engine 
(DRE)  concept  (Figure  1),  the  oxidizer  is  air  while  the 
fuel  is  pre-combusted  gas  produced  by  a  fuel-rich,  end 
burning,  Glycidyl  Azide  Polymer  (GAP)  solid  fuel-gas 
generator  (SFGG).  The  requirements  for  improved 
fuels  are  guided  by  three  primary  goals:  1)  increased 
performance  (i.e.  longer  range,  higher  velocity,  and 
thrust  control);  2)  joint  service  requirements  for 
insensitive  munitions;  and  3)  minimum  signature 
requirements. 
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The  evaluation  and  selection  of  a  high  energy  fuel 
additive  is  an  important  characteristic  when 
investigating  combustion  efficiency.  In  practice,  this 
requires  that  the  ignition  and  combustion  times  of  solid 
particulate  fuel  additives  match  residence  times  within 
the  combustor.  For  example,  the  ignition  of  carbon  is 
substantially  simpler  to  achieve  as  opposed  to  boron 
which  makes  carbon  an  attractive  high  energy 
additive.^^  This  practical  difficulty  observed  with  the 
higher  energy  density  materials  (i.e.  boron)  has  lead  to 
the  development  of  thin  magnesium^  and  GAP 
coatings^  which  alter  the  ignition  kinetics  and  reduces 
combustion  times. 

Fullerenes  (Csc/Cjo  molecules)  have  produced 
wide  interest  because  of  the  many  suggested  potential 
applications  including  super-conductors,  lubricants, 
catalysts,  high  energy  fuels,  polymers  and  biomaterials. 
Fullerenes  were  first  detected  in  carbon  vapor  produced 
by  laser  evaporation  of  graphite  in  1985.^  The 
spherical  nature  of  the  structure,  which  has  no  edge 
atoms  vulnerable  to  reaction,  was  proposed  to  explain 
the  observed  high  stability  of  certain  Carbon  clusters 
relative  to  that  of  others  at  high  temperature  and  in  the 
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presence  of  an  oxidizing  gas.  The  proposed  structure 
remained  unconfirmed  until  1990^,  when  samples  large 
enough  for  spectroscopy  were  produced  by 
vaporization  of  graphite  rods  with  resistive  heating 
under  an  inert  atmosphere.  This  production  method 
was  quickly  adopted  by  several  research  ^oups  and 
small  companies.  The  interest  rapidly  expanded  to 
include  fiillerene  based  entities,  i.e.  fiillerenes  with 
metals,  hydrogen,  oxygen,  methyl,  phenyl,  or  other 
ftmctional  groups  or  cross-links.  Physical  properties  of 
the  novel  fullerenes  are  currently  under  investigation  by 
various  researchers.^’^ 

With  this  information,  propulsion  system 
developers  have  begun  development  of  DRE 
technologies  which  can  be  deployed  as  next-generation 
propulsion  systems.  Fullerenes  may  be  the  best  solid 
fuel  additive  for  air-breathing  propulsion  systems  that 
require  minimum  signature.  As  a  pure  fuel,  carbon  has 
a  favorable  density  range  of  1.6  to  1.9  gm/cm^  and  has 
excellent  gravimetric  and  volumetric  heating  values  of 
7.8  kcal/gm  and  14kcal/cm^  respectively.  The  Ceo 
molecule  has  a  density  of  1.67  gm/cm  and  bond 
energy  of  7  eV/atom  whereas  graphite  and  diamond  are 
both  less  than  this  value.  While  the  heat  of  formation 
for  carbon  black  is  zero  and  for  diamond  it  is 
insignificant,  the  heat  of  formation  for  is  +757 
cal/gm.  The  complete  combustion  by-product  of  a 
carbon/oxygen  fiiel/oxidizer  system  is  carbon  dioxide- 
the  epitome  of  minimum  signature  with  no  solids  or 
water  vapor  to  form  a  condensation  trail.  Thus,  one  of 
the  major  goals  of  the  carbon/ducted  air  system  is  to 
achieve  high  combustion  efficiency. 

Carbon  black  agglomerates  and  fiillerene 
particles  are  investigated  as  high  energy  fuel  additives 
for  a  DRE  concept.  A  SFGG  decomposes  to  supply 
heated  fuel-rich  gases  and  carbon  particles  to  the  rocket 
combustor.  Within  this  combustor  the  solid  carbon 
particles  through  their  exothermic  heat  release  provide 
additional  thrust  to  the  rocket.  The  salient  parameters 
which  have  been  identified  in  the  literature  to  have 
significant  effects  on  combustion  rates  are  the 
following:  1)  particle  size,  2)  ambient  temperature, 
3)  oxidant  concentration,  4)  homogenous  versus 
heterogeneous  ignition  and  combustion, 
5)  agglomeration,  6)  particle  surface  morphology  and 
porosity,  7)  convective  environments,  8)  chemical 
kinetic  mechanisms,  9)  particle  surface  temperature, 
and  10)  diffusion  processes.  Issues  one  through  six  are 
planned  for  evaluation  in  the  current  facility  at  the 


Propulsion  Research  Center  at  the  University  of 
Alabama  in  Huntsville. 

Experimental  Methods 

A  flat  flame  burner  apparatus  has  been 
assembled  to  approximate  the  thermal  and  chemical 
environment  of  the  conditions  of  DRE  combustor.  The 
flat  flame  burner  method  has  been  selected  because  the 
flame  is  laminar,  steady,  easily  controlled,  and 
approaches  chemical  equilibrium  conditions  which  can 
be  computed  using  a  chemical  equilibrium  code.^’  In 
the  actual  propulsion  system,  a  GAP"based  gas 
generator  provides  heated  fuel-rich  gases  and  carbon 
particles.  The  process  is  simulated  by  injecting  carbon 
particles  or  agglomerates  into  the  premixed  laminar 
flame.  Qualitative  ignition  and  combustion  results  for 
both  carbon  black  and  fiillerene  (C«/C7o  molecules) 
agglomerates  are  reported.  The  thermal  and  chemical 
environment  of  the  flame  can  be  altered  by  varying  the 
fiiel/air  ratio  to  investigate  conditions  under  which 
ignition  and  complete  combustion  of  the  particulates 
can  be  optimued.  An  optical  imaging  system  is  used  to 
record  the  flame  radiation  of  the  burning  particles. 
Ignition  and  combustion  times  are  determined  from  the 
streak  trails  recorded  by  the  imaging  system.  Prior  to 
testing,  some  candidate  carbon  particles  were 
previously  examined  for  their  size,  surface  morphology 
and  agglomerate  behavior  using  an  Atomic  Force 
Microscope  (AFM).^ 

Flat  Flame  Burner  Apparatus 

A  schematic  of  the  laminar  flat  flame  burner 
and  associated  equipment  is  shown  in  Figure  2.  The 
flat  flame  burner  is  manufactured  by  McKenna 
Products,  Incorporated,  and  consists  of  a  6  cm  diameter 
sintered  porous  plate,  a  1.4  cm  wide  sintered  porous 
shroud  ring,  and  a  3  mm  diameter  particle  inlet  tube  in 
the  center.  The  flat  flame  burner  is  used  to  generate 
laminar  premixed  flames  with  near  uniform 
temperature  and  uniform  species  concentrations  in  the 
flame  region.  The  flame  temperature  is  known  to 
decrease  slightly  with  height  above  the  burner.^  The 
nitrogen  shroud  flow  is  used  to  prevent  the  diffusion  of 
additional  oxidizers  from  the  ambient  air  into  the  center 
of  the  flame  and  to  prevent  flame  destabilization.  The 
flow  rates  and  pressures  of  all  gas  flows  are  monitored 
by  rotometers  and  pressure  gauges  to  regulate  known 
reactant  fuel/air  ratios  and  to  ensure  steady  state 
behavior.  Water  from  standard  city  water  lines  cools 
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the  burner.  The  flat  flame  burner  apparatus  is  similar  to  system  records  ignition  and  combustion  distances  based 
that  used  by  several  researchers.^^’  ’  upon  the  exposure  time  of  the  intensifler. 


Tht  carbon  agglomerates  are  injected 
individually  into  the  premixed  flame  through  the 
centered  tube.  A  fluidized  bed  is  loaded  with  either 
carbon  black  agglomerates  (ranging  in  size  from  45  pm 
to  63  pm)  or  fullerene  particles  (mean  size  of  35  pm) 
and  100  micron  silica  spheres.  The  experimentation  is 
performed  with  a  weight  loading  ratio  of  50:1  of  silica 
to  carbon.  The  fluidized  bed  consists  of  a  vertical  glass 
tube  of  1.5  inches  in  diameter  and  18  inches  in  length. 
A  low  nitrogen  gas  flow  is  passed  through  the  bed  to 
adequately  levitate  the  carbon  particles;  however,  the 
nitrogen  flow  rate  is  insufiflcient  to  disturb  the  silica 
bed  particles.  An  additional  nitrogen  flow  is  provided 
to  enable  control  of  the  exit  carrier  gas  velocities  of  the 
particles  as  they  enter  the  flame.  The  convective 
prop)erties  of  the  ignition  attributes  and  particle  burning 
can  thereby  be  investigated  with  this  design.  The 
fluidized  bed  is  supported  between  two  end  flanges  and 
an  outer  thick-walled  plastic  tube  to  ensure  safety  and 
rigidity.  The  nitrogen  gas  enters  the  fluidized  bed 
through  a  glass  porous  plug  in  order  to  ensure  cross- 
sectional  uniformity.  The  silica  and  carbon  particle 
mixture  occupy  approximately  1.5  inches  at  the  bottom 
of  the  fluidized  bed.  At  the  base  of  the  flat  flame 
burner,  a  hypodermic  needle  of  275  microns  iimer 
diameter  has  been  positioned  inside  of  a  tee  to  further 
reduce  particle  supply  rates  to  the  flame. 

Optical  Measurement  System 

An  optical  system,  shown  by  Figure  3,  is  used 
to  accurately  measure  ignition  and  combustion  times  of 
the  particles.  The  technique  is  based  upon  the  detection 
of  flame  radiation  emitted  by  the  burning  carbon 
particles  primarily  in  the  visible  wavelength.  A 
Princeton  Instruments  Image  Intensified  CCD  (ICCD) 
camera  and  image  processing  system  equipped  with  a 
NIKON  f/1.4  105  mm  Micro  lens  is  used  to  visualize 
and  record  the  burning  particle  streaks  inside  the 
laminar  flame.  Particles  are  injected  into  the  central 
tube  and  detected  upon  ignition.  With  the  association 
of  the  1:1  ratio  image  replication  lens,  the  ICCD 
camera  clearly  captures  the  20.0  mm  by  13.2  mm 
region.  Larger  fields  of  interest  may  be  captured  by 
reducing  the  macro  capabilities  of  the  lens.  The  current 
image  intensifier  and  array  can  be  gated  to  a  minimum 
exposure  time  of  5  milliseconds.  With  upgrades,  the 
intensifier  system  can  be  gated  to  6  nanoseconds.  The 


Test  Conditions 

The  flame  conditions  simulate  the  mixing 
region  in  the  combustor  where  the  air  and  SFGG 
products  mix.  The  SFGG  effluent  has  been  estimated 
using  a  chemical  equilibrium  calculation  and  are 
summarized  in  Table  1.  The  computed  adiabatic  flame 
temperature  is  1465  K.  It  is  noted  that  the  actual 
composition  could  differ  from  that  computed  using 
chemical  equilibrium  assumptions  since  chemical 
equilibrium  is  rarely  observed  in  actual  systems.  Upon 
mixing  with  air,  the  solid  and  gaseous  combustibles 
react  and  produce  heat.  For  the  present  study  the 
conditions  of  Table  1  are  not  to  be  reproduced  but  are 
to  be  addressed  in  subsequent  research.  A  wet  gas 
flame,  similar  to  practical  conditions,  could  be 
produced  with  prembced  propane,  oxygen,  and  nitrogen 
flames.^ 


Major  Specie 

Mole  Fraction 

N2 

2230 

C(s) 

28.47 

CO 

14.95 

C02 

0.13 

CH4 

2.15 

H2 

32.10 

H2O 

0.71 

Table  1  Computed  Equilibrium  Fuel  Rich  Solid 
Propellant  Combustion  Products 

Future  efforts  will  focus  on  examining  complete  flame 
conditions  which  represent  the  regions  where  SFGG 
products  and  air  mix. 

For  the  present  research,  methane/air  flames 
have  been  used  to  evaluate  the  ignition  and  combustion 
of  the  carbon  black  and  fullerenes.  The  flow  rates  of 
methane,  air,  shroud  nitrogen,  and  fluidized  bed  carrier 
nitrogen  gas  flows  for  the  flames  reported  are 
summarized  in  Table  2.  The  flame  conditions,  as 
shown  by  Table  3,  have  been  estimated  using  chemical 
equilibrium  computations.  Actual  flame  conditions  are 
expected  to  be  different  since  chemical  equilibrium  is 
probably  not  completely  achieved;  however,  chemical 
equilibrium  is  a  good  first  estimate  of  the  representative 
species  concentrations.  It  is  assumed  that  heat  losses 
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cause  the  observed  flame  temperatures,  in  small  scale 
frame  hardware,  to  be  about  200  K  less  than  computed 
adiabatic  flame  temperatures.  The  chemical 
equilibrium  conditions  in  the  flame  ha^  been 
computed  using  a  modified  PC  version  of  tr»^  NASA 
Gordon  McBride  Chemical  equilibrium  code  known  as 
ODETRAN. 

Carbon  Particle  Samples 

The  high  density  fuel  additives  presently  being 
considered  for  SFGG  are  carbon  black  and  fullerenes  or 
buckyballs.  The  carbon  black  samples  have  been 
provided  by  R.  T.  Vanderbilt  Chemical  Company 
(Product  N991)  with  a  reported  individual  mean  size  of 
approximately  270  nm.  This  was  verified  using  Atomic 
Force  Microscopy  (AFM).^  The  individual  particles  are 
near  spherical,  non-porous,  and  approximately  300  nm, 
diameter,  in  agreement  with  the  product  data  sheet.  A 
fractal  surface  dimension  of  the  surface  represented 
was  taken  and  was  found  to  be  2.3.  A  surface  fractal 
dimension  of  2.0  corresponds  to  a  flat  plane  and  a 
surface  fractal  dimension  of  3.0  corresponds  to  a 
completely  rough  surface.  Based  upon  these 
dimensions  the  carbon  blacks  appear  to  be  smooth. 
Highly  magnified  images  of  a  single  particle  verified 
that  the  particles  are  smooth  and  non-porous. 
Agglomerates  of  these  particles  were  vibration  sieve- 
separated  into  the  size  range  of  45  to  63  pm,  loaded 
into  the  fluidized  bed  and  injected  into  the  flame. 

The  Fullerene  sample  presently  considered  is  a 
85:15  mixture  of  C^o  and  C70  molecules,  respectively. 
The  fullerenes  have  been  manufactured  and  provided 
by  the  U.S.  Army  Missile  Command  (MICOM).  These 
C6(/C7o  particles  have  been  manufactured  by  the 
application  of  a  high  intensity  voltage  and  current 
across  a  graphite  rod  in  helium.  AFM^  indicated  that 
the  C60/C70  particles  appear  to  be  rougher  on  the  surface 
and  to  have  a  higher  degree  of  scatter  in  their  particle 
size.  The  actual  size  of  the  particles  is  inconclusive 
from  the  AFM  images;  however,  the  manufacturer 
reports  a  mean  size  of  35  pm.  High  magnification  of 
the  surfaces  indicated  the  materials  were  non-porous. 

Results  and  Discussion 

Ignition  and  combustion  times  for  the  fullerene 
particles  are  more  rapid  than  the  carbon  black  samples 
considered.  The  ignition  and  combustion  of  the  carbon 
black  and  fullerene  samples  were  approximately  2.5 


millisecond,  and  1  millisecond,  respectively.  The 
fullerene  combustion  process  does  not  resemble 
classical  combustion  of  solid  spheres  such  as  carbon 
black.  The  carbon  black  combustion,  as  shown  by 
Figure  4,  exhibits  burning  characteristics  independent 
to  flame  temperature  over  the  given  test  conditions.  The 
combustion  of  the  fullerene  particle  exhibits  rapid 
deflagration  as  shown  by  Figure  6.  The  effect  of 
oxidant  concentration,  i.e.:  fuel  equivalence  ratio, 
caimot  be  determined  from  the  present  results  since 
both  the  flame  temperature  and  oxidant,  O2, 
concentration  are  not  measured  independently. 
Although  the  oxidant  concentration  in  leaner  flames  is 
higher,  the  lower  flame  temperature  delays  the  ignition 
time.  Figures  5  and  7  display  carbon  black  and 
fullerene  ignition  times,  respectively.  The  combustion 
time  for  the  carbon  black  is  rather  constant  over  the 
range  of  flame  conditions;,  on  the  other  hand,  the 
fullerene  particles  combust  in  an  rapid,  undirected 
fashion.  The  fullerene  approximated  combustion  times 
for  the  conditions  considered  are  shown  by  figure  8. 

Summary 

This  paper  describes  the  ongoing  research  into 
the  ignition  and  combustion  of  particulate  carbon, 
focusing  on  high  combustion  efficiency.  A  test 
apparatus,  including  a  flat  flame  burner  and  optical 
measurement  system,  has  been  constructed  to 
characterize  the  ignition  and  combustion  of  carbon 
black  and  fullerenes.  This  research  provides  optical 
data  of  the  reacting  particulate  carbon  forms  in  an 
oxidizer  rich  environment. 

Future  experiments  will  focus  on  laminar, 
propane/oxygen  flames,  smaller  carbon  black 
agglomerations,  and  extensive  fullerene  particle 
investigations. 
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Figure  1 .  Schematic  of  the  Ducted  Rocket  Engine. 
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Figure  2.  Schematic  ofFlat  Flame  Burner  Apparatus. 


Figure  3.  Schematic  of  Optical  Measurement  System 
for  Ignition  and  Combustion  Times. 


Fuel 

Equivalence 

Ratio 

Adiabatic 

Flame 

Temperature 

(•K) 

Methane  Flow 
Rate  (Standard 
cmVsec) 

Air  Flow  Rate 
(Standard 
cmVsec) 

Guard 

Nitrogen  Flow 
Rate  (Standard 
cmVsec) 

Nitrogen 
Fluidized  Bed 
Flow  Rate 
(Standard 
cmVsec) 

074 

1896 

21.6 

308J 

98.3 

983 

2074 

25.1 

308.3  n 

98.3 

98.3 

0.97 

2203 

28.4 

308.3 

98J 

983 

1.07 

2226 

31.5 

308J3 

98.3 

98.3 

Table  2.  Summary  of  the  Gas  Flow  Rates  for  each  Premixed  Flame. 
Standard  Flow  Rate  refers  to  Atmospheric  Pressure  and  Temperature. 
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Table  3.  Chemical  Equilibrium  Conditions  for  the  Flames  Investigated, 


Figure  4.  Streak  photograph  of  carbon  black  in  methane-air  flame. 

O  =  0.85  Particle  velocity:  617  cm/s 
Photograph  is  rotated  90®  CW  Actual  trajectory  is  vertical  upward. 
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Time  (milliseconds) 


Carbon  Black  Agglomerates 
Size:  45-66  microns 


Figure  5.  Equivalence  Ratio  vs.  ignition  and  Combustion  Times. 
Carbon  Black  Agglomerations  Size:  45-66  microns. 


Figure  6.  Streak  photograph  of  fullerene  in  methane-air  flame. 

0  =  0.85  Particle  velocity:  1342  cm/s 
Photograph  is  rotated  90°  CW.  Actual  trajectory  is  vertical  upward. 
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Time  (microseconds) 


V 


Fullerene  Particles 
Nominal  Size:  35  microns 


Equivalence  Ratio 


Figure  7.  Equivalence  Ratio  vs.  Ignition  Times. 
Fullerene  Particles  Nominal  Size:  35  microns. 
Time  Scale:  Milliseconds 


Fullerene  Particles 
Nominal  Size:  35  microns 


Equivalence  Ratio 


Figure  8.  Equivalence  Ratio  vs.  Combustion  Times. 

Fullerene  Prides  Nominal  Size:  35  microns. 
Time  Scale:  Microseconds 
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Abstract 

The  results  of  an  assessment  of  selected  data 
reduction  methods  for  performance  determination  in 
connected-pipe  ramjet  testing  from  an  uncertainty 
analysis  viewpoint  is  presented.  The  study  identifies 
and  reviews  4  distinct  methods  of  determining 
characteristic  exhaust  velocity  (C*),  efficiency  based 
on  C*,  specific  impulse  (Isp)  and  efficiency  based  on 
Isp,  and  16  distinct  methods  for  determining  thermal 
efficiency.  The  general  uncertainty  analysis  of  the 
performance  parameters  for  a  case  study  is  presented 
for  conditions  of  one  percent  imcertaindes  in  the 
input  variables  and  for  reasonable  estimates  of  the 
uncertainty  of  input  variables  for  a  liquid  fuel  ramjet 
and  a  ducted  rocket  engine.  The  study  shows  that  a 
wide  range  of  values  and  uncertainties  for  the 
performance  parameters  can  be  calculated  using  the 
different  data  reduction  methods  with  identical  input 
parameters.  Also,  the  relative  influence  of  the 
uncertainty  of  each  input  parameter  on  the  result 
uncertainty  is  presented  for  all  identified  data 
reduction  methods. 

Introduction 

Connected-pipe  testing  is  used  for  performance 
determination  and  fundamental  combustion  studies  in 
ramjet  and  scramjet  engines^*^.  A  schematic  of  a 
typical  cormected-pipe  facility  is  shown  in  Figure  1. 
In  cormected-pipe  testing,  the  air  supply  is  connected 
directly  to  the  ramjet  combustor  and,  therefore, 
connected-pipe  testing  considers  only  the  combustor 
performance  and  no  aerodynamic  or  inlet  effects. 
By  considering  only  the  ramjet  combustor,  die  air 
supply  requirement  and  equipment  necessary  for 
testing  is  minimized  making  cormected-pipe  testing 
the  most  cost-effective  method  for  evaluation  and 
development  of  ramjet  engines  prior  to  free  jet  and 
flight  testing^. 

"‘Graduate  Research  Assistant,  Student  Member  AlAA 
tEminent  Scholar  in  Propulsion  and  Professor,  Senior 
Member  AIAA 
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There  exist  various  data  reduction  methods  for 
each  of  the  performance  parameters  associated  with 
ramjet  testing  in  connected-pipe  facilities.  An 
AGARD  working  group  was  tasked  to  review  and 
report  currently  accepted  methods  for  performance 
determination  in  cormected-pipe  testing.  The 
subsequent  report^  issued  by  AGARD  included  a 
concise  listing  of  data  reduction  methods  for 
performance  determination  that  are  used  by  the 
international  technical  community.  The  various  data 
reduction  methods  yield  a  range  of  values  for  each  of 
the  performance  parameters  (such  as  thermal 
efficiency,  which  ranged  from  86%  to  96%  in  the 
sample  case  presented  in  this  paper  using  16  different 
data  reduction  methods).  The  difference  in  the 
performance  parameter  values  that  can  be  determined 
by  the  use  of  the  various  data  reduction  methods  is 
die  impetus  for  the  assessment  of  data  reduction 
methods  from  an  uncertainty  analysis  standpoint  as 
presented  in  this  paper. 

The  generic  combustor  shown  in  Figure  1 
includes  station  designations  that  have  been 
standardized^  in  order  to  simplify  the  reporting  of 
test  results.  The  station  identification  numbers  are 
used  as  subscripts  in  this  paper  to  denote  station. 
Station  2  corresponds  to  the  post  compression  inlet 
conditions.  Station  3  indicates  upstream  combustor 
conditions.  Station  4  is  used  to  indicate  the 
downstream  combustor  conditions.  Station  5 
corresponds  to  the  nozzle  throat  and  station  6  is  used 
to  indicate  the  nozzle  exit  plane. 

For  the  appropriate  simulation  of  vehicle  flight 
conditions,  the  air  must  be  supplied  at  the  stagnation 
temperatures  and  pressures  that  are  to  be  encountered 
during  flight.  In  order  to  supply  the  combustor  with 
the  high  temperature  air  necessary  to  simulate  the 
conditions  produced  by  a  supersonic  compression 
inlet,  a  vitiated  heater  is  used  to  increase  the  air 
temperature.  A  vitiated  heater  uses  combustion  of  a 
fuel  added  to  the  air  flow  to  increase  the  temperature. 
Also,  oxygen  is  added  to  the  flow  field  to  offset  the 
consumption  of  oxygen  by  the  combustion  of  the 
vitiator  fuel.  There  are  two  approaches  to  oxygen 
replenishment:  1)  make-up  oxygen  is  added  to 


account  for  the  mass  of  oxygen  consumed  by  the 
combustion  of  the  vitiator  fuel,  and  2)  make-up 
oxygen  is  added  to  preserve  the  volumetric  content  of 
oxygen  in  the  supply  air.  Both  of  these  vitiator 
methods  are  discussed  in  reference  [3].  By  adding 
the  fuel  and  makeup  oxygen,  the  composition  of  the 
oxidizer  supplied  to  the  combustion  chamber  is  no 
longer  that  of  air  because  it  includes  combustion 
products  from  the  vitiated  heater.  The  composition 
can  generally  be  considered  to  be  all  of  the 
constituents  (air,  vitiator  fuel,  and  makeup  oxygen)  at 
equilibrium  at  the  static  temperature  and  pressure  of 
the  inlets  to  the  ramjet  combustor. 

The  review  and  assessment  of  the  data 
reduction  methods  presented  in  this  paper  include  a 
description  of  the  methods  for  selected  performance 
parameters,  an  assessment  of  the  propagation  of 
uncertainty  for  each  input  parameter,  and  a  general 
uncertainty  analysis  for  each  of  the  data  reduction 
methods  for  the  selected  performance  parameters. 

Data  Reduction  Methodolosv 

The  various  data  reduction  methods  for  selected 
performance  parameters  associated  with  ramjet 
testing  in  connected-pipe  facilities  are  discussed  in 
this  section.  These  data  reduction  methods  are  taken 
from  the  AGARD  advisory  report  previously 
mentioned. 

Chemical  Equilibrium  Combustion  Codes 

The  majority  of  methods  for  performance 
determination  utilize  output  from  chemical 
equilibrium  combustion  (CEC)  codes.  The  codes, 
which  were  developed  for  classical  rocket  motor 
performance  and  species  determination,  are  based  on 
a  zero  velocity  combustion  model  since  the  Mach 
number  in  classical  rocket  motors  is  very  small.  To 
account  for  velocity  in  the  combustion  chamber  of  a 
ramjet,  the  accepted  practice  is  to  consider  stagnation 
flow  properties  as  inputs  to  the  code  instead  of  static 
flow  properties^. 

There  are  several  different  CEC  codes  in  use. 
The  two  most  frequently  used  codes  are  NASA 
CET89^’5  and  the  PEP  code^.  Previous 
investigations  have  shown  no  significant  difference 
in  the  performance  values  calculated  using  the  NASA 
CET89  and  PEP  codes^. 

Performance  Determination 

Use  of  the  different  recommended^  data 
reduction  methods  with  identical  inputs  produces 
different  values  for  the  performance  parameters.  The 


variation  in  these  determined  values  using  the  same 
input  parameters  is  one  of  the  primary  reasons  for  the 
undertaking  of  this  study  to  assess  the  different 
methods  of  performance  determination. 

The  performance  parameters  chosen  as  the 
subject  of  the  uncertainty  analysis  assessment 
presented  in  this  article  are  characteristic  exhaust 
velocity  (C*),  efficiency  based  on  C*  (ilc*)>  vacuum 
specific  impulse  (Isp),  efficiency  based  on  Isp 
and  thermal  efficiency  In  this  study,  4  distinct 

methods  are  identified  for  the  determination  of  each 
of  the  parameters  C*,  Isp,  and  qj^p  .  The 
thermal  efficiency  (q^j)  can  be  determined  from  4 
different  equations,  two  using  C*  as  input  and  two 
using  Isp  as  input,  resulting  m  16  distinct  data 
reduction  methods  for  the  determination  of  q^j. 

The  different  C*  and  Isp  methods  are  based  on 
different  methods  for  determining  the  total  pressure 
at  station  4  (pt4)-  A  diagram  illustrating  the 
relationship  between  ,  C*,  Isp,  and  q^j  is  shown 
in  Figure  2. 

The  4  methods  for  determining  Pt4  identified  in 
this  study  are 

MgthodJi 

Pt4  =  Pt4CEC 

where  Viacec  determined  by  a  CEC  code  using  the 
isentropic  flow  relations. 

Method  2: 

+  (2) 

where  P4  is  the  measured  static  pressure  at  station  4, 
M4  is  the  Mach  number  at  station  4,  and  y  is  the 
isentropic  exponent  which  relates  the  properties  of 
station  4  and  station  5  as  an  isentropic  process.  The 
value  of  y  is  obtained  through  the  relation^ 

.  '"(P4/P!) _ ^  (3, 

HP4/P5)-KT4/T5) 

where  P5  ,  T4  and  T^  can  be  found  as  output  from  the 
CEC  codes. 

The  difference  in  method  1  and  method  2  is  that 
in  method  1  the  code  uses  an  isentropic  exponent  that 
is  not  the  same  as  the  process  isentropic  exponent,  y, 
used  in  method  2.  The  use  of  the  process  y  as 
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defined  in  Equation  (3)  is  the  recommendation  of  the 
AGARD  working  group  report  since  it  exactly  relates 
the  properties  of  the  end  states  (combustor  exit  and 
nozzle  throat)  by  an  isentropic  process. 

The  process  y  defined  above  is  not  calculated 
directly  by  all  codes.  The  value  of  the  process  y  lies 
within  the  range  of  values  for  the  isentropic  exponent 
based  on  frozen  flow  and  the  isentropic  exponent 
based  on  equilibrium  flow.  For  simplicity,  in  this 
paper  y  always  refers  to  tiie  process  isentropic 
exponent  based  on  the  flow  properties  at  station  4 
and  station  5  as  defined  in  Equation  (3). 


Method  3: 


Fs 

Pt4=  — 


'  Pamb^S 


(1  +  YCD5)\5'  -  ' 

where  is  the  ambient  pressure,  A5  is  the  nozzle 
throat  area  (station  5),  CD5  is  the  nozzle  discharge 
coefficient,  and  F5  is  the  stream  thrust  as  determined 
by 

F5  =  FlC  FpL  “  Ab(pb  -  Pamb)  (^) 


where  Flc  is  the  load  cell  measurement,  FpL  is  the 
preload  on  the  load  cell,  A^,  is  the  nozzle  base  area, 
and  p^  is  the  base  pressure. 


where  ±4  is  the  total  mass  flow  rate  exiting  the 
combustor  and  is  from  one  of  Equations  (1),  (2), 
(4),  or  (6),  thus  giving  four  ways  to  determine  C*. 

The  efficiency  based  on  C*  is  determined  from 

C* 

Tic*  =7:7-  (9) 

'-'th 

where  C*  is  determined  from  Equation  (8). 

The  determination  of  Isp  is  based  on  the 
equation 


Isp  = 


Pt4As  ^  ^  1 

ihs  Vy  +  lJ 


y 

y-l 


•(l  +  yCDj)  (10) 


where  is  the  mass  flow  rate  at  the  nozzle  throat 
(assuming  that  the  mass  flow  in  the  combustor  is 
expanded  through  the  nozzle,  015  =  014 ),  and  p^^  is 
from  one  of  Equations  (1),  (2),  (4),  or  (6),  thus  giving 
four  ways  to  determine  Isp.  It  should  be  noted  that, 
when  using  Equation  (10)  for  determining  the  Isp 
based  on  the  total  pressure,  p^^  ,  determined  by  using 
Equation  (4),  the  influence  of  p^^  is  actually 
eliminated  resulting  in 

Isp  ^hllss^  (11) 

TTTs 


Method  4: 


Pt4  = 


c; 


ispth; 


^5  Pamb-^5 


CDsAs 


(6) 


where  C*  is  the  maximum  theoretical  C*  calculated 
by  the  chemical  equilibrium  code  and  Isp^,  is 
determined  by  the  code  from 


The  efficiency  based  on  Isp  is  determined  by 

(12) 


Isp 

Tlkp  =  7-^ 

Ispa, 


where  Isp^^  is  the  theoretical  maximum  value  for  Isp 
determined  by  Equation  (7)  and  Isp  is  determined 
from  Equation  (10)  or  (11). 

The  determination  of  thermal  efficiency  is 
based  on  the  equation 


(7) 

ms 

where  ihs  is  the  mass  flow  rate  at  the  nozzle  throat 
and  C5  is  the  velocity  of  exhaust  gases  at  the  nozzle 
throat.  Assuming  that  the  nozzle  is  choked,  this 
corresponds  to  the  speed  of  sound  which  is  an  output 
of  the  CEC  codes. 

The  determination  of  C*  is  based  on  die 
equation 

C*^Pt4A£D^  (8) 

m4 


t1at  = 


"^14, exp  '^12 
Tt4,th  “  Tt2 


(13) 


where  Tt4,exp  ,  the  total  temperature  at  station  4,  is 
calculated  using  experimental  measurements  and 
T^4  is  the  theoretical  total  temperature  at  station  4 
(adiabatic  flame  temperature  based  on  stagnation 
flow  properties)  as  determined  by  the  use  of  a  CEC 
code.  Direct  measurements  to  determine  the  cross- 
section  average  value  of  T^4  are  generally  not 
performed  due  to  the  difficulty  in  obtaining  good 
results^.  For  this  reason,  the  direct  measurement  of 
Tt4  exp  is  not  considered  in  this  paper.  Tt4  gxp  can  be 
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determined  from  four  equations,  two  using  C*  as 
inputs  and  two  using  Isp  as  inputs.  Therefore,  each 
equation  for  Wp  represents  four  different  data 
reduction  methods  for  as  shown  in  Figure  2. 
The  equations  for  the  determination  of  Tt4  are 


of  bias  and  precision  uncertainty  components. 
Consider  the  uncertainty  of  a  result,  r,  of  a  data 
reduction  method  with  J  input  variables  Xj ,  such  tiiat 

r  =  r(xi,X2,...,Xj)  (18) 


and 


and 


Ill 

2  Ir-i 


t4.exp  jJ 


f  /-■*2  ^ 

^exp 


\  ^4  J 


'^exp 


^t4,exp 


‘’t4,exp 


V  y 


2(y  +  l) 


M4,th 


^Ispexp"' 


\  ^4  J 


(14) 


(15) 


(16) 


and 


^t4,cxp 


risp_ 

I  Isptf, 


>2 


*  ‘^t4,th 


(17) 


where  R4  is  the  gas  constant  as  determined  by  the  use 
of  a  CEC  code. 

TTie  diagram  provided  in  Figure  2  illustrates  the 
relationship  between  the  various  data  reduction 
methods.  Each  method  for  C*  and  Isp  corresponds  to 
a  method  to  detennine  Pt4  .  The  methods  for  use 
the  four  different  methods  to  determine  pj4  in  each  of 
the  four  different  equations  for  1^4  resulting  in  16 
methods  for  The  subscripts  used  in  the 

uncertainty  analysis  section  correspond  to  Figtire  2 
and  are  provided  for  the  identification  of  specific 
data  reduction  methods  for  each  performance 
parameter. 


The  overall  uncertainty  in  r  can  be  determined  by 


Ur  = 


^vaxj; 


(19) 


where  Uv  is  the  uncertainty  in  the  input  variable  X:  , 

and  where  the  partial  derivatives  are  sometimes 
referred  to  as  “sensitivity  coefficients.”  The  interval 
r  ±  contains  the  true  (but  unknown)  value  of  r 

about  95  times  out  of  100. 

An  uncertainty  analysis  was  performed 
considering  the  equations  for  C*,  r\Q*,  Isp,  tij^p  and 
■Hat  of  Equation  (18).  The  required 

sensitivity  coefficients  were  numerically 
approximated  by  perturbing  each  input  variable,  in 
sequence,  by  1%  and  determining  the  perturbed  value 
of  the  result.  This  allowed  determination  of  Ar/Ax^ 
for  each  result  and  each  input  variable  for  the 
particular  nominal  values  in  Ais  study.  The  PEP 
code  was  used  to  perform  all  CEC  code  runs. 

The  results  of  the  study  include  the  values  of 
performance  parameters  determined  using  the 
different  data  reduction  methods,  sensitivity 
coefficients  presented  in  the  form  of  an  “uncertainty 
magnification  factor”  (UMF),  an  uncertainty  in  r 
based  on  a  1%  uncertainty  in  all  input  parameters, 
and  an  uncertainty  in  r  based  on  “reasonable” 
estimates  of  uncertainties  in  input  parameters. 

The  sensitivity  coefficients  are  presented  in  the 
form  of  an  “uncertainty  magnification  factor”  (UMF) 
defined  as 


Uncertainty  Analysis  of  Data  Reduction  Method 
In  order  to  assess  the  uncertainty  behavior 
associated  with  the  different  data  reduction  methods, 
an  uncertainty  analysis  of  a  nominal  case  was 
performed.  The  nominal  values  of  the  case  study 
used  in  this  tmcertainty  analysis  were  taken  from  the 
case  presented  for  example  calculations  in  reference 
[3]  and  are  presented  in  Table  1. 

In  this  study,  a  general  uncertainty  analysis*^ 
was  performed,  and  therefore,  the  uncertainties  of  the 
input  variables  are  not  considered  separately  in  terms 


UMF  =  ^-—  (20) 

r  Ax; 

where  r  represents  the  result  (perfonnance  parameter) 
and  Xj  represents  an  input  variable.  The  usefulness  of 
the  UMF  value  is  that  it  illustrates  the  influence  of 
the  uncertainty  in  the  input  variable  Xj  as  it 
propagates  through  the  data  reduction  method  into 
the  result.  If  the  magnification  factor  is  less  than  1, 
this  indicates  that  the  influence  of  the  imcertainty  in 
the  input  variable  diminishes  as  the  uncertainty  is 
propagated  to  the  result.  If  the  UMF  value  is  greater 
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than  1,  the  influence  of  the  uncertainty  of  the  input 
variable  is  magnified  as  it  is  propagated  to  the  result 
(a  one  percent  uncertainty  in  the  input  variable 
accounts  for  greater  than  1%  uncertainty  in  the  final 
result).  In  this  study,  all  UMF  values  are  presented 
as  positive  numbers.  Since  the  overall  uncertainty 
equation  (Equation  (19))  uses  the  square  of  the 
sensitivity  coefficient,  sign  has  no  impact  on  the 
propagation  of  uncertainty  to  the  result  for  the 
general  uncertainty  analysis  considered  in  this  case. 

The  results  of  the  uncertainty  analysis  are 
summarized  in  Tables  2-6.  Each  table  corresponds  to 
a  specific  performance  parameter  and  affords 
evaluation  of  the  different  methods  used  in  the 
determination  of  the  performance  parameter.  In  the 
upper  row  of  values  in  these  tables,  the  calculated 
values  for  the  performance  parameter  are  given 
illustrating  the  range  of  values  Aat  can  be  attained  by 
simply  using  the  different  methods.  The  UMF  values 
fill  the  body  of  the  tables  and  can  be  cross-referenced 
to  the  input  variable  and  the  performance  parameter 
method  to  which  they  apply. 

Also  included  in  the  tables  of  the  uncertainty 
analysis  are  three  uncertainty  values  for  the  result. 
The  uncertainty  value  designated  by  Ur(l%)  is  an 
uncertainty  value  for  the  result  based  on  an  assumed 
1  percent  uncertainty  in  all  of  the  input  variables. 
The  other  uncertainty  values  presented  are  based  on 
reasonable  estimates  for  the  overall  uncertainty  in  the 
input  variables  for  the  cases  of  a  liquid  fuel  ramjet 
and  a  ducted  rocket  engine  and  are  designated  by 
Ur(LFR)  and  Ur(DRE),  respectively.  The  estimated 
uncertainties  in  the  input  variables  are  shown  in 
Table  7.  The  only  difference  in  estimates  of  the 
uncertainties  of  the  input  variables  is  the  estimated 
uncertainty  in  the  enthalpy  of  formation  of  the  fuel. 
The  larger  uncertainty  estimate  for  the  ducted  rocket 
engine  case  is  due  to  the  range  of  values  reported  for 
the  enthalpies  of  formation  for  gas  generator  fuel 
formulations.  For  example, for  fuels  based  on 
glycidyl  azide  polymers  reported  values  that  vary  by 
up  to  ±50  kcal/kg. 

Figures  3  through  7  show  a  graphical  summary 
of  the  results  of  this  study.  The  uncertainties  used  to 
construct  the  uncertainty  bands  were  the  Ur(LFR) 
values  presented  in  Tables  2  through  6. 

The  case  study  presented  in  this  paper 
represents  that  of  a  ramjet  under  a  single  set  of 
nominal  operating  conditions.  The  results  of  the 
uncertainty  analysis  presented  in  this  case  do  not 
necessarily  indicate  the  behavior  of  the  uncertainty 
under  different  regimes  of  operating  conditions. 
Furthermore,  the  uncertainty  assessment  of  data 
reduction  methodologies  in  this  study  did  not 


consider  uncertainties  based  on  assumptions  in  the 
different  methods  or  in  the  internal  computations  of 
the  CEC  codes  (such  as  y,  C5 ,  R4,  etc.). 


Summary_and  Conclusions 

There  exist  numerous  data  reduction  methods 
for  performance  determination  of  ramjet  connected- 
pipe  testing.  Identified  in  this  study  are  4  distinct 
methods  for  C*,  4  distinct  methods  for  Isp,  and  16 
distinct  methods  for  thermal  efficiency  The 

results  in  Figures  3-7  indicate  that  the  methods  of 
determination  of  r\^j  produce  a  larger  range  of 
values  than  do  the  methods  for  determining  the  other 
performance  parameters  and  that  has  the  largest 
uncertainty  of  all  of  the  performance  parameters.  It 
should  be  noted  that  the  other  efficiency  performance 
parameters  ('n(;>*  and  qj^p)  are  limited  to  a  range  of 
values  that  cannot  approach  a  lower  limit  near  zero 
since  the  air  flow  in  the  combustor  always  produces 
values  of  C*  and  Isp  that  do  not  approach  zero. 
Therefore,  while  the  uncertainties  of  and  qj^p  for 
the  sample  case  in  this  study  are  lower  than  that  of 
,  the  resolution  of  reasonable  values  of  q^*  and 
qjgp  is  less  than  that  of  reasonable  values  of  r\^j . 

The  data  reduction  methods  reviewed  in  this 
study  can  be  divided  into  the  categories  of  thrust- 
based  calculations  and  non-thrust  calculations  as 
shown  in  Figure  2.  The  uncertainty  values 
determined  (subject  to  the  assumed  uncertainties  in 
the  input  variable)  showed  that  the  thrust  based 
methods  typically  provided  results  with  lower  values 
of  imcertainty  than  the  non  thrust  calculations.  This 
conclusion  is  particularly  true  for  the  q^j  values  in 
Table  6. 

The  different  uncertainty  calculations  for  the 
liquid  fuel  ramjet  and  ducted  rocket  case  illustrate  the 
insensitivity  of  the  result  uncertainty  to  the 
imcertainty  in  the  fuel  enthalpy  of  formation  for  the 
case  study  (a  result  that  is  not  intuitively  obvious). 
This  can  also  be  illustrated  by  the  UMF  values  in  the 
table  for  the  fuel  enthalpy  of  formation.  The  UMF 
values  provide  insight  to  which  parameters  will 
influence  the  overall  uncertainty  of  the  result  the 
most.  This  can  be  viewed  as  an  indication  of  where 
to  focus  efforts  to  decrease  or  minimize  the  input 
variable  uncertainty  such  that  the  result  uncertainty 
may  be  decreased. 


5 

American  Institute  of  Aeronautics  and  Astronautics 


References 

ijimnat,  Y.  M.,  “Recent  Developments  in  Ramjets, 
Ducted  Rockets  and  Scramjets,”  Prog,  of  Aerospace 
Sci.,  Vol.  27, 1990,  pp.  201-235. 

^Dunsworth,  L.  C.  and  Reed,  G.  J.,  “Ramjet  Engine 
Testing  and  Simxilation  Techniques,”  Journal  of 
Spacecraft  and  Rockets,  Vol.  16,  No.  6,  pp.  382-388, 
1979. 

^AGARD  Advisory  Report  323,  Experimental  and 
Analytical  Methods  for  the  Determination  of 
Connected-Pipe  Ramjet  and  Ducted  Rocket  Internal 
Performance,  August,  1994 


'^McBride,  B.  J.,  Computer  Program  for  Cedculation 
of  Complex  Chemical  Equilibrium  Compositions, 
Rocket  Performance,  Incident  and  Reflected  Shocks, 
and  Chapman-Jouguet  Detonations,  NASA-SP-273, 
Revision  July  15, 1986. 

^McBride,  B.  J.,  CEr89  -  Chemical  Equilibrium  with 
Transport  Properties,  Lewis  Research  Center,  1989. 

^Cruise,  D.  R.,  Theoretical  Computations  of 
Equilibrium  Compositions,  Thermodynamic 
Properties,  and  Performance  Characteristics  of 
Propellant  Systems,  Naval  Weapons  Center,  China 
Lake,  CA,  April  1979,  NWC-TM-5164. 

’Coleman,  H.  W.  and  Steele,  W.  G.,  Experimentation 
and  Uncertainty  Analysis  for  Engineers.  John  Wiley 
&  Sons,  1989. 


Table  1.  Nominal  Values  for  Case  Study 


Input  Variable 

Nominal  Value 

ihgir  (air  flow  rate) 

6.692  kg/s 

®fuei  (fuel  flow  rate) 

0.311  kg/s 

P4  (static  press,  at  station  4) 

568800  Pa 

P2  (static  press,  at  station  2) 

650200  Pa 

Phase  pressure) 

78065  Pa 

(ambient  pressure) 

101300  Pa 

Tj2  (total  temp,  at  station  2) 

606  K 

Hgjei  (fuel  enflialpy  of  fonnation) 

-482  kcal/kg 

A4  (combustor  exit  area) 

0.022698  m2 

A5  (nozzle  throat  diameter) 

0.012668  m2 

Abase  (ttozzle  base  area) 

0.004304  m2 

CD5  (nozzle  discharge  coef.) 

1.0 

Flc  measurement) 

13400N 

FpL  (preload  measurement) 

5000  N 
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Table  2.  UMF  values  for  C* 


METHOD 

12  3  4 

InputN. 

Variabl^S^  €"‘(111/5) 

1112 

1110 

1119 

1137 

mair 

0.95 

0.94 

0.95 

0.68 

rilfiiei 

0.04 

0.05 

0.03 

0.35 

P4 

0.98 

1.00 

0.01 

0.37 

P2 

0.00 

0.00 

0.00 

0.00 

Pbase 

0.00 

0.00 

0.03 

0.02 

Pamb 

0.00 

0.00 

0.18 

0.18 

Tc 

0.00 

0.00 

0.00 

0.05 

0.00 

0.00 

0.01 

0.00 

A4 

0.33 

0.00 

0.01 

0.14 

^5 

1.16 

1.44 

0.13 

0.29 

^base 

0.00 

0.00 

0.01 

0.01 

CD5 

1.00 

1.00 

0.45 

1.00 

Flc 

0.00 

0.00 

1.37 

1.37 

FpL 

0.00 

0.00 

0.51 

0.51 

Ur(l%)  (%) 

2.1 

2.2 

1.8 

2.0 

Ur(LFR)  (%) 

2.0 

2.1 

1.2 

1.5 

Ur(DRE)  (%) 

2.0 

2.1 

1.2 

1.5 

Table  3.  UMF  values  for  ti^* 


METHOD 

12  3  4 

Input'V 

VariableK^TiQ.  (%) 

95.1 

94.9 

95.6 

97.1 

0.67 

0.66 

0.67 

0.40 

^fuel 

0.57 

0.58 

0.56 

0.18 

P4 

0.96 

0.98 

0.01 

0,39 

P2 

0.00 

0.00 

0.00 

0.00 

Phase 

0.02 

0.02 

0.05 

0.04 

Pamb 

0.00 

0.00 

0.18 

0.18 

Tt2 

0.10 

0.10 

0.00 

0.00 

0.01 

0.00 

A4 

0.34 

0.02 

0.01 

0.12 

^5 

1.16 

1.44 

0.13 

Hggl 

Abase 

0.00 

0.00 

0.01 

0.01 

CD5 

■jjQjl 

0.00 

0.56 

0.00 

Flc 

0.00 

0.00 

1.37 

1.37 

FpL 

0.00 

0.00 

0.51 

0.51 

Ur(l%)  (%) 

1.8 

2.0 

1.8 

1.6 

Ur(LFR)  (%) 

1.6 

1.8 

1.2 

1.0 

Ur(DRE)  (%) 

1.6 

1.8 

1.2 

1.0 

Table  4.  UMF  values  for  Isp 


METHOD 

12  3  4 

InputN. 

VariablesSjsp{Ns/kg) 

1389 

1387 

1397 

1419 

“air 

0.94 

0.94 

0.95 

0.67 

“fuel 

0.05 

0.06 

0.04 

0.34 

P4 

0.97 

0.99 

0.00 

0.38 

P2 

0.00 

0.00 

0.00 

0.00 

Phase 

0.01 

0.01 

0.03 

0.02 

Pamb 

0.00 

0.00 

0.18 

0.18 

T.2 

0.00 

0.00 

0.00 

0.04 

^fuel 

0.01 

0.01 

0.00 

0.00 

A4 

0.33 

0.01 

0.00 

0.13 

^5 

1.16 

1.44 

0.13 

0.29 

A 

^basc 

0.00 

0.00 

0.01 

0.01 

CD5 

0.56 

0.56 

0.00 

0.56 

Flc 

0.00 

0.00 

1.37 

1.37 

FpL 

0.00 

0.00 

0.51 

0.51 

Ur(l%)  (%) 

1.9 

2.1 

1.8 

1.8 

Ur(LFR)  (%) 

1.8 

1.9 

1.1 

1.2 

Ui<DRE)  {%) 

1.8 

2.0 

1.1 

1.2 

Tables.  UMF  values  for 


METHOD 

12  3  4 

96.6 

96.4 

97.1 

98.7 

^air 

0.39 

0.39 

0.40 

0.12 

*^fuel 

0.19 

0.19 

0.18 

0.20 

P4 

0.58 

0.60 

0.39 

0.76 

P2 

0.00 

0.00 

0.00 

0.00 

Phase 

0.01 

0.01 

0.04 

0.02 

Pamb 

0.00 

0.00 

0.18 

0.18 

Ta 

0.06 

0.06 

0.06 

0.02 

Kftiel 

0.00 

0.01 

0.00 

0.00 

A4 

0.21 

0.11 

0.12 

0.25 

^5 

0.73 

1.01 

0.29 

0.72 

^base 

0.00 

0.00 

0.01 

0.01 

CDj 

0.56 

0.56 

0.00 

0.56 

Flc 

0.00 

0.00 

1.37 

1.37 

^PL 

0.00 

0.00 

0.51 

0.51 

Ur(l%)  (%) 

1.2 

1.4 

1.6 

1.9 

Ur(LFR)  (%) 

1.1 

1.4 

Ur(DRE)  (%) 

1.1 

1.3 

1.0 

1.4 
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Table  7.  Estimated  Uncertainties  of  Input  Variables 


li^ut 

Uncertainty,  U^. 

Vari^le 

Liquid  Fuel 
Ramjet 

Ducted  Rocket 
Engine 

rilair  (kg/s) 

0.06 

0.06 

mfiid  (kg/s) 

0.0016 

0.0016 

P4  (Pa) 

6900 

6900 

P2  (Pa) 

6900 

6900 

Phase  (Pa) 

1000 

1000 

Pamb  (Pa) 

500 

500 

Tc  (K) 

5 

5 

Hfuj,  (kcal/kg) 

5 

50 

A4  (in2) 

0.000136 

0.000136 

Aj  (in2) 

0.00010134 

0,00010134 

Abase  (“^) 

0.00005814 

0.00005814 

CDj 

0.01 

0.01 

Flc  (N) 

67 

67 

PpL  (N) 

25 

25 

n  Ramjet 

Generatnr  ^ 


2 


Figure  1.  Schematic  of  Connected-P^e  Testing 
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C*  (m/s) 


Non-Thrust  Calculations 


Thrust-Based  Calculations 


/ 

[C*],  Psplj  hlspli 


hAxli,! 


/ 

[C*]2  lnc‘)2  PSPI2 

I 


[C*]3  hc’ls  PSPI3  Illlspls 


I 

[^ATli,3 


[C*]4  hc»]4  Psp]4  l%p]4 


I 

hAT]i,4 


i  =  1  to  4,  representing  the  four  equations  for  7^4 

1  Equation  (14) 

2  ^  Equation  (15) 

3  ^  Equation  (16) 

4  =>  Equation  (17) 

Figure  2.  Relationship  Between  Data  Reduction  Methods 


Analysis  Results  for  C*  Methods 


Figures.  Results  for  C*  Methods 
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Isp  (N*s/kg) 
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Thermal  Efficiency  (%)  I  I  Isp  Efficiency  (%) 


Analysis  Results  for  Isp  Efficiency  Methods 


Figure  6.  Results  for  t^j^p  Methods 


100  4- 


96  + 


92  4- 


Anaiysis  Results  for  Thermal  EfTiciency  Methods 


Methods 


Figure  7.  Results 
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Abstract 

An  experimental  program  evaluating  ducted 
rocket  engine  propellant  performance  in  a  connected- 
pipe  facility  is  underway.  An  overview  of 
connected-pipe  testing  and  the  facility  used  in  this 
study  is  presented.  Sample  test  results  from  the 
current  study  are  reported.  Also,  further 
modifications  of  the  facility  to  enhance  ramjet  and 
ducted  rocket  research  and  development  programs  at 
the  U.S.  Army  Missile  Command  are  discussed. 

Nomenclature 
Initial  mass  of  gas  generator 

°^final  generator 

Po  Gas  generator  pressure 

t  Time 

IntrQdufftjpn 

This  article  contains  an  overview  of 
connected-pipe  testing,  the  connected-pipe  facility 
located  at  the  Propulsion  Directorate  of  the  U.S. 
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Army  Missile  Command  (MICOM),  Redstone 
Arsenal,  Alabama  and  sample  test  results  of  ducted 
rocket  propellant  development  studies  that  are 
underway. 

Connected-Pipe  Testing 

Connected-pipe  testing  is  used  for  performance 
determination  and  fundamental  combustion  studies  in 
ramjet  and  scramjet  engines^»^.  In  connected-pipe 
testing,  the  air  supply  is  connected  directly  to  die 
ramjet  combustor.  Therefore,  connected-pipe  testing 
considers  only  the  combustor  performance  and  no 
aerodynamic  or  inlet  effects.  By  considering  only 
the  ramjet  combustor,  the  air  supply  requirement  and 
equipment  necessary  for  testing  is  minimized  making 
connected-pipe  testing  the  most  cost-effective 
method  for  evaluation  and  development  of  ramjet 
engines  prior  to  free  jet  and  flight  testing^. 

For  the  appropriate  simulation  of  vehicle  flight 
conditions,  the  air  must  be  supplied  at  the  stagnation 
temperatures  and  pressures  that  are  to  be  encountered 
during  flight.  In  order  to  supply  the  combustor  with 
the  high  temperature  air  necessary  to  simulate  the 
conditions  produced  by  a  supersonic  compression 
inlet,  a  vitiated  heater  is  used  to  increase  the  air 
temperature.  A  vitiated  heater  uses  combustion  of  a 
fuel  added  to  the  air  flow  to  increase  the  temperature. 
Also,  oxygen  replenishment  is  added  to  the  flowfield 
to  offset  the  consumption  of  oxygen  by  the 
combustion  of  the  vitiator  fuel.  There  exist  two 
methods  of  oxygen  replenishment:  1)  make-up 
oxygen  is  added  to  account  for  the  mass  of  oxygen 
consumed  by  the  combustion  of  the  vitiator  fuel,  and 
2)  make-up  oxygen  is  added  to  preserve  the 


volumetric  content  of  oxygen  in  the  supply  air.  Both 
of  these  vitiator  methods  are  discussed  in  reference 
[3].  By  adding  the  fuel  and  makeup  oxygen,  the 
composition  of  the  oxidizer  supplied  to  the 
combustion  chamber  is  no  longer  that  of  air  because 
it  includes  combustion  products  from  the  vitiated 
heater.  The  composition  can  generally  be  considered 
to  be  all  of  the  constituents  (air,  vitiator  fuel,  and 
makeup  oxygen)  at  equilibrium  at  the  static 
temperature  and  pressure  of  the  inlets  to  the  ramjet 
combustor. 

Data  reduction  methods  for  performance 
determination  in  connected-pipe  facilities  vary  and 
are  the  subject  of  recent  studies^»^. 


upstream  of  the  vitiated  heater  in  order  to  promote 
more  complete  combustion  in  the  vitiator. 

A  high  pressure  air  ejector  is  used  to  provide 
subatmospheric  exhaust  pressure  conditions.  The 
high  pressure  air  ejector  allows  the  nozzle  base 
pressure  to  be  maintained  at  pressures  as  low  as  4 
psia. 

The  resulting  flight  simulation  envelope 
expressed  in  terms  of  altitude  and  Mach  number, 
based  on  the  temperatures  and  air  flow  rates 
achievable  by  the  facility  is  shown  in  Figure  4. 

Propellant  Comparison  Studies 


Ducted  Rockets 

Ducted  rockets  are  characterized  by  a  fuel 
supply  to  the  ramjet  combustor  from  a  fuel  rich  gas 
generator.  Air  is  supplied  through  side  inlets.  A 
schematic  of  a  generic  ducted  rocket  is  shown  in 
Figure  1.  In  a  ducted  rocket,  the  rocket  body  is 
boosted  to  supersonic  velocities  by  using  a  solid 
propeUant.  Upon  completion  of  the  boost  phase,  the 
inlets  are  blown  open  and  a  nozzle  insert  ejected 
allowing  use  of  the  boost  phase  combustion  chamber 
to  be  the  ramjet  combustor. 


MirOM^s  Connected-Pipe  Facility 

The  connected-pipe  facility  located  at  the 
MICOM  Propulsion  Directorate  has  recently 
undergone  a  comprehensive  renovation  and  is 
currently  in  use  for  the  evaluation  of  ducted  rocket 
gas  generator  fuels.  MICOM’ s  coimected-pipe 
facility  is  a  comprehensive  air  breathing  combustion 
development  tool  affording  the  ability  to  perform 
ducted  rocket  fuel  comparison  studies  and  other 
ramjet  testing  including  the  ability  to  perform 
trajectory  simulation  runs.  A  schematic  and  a 
photograph  of  the  facility  are  shown  in  Figures  2  and 
3,  respectively. 

The  facility  utilizes  a  500  cubic  feet  air  stor^e 
system  at  a  maximum  pressure  of 2500  psi  for  supply 
air  which  can  be  varied  at  flow  rates  up  to  10 
Ibm/second.  The  large  air  capacity  affords  long  run 
times  exceeding  6  minutes  in  duration,  even  at  the 
largest  air  flow  rate. 

The  vitiated  heater  uses  gaseous  hydrogen  as  a 
fuel  and  can  provide  inlet  temperatures  up  to  1500 
Rankine  at  iht  maximum  air  flow  rate  of  10 
Ibm/second.  The  replenishment  oxygen  is  added 


The  initial  test  series  for  the  evaluation  of  two 
proprietary  candidate  formulations  of  ducted  rocket 
gas  generator  fuels  has  been  completed  and  further 
testing  and  propellant  evaluation  studies  are 
underway. 

A  representative  sample  of  test  results  from  a 
ducted  rocket  test  is  provided  in  Figures  5,  6,  7,  8,  9 
and  10.  These  figures  are  from  a  test  simulating 
conditions  2  kilometers  altitude  at  a  Mach  number  of 
2.8.  The  figures  related  to  propellant  performance 
(gas  generator  pressure,  rambumer  pressure,  and 
thrust)  are  presented  without  a  scale,  but  are  included 
to  illustrate  the  stability  of  the  gas  generator  and 
rambumer  combustion  during  the  ducted 
rocket/ramjet  test  run.  The  facility  operation  is 
shown  in  the  remaining  figures  to  illustrate  the  steady 
operation  of  the  facility  during  test  conditions. 

The  gas  generator  pressure  profile  relates 
directly  to  the  mass  flow  rate  of  fuel  in  a  ducted 
rocket  by 


™fuel  ■ 


initial  ^final  (0 


(1) 


where  and  are  the  initial  and  final  mass 

of  the  gas  generator  cartridge,  Pg(t)  is  the  gas 


generator  pressure  at  time  t 


and  file  integral  |Pg(t)dt 


is  evaluated  over  the  time  interval  of  the  gas 
generator  bum.  Equation  (1)  is  valid  assuming  the 
flow  is  always  choked  and  the  characteristic  exhaust 
velocity  (C*)  fixim  the  gas  generator  remains 
constant.  Therefore,  since  the  change  in  mass  is  a 
constant  and  the  pressme-time  integral  produces  a 
constant  value,  the  flow  rate  curve  has  the  same 
profile  as  the  gas  generator  curve. 

The  performance  of  the  ducted  rocket  bum  is 
evaluated  in  the  form  of  various  performance 
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parameters  such  as  characteristic  exhaust  velocity, 
(C*),  specific  impulse  (Isp)  and  thermal  efficiency 
(rith).  Due  to  the  proprietary  nature  of  this  study,  no 
such  parameters  are  reported  in  this  paper.  A 
comprehensive  list  of  performance  parameters  and 
methods  of  deteraiination  are  included  in  references 
[3]  and  [4]. 

Future  Facility  Modifications  and  Studies 

In  order  to  support  the  current  propellant 
studies  and  additional  research  topics,  the  connected- 
pipe  facility  will  be  modified  by  the  addition  of 
methane  injection  and  optical  diagnostics  The 
methane  injection  system  will  provide  the  unique 
capability  to  inject  a  well  known  gaseous  fuel  into 
the  ramjet  combustor  to  evaluate  the  best  case 
steadiness  and  repeatability  of  the  facility.  This  will 
afford  insight  into  further  propellant  comparison 
studies  by  providing  a  facility  “calibration”  of  the 
steadiness  of  the  facility  during  well  knovm  operating 
conditions  and,  therefore,  the  combustion 
unsteadiness  due  to  gas  generator  propellant  can  be 
better  characterized.  Also,  a  comprehensive 
uncertainty  analysis^  of  future  experimental  work 
will  be  aided  by  the  facility  history  and  repeatability 
information  gained  from  methane  testing. 

Optical  access  ports  for  non-intrusive  optical 
diagnostics  to  evaluate  the  reacting  flowfield  will 
also  be  added.  The  purpose  of  these  studies  will  be 
to  observe  the  reacting  flow  mixing  that  occurs  in 
these  complex  combustors. 
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Figure  1.  Typical  Ducted  Rocket  Flight  Vehicle 
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Figure  3.  Photogr^h  of  MICOM  Connected-Pipe  Facility 
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